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Preface 


The  purpose  of  this  study  was  to  support  the  Engine 
Structural  Integrity  Program  and  the  Ret irenent-for-Cause 
maintenance  program.  Both  programs  have  been  initiated  by 
the  USAF  to  extend  the  useful  service  life  of  future  and 
present  engine  components.  Successful  implementation  of 
both  programs  requires  accurate  analytical  methods  for 
predicting  crack  growth.  I  personally  found  it  rewarding  to 
develop  a  analytical  technique  for  predicting  sustained-load 
crack  growth  after  overloads  in  engine  components. 

This  study*  however*  would  not  have  been  possible 
without  all  the  help  and  support  I  received.  I  wish  to 
express  my  sincere  thanks  and  gratitude  to  Dr.  T.  Nicholas 
AFWAL/MLLN  and  his  department  for  the  use  of  their 
facilities  and  their  time  and  efforts  in  assisting  me.  In 
particular  I  would  like  to  thank  Mr.  G.  Ahrens  and  Mr.  W. 
Goddard,  UORI*  for  their  help  setting  up  the  experimental 
test  apparatus.  I  also  thank  Mr.G.  Hartman  and  Mr.  D. 
Johnson  for  helping  me  overcome  the  difficulties  learning  a 
new  computer  system.  In  addition,  I  appreciated  Capt  K. 
Harms  explaining  his  previous  work  to  me.  I  also  greatly 
appreciated  the  overall  guidance  and  support  my  advisor 
Major  G.  K.  Haritos*  AFIT/CNY*  provided  during  this  study. 

I  especially  thank  my  wife  Victoria  who  provided 
immeasurable  support  and  encouragement  when  I  needed  it. 
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This  study  investigates  methods  of  modeling  the  effects 
of  overloads  on  high-temperature  sustained-load  crack 
growth.  In  addition  to  a  model  previously  developed  for 
this  specific  problem,  a  computer  program  developed  for 
low-temperature,  high-frequency  cyclic  load  applications  was 
evaluated.  Sustained-load  hold  times  were  converted  to 
equivalent  fatigue  cycles  to  analyze  a  load  spectrum, 
consisting  of  sustained-load  with  periodic  overloads.  The 
CRACKS  crack  growth  program  was  used  with  the  Mheeler  and 
Willenborg  models  used  to  account  for  crack  growth 
retardation  due  to  overloads. 

Predictions  were  compared  with  experimental  test  data 
generated  on  specimens  of  Inconel  718  at  650  C  with  periodic 
overloads  of  either  20  or  50  percent.  Crack  measurements 
were  made  using  a  electric  potential  system.  The 
application  of  the  electric  potential  system  to  crack  growth 
measurement  following  overloads  was  extensively  evaluated. 

It  was  concluded  that  the  system  had  to  be  recalibrated 
after  each  overload  due  to  a  sudden  advancement  in  crack 
length . 

The  retardation  models  were  found  to  require  empirical 
parameters  that  depend  upon  the  stress  intensity  level  for 

Vi  i 


each  overload  application.  Using  relationships  developed 
for  these  paraneters,  the  CRACKS  program  using  the  Mheeler 
model  was  found  to  be  capable  of  predicting  the 
t I me-to-fal 1 ure  for  sustained-loading  with  periodic 
overloads  within  20  percent  of  test  data.  The  Wlllenborg 
model  was  found  to  be  Inapplicable  to  this  problem  because 
It  depends  solely  on  stress  ratio  which  has  no  physical 
meaning  for  sustained-loading.  The  Wheeler  model,  on  the 
other  hand,  could  generally  be  applied  to  sustained-load 
crack  growth  using  equivalent  fatigue  cycles.  In 
conjunction  with  the  CRACKS  computer  program,  this  could 
provide  a  powerful  new  method  for  evaluating  crack  growth 
under  general  engine  mission  spectra  Including  the  effects 


of  overloads. 


PREDICTING  THE  EFFECTS  OF  OVERLOADS 


ON  SUSTAINED-LOAD  CRACK  GROWTH 
IN  A  HIGH-TEMPERATURE  SUPERALLOY 


I .  Introduct Ion 


Background 

Design  trends  for  modern  engines  have  emphasized 
increased  performance  with  higher  thrust/we ight  ratios, 
while  also  requiring  Improved  engine  durability  and 
maintainability.  Engine  components  must  therefore  be 
designed  to  endure  the  increasing  severity  in  operating 
conditions.  Until  recently,  the  method  used  to  predict  the 
useful  service  life  of  engine  components  was  based  on 
statistical  life  models.  This  method  was  very  conservative 
and  would  retire  an  entire  population  of  engine  disks  when 
it  was  predicted  that,  statistically,  I  in  1000  disks  would 
develop  a  0.03-inch  fat igue- induced  crack  ill.  Although  this 
method  helped  preclude  disk  failures,  significant  useful 
life  remained  in  the  other  999  disks  retired.  Estimates 
placed  this  residual  useful  life  at  greater  than  10  more 


Under  the  new  “Ret irement-for-Cause“  concept  Initiated 
by  the  USAF,  this  additional  useful  life  can  be  utilized  by 
adopting  an  inspection  criterion  applied  to  conponents  after 
a  specific  period  of  time.  If  the  inspected  components  pass 
this  criterion,  they  may  be  returned  to  service.  The 
criterion  is  based  on  fracture  mechanics  calculations  to 
determine  what  minimum  crack  size,  if  undetected,  would  grow 
to  failure  before  the  next  inspection. 

Fracture  mechanics  is  also  the  basis  for  predicting 
crack  growth  as  required  by  the  Engine  Structural  Integrity 
Program  (ENSIP)  specification  (31.  The  ENSIP  specification 
requires  a  damage  tolerant  design  approach  be  applied  to 
structural  critical  components  on  all  future  USAF  engine 
designs.  Under  this  approach,  initial  flaws  or  defects  are 
assumed  to  exist  in  the  components  when  they  enter  service. 
Analysis  of  components  along  with  verification  testing  must 
demonstrate  that  the  initial  flaws  will  not  grow  to  a 
catastrophic  size  within  the  design  lifetime  of  the 
component . 

It  is  clear  that  successful  Implementation  of  the 
Ret irement-for-Cause  program  and  the  ENSIP  design  approach 
depends  upon  technical  capability  in  two  key  areas.  The 
first  is  a  Nondestructive  Evaluation  (NDE)  procedure  used  to 
determine  the  largest  initial  flaw  size  existing  in  a 
component  after  an  inspection.  This  initial  flaw  size  is 
then  assumed  to  exist  in  all  components.  The  second  is  the 


capability  to  use  analytical  prediction  models  to  accurately 
predict  the  crack  propagation  from  the  initial  crack  size  to 
failure . 

There  are  numerous  crack  growth  rate  prediction  models 
with  varying  degrees  of  complexity.  Most  of  these  models 
have  been  developed  for  application  to  airframe  components 
under  typical  airframe  spectra  including  large  numbers  of 
cycles  with  periodic  overloads.  For  engine  applications, 
the  spectra  are  simpler,  involving  fewer  cycles  and  only 
occasional  overloads.  No  complex  interaction  models  have 
been  developed  for  engine  spectrum  loading  which  involves 
both  cyclic  and  sustalned'loading. 

Simple  crack  growth  models  usually  calculate  the  growth 
cycle-by-cycle  by  integrating  a  crack  growth  rate  equation. 
More  complex  models  used  in  airframe  analysis  include 
retardation  routines  to  account  for  the  decrease  in  growth 
rate  following  a  peak  overload  cycle.  These  models  do  not 
address  crack  growth  under  sustained  (hold-time)  loading 
which  is  present  in  a  typical  engine  spectrum.  CRACKS  (41 
which  represents  the  state-of-the-art  in  airframe  spectrum 
crack  growth  analysis  is  a  complex  prediction  program  used 
to  analytically  calculate  crack  growth  under  large  spectrums 
of  cyclic  loading.  This  program  includes  the  Wheeler  and 
Wlllenborg  models  for  predicting  retardation  effects. 


This  thesis  explores  one  aspect  of  the  complex 
crack-growth-rate  prediction  problem.  Specifically,  this 
thesis  explores  the  applicability  of  existing  crack-growth 
retardation  models,  developed  for  high-frequency, 
low-temperature  airframe  applications,  to  high-temperature 
sustained-load  crack  growth  retardation.  Procedures  will  be 
developed  to  convert  sustained-load  time  to  equivalent 
fatigue  cycles  so  that  classical  Wheeler  and  Willenborg 
retardation  models  can  be  used  in  the  CRACKS  computer 
program  to  predict  sustained-load  crack  growth  rates 
following  overloads.  The  retardation  models  will  be  applied 
to  data  obtained  from  previous  experimental  work  as  well  as 
new  experimental  proof  tests  to  verify  each  model's 


In  this  study  three  plastic  zone  retardation  models 


were  examined.  First  was  the  Wheeler  model  151  which 
reduces  the  crack  growth  rate  da/dn  to  account  for 
retardation.  Second  was  the  Wlllenborg  model  [61  which 
accounts  for  retardation  by  reducing  the  maximum  and  minimum 
stress  intensity  factors.  Also,  the  minimum  stress 
intensity  factor,  if  negative,  is  truncated  at  zero. 

Finally  was  the  Overload  model  developed  by  K.  Harms.  T. 
Weerasooriya.  and  T.  Nicholas  (71  (81  which  accounts  for 
retardation  by  reducing  the  stress  intensity  factor  K.  to  a 
lower  effective  value  9.fter  each  overload.  The  theory 

behind  each  of  these  models  is  discussed  in  the  following 
sections. 


Wheeler  Model 

The  basis  for  the  Wheeler.  Wlllenborg,  and  Overload 
models  is  that  an  overload  cycle  produces  an  extended 
plastic  zone  that  retards  crack  growth.  Thus,  in  figure  1, 
r^  represents  the  plastic  zone  due  to  sustained-loading  and 
represents  the  plastic  zone  due  to  an  applied  overload. 

As  long  as  a^  is  less  than  a^.  the  sustained-load  crack  is 
growing  in  the  overload  plastic  zone  at  a  slower  or  retarded 
rate.  The  Wheeler  model  predicts  this  retardation  effect  by 
reducing  the  crack  growth  rate  while  growing  through  the 


overload  plastic  zone.  The  crack  growth  rate  da/dn  is 
reduced  by  multiplication  with  a  retardation  parameter  C^. 
This  parameter  depends  on  the  ratio  of  the  current 
sustained-load  plastic  zone  size  to  the  previous  overload 
plastic  zone  size  raised  to  a  shaping  exponent  m.  If  the 
sustained-load  plastic  zone  grows  past  the  prior  overload 
plastic  zone,  no  retardation  is  predicted.  In  terms  of  the 
symbols  used  in  figure  1,  the  retardation  parameter  is 
defined  as, 

for  a^  <  ag  <1> 

for  aj>  Sg  (2) 

where  r  *  extent  of  current  yield  zone 
P 

a2  -  Sj  +  r  =  distance  from  crack  tip  to  elastic 
^  plastic  interface 

m  =  shaping  exponent 

The  shaping  exponent  m  is  used  to  calibrate  the  retardation 

model  with  experimental  data.  Generally,  m  has  been  found 

to  be  a  material  dependent  constant. 

Once  the  shaping  exponent  is  defined  for  the  given 

material,  the  value  of  C  is  substituted  into  equation  3. 

P 

da  da 

—  (retarded)  =  C  "  <non-retarded)  (3) 

dn  ^  dn 
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The  Wheeler  model  accounts  for  retardation  by  substituting 
da/dn  (retarded)  for  da/dn  (non^retarded)  while  the  crack  is 
growing  through  the  overload  plastic  zone.  The  numerical 
procedures  used  to  implement  the  Wheeler  model  are  contained 
in  the  CRACKS  program. 


mUsn>?9rq 

The  Willenborg  model  like  the  Wheeler  model  is  based  on 

yield  zone  analyses.  The  retardation  effect  is  accounted 

for  by  a  reduction  In  stress  intensity  factor.  The 

reduction  factor  is  calculated  by  first  finding  an 

equivalent  stress,  e  ,  that  will  produce  a  plastic  zone 

ap 

from  the  current  crack  tip  location  to  the  edge  of  the 

overload  plastic  zone.  Refering  to  figure  I,  this 

corresponds  to  a  plastic  zone  of  radius  of  r^  -  Aa.  Second, 

a  stress  reduction  factor,  r  .,  is  calculated  by 

red 

subtracting  the  currently  applied  maximum  stress  from 

max 

r 

ap 


r  .  =  •■ 
red  ap 


9 

max 


(4) 


When  the  crack  growth  amount  Aa  plus  the  current  plastic 


value  of  e  .is  set  equal  to  zero,  since  the  crack 
reel 

propagation  is  no  longer  retarded.  Third,  effective  values 
of  the  currently  applied  stresses  are  calculated  by: 


-8- 


-•j  -.j.  -.f  -- 


(5) 


”  %ax  *^red 


^%ln^eff  ~  *^mln  *^red 


If  either  effective  stress  is  less  than  zero,  it  is  set 
equal  to  zero.  Finally,  using  the  effective  ■aziaua  and 
miniaua  stress  values,  effective  stress  intensity  factors 
are  calculated.  This  produces  the  following  retardation 
relationships. 


[K  ,K  .  ]  =  [K  -K  ,  ,  K  .  -K  .1 

max  Bin  eff  sax  red  a  in  red 


where  K  .  =  K  -  K 

red  ap  nax 


for  aj  <  a2 


K  .  =  0  No  retardation  for  a,  >  a. 
red  1  “  2 


After  each  crack  growth  increaent  a  new  value  of  r  is 

ap 

calculated.  The  corresponding  new  value  of  calculated 

using  equation  <4),  Is  substituted  into  equation  (7), 
yielding  the  new  values  for  the  next  growth  increaent. 


This  aodel  uses  a  linear  cuaulatlve  daaage  concept  to 
sum  the  growth  contributions  of  a  single  overload  fatigue 
cycle  and  growth  due  to  sustained  load.  The  basis  for  the 
Overload  aodel,  like  the  Wheeler  and  Wlllenborg  models,  is 
that  an  overload  cycle  produces  an  extended  plastic  zone 
that  retards  crack  growth.  The  plastic  zone  concept  is 
Illustrated  schematically  in  figure  1.  The  sustained  load 
and  overload  fatigue  cycle  have  stress  intensity  factors 


o 


o 


< 


K  and  K  and  plastic  zone  radii  denoted  by  r  and  r 
s  n  P  P 

respectively.  The  distances  fron  the  center  of  the  crack  to 
the  edge  of  the  plastic  zones  due  to  sustained  loading  and  a 
fatigue  overload  cycle  applied  when  the  crack  length  was  a 

0 

are  labeled  a,  and  a..  For  a  crack  advancement  Aa  from  a  , 
the  crack  tip  will  be  in  an  overload  plastic  zone  until  Aa 
-  ^2  '  ^1*  Ini  this  plastic  zone  the  growth  rate  will 

be  retarded.  The  Overload  model  uses  a  reduced  value  of 
stress  Intensity  factor,  account  for  the  retarded 

growth  rate.  For  modeling  purposes,  taken  In  the 

form : 

~  ^  1  -  a  eip  <-d  Aa)  1  (8) 

where  =  effective  <reduced)  stress  Intensity  factor 

K  =  sustained-load  stress  Intensity  factor 
s 

«,  d  s  modeling  parameters 
Aa  =  incremental  crack  extension 

The  parameter  d  is  chosen  such  that  steady-state  crack 
growth  will  resume  after  the  crack  has  traversed  the 
overload  plastic  zone.  Mathematically,  it  is  desired  to 
have  approach  when  Aa  approaches  (Tp  -  This 

is  accomplished  by  letting 

d  Aa  =  w  <  9  > 

When  Aa  =  r  -  r  ,  the  resulting  value  of  K  .,/K  approaches 
p  p  ef  f 


-10- 


unity  to  within  one  percent  as  shown  in  figure  2.  The  plane 
stress  plastic  zone  sizes  for  the  overload  cycle  and 


sustained  load  are  given  by> 


C  K,  /  . 

[  K  /  r  '^  /  w 

s  y 


( 10) 


where  r 

y 

material . 
equat ion 


Is  the  uniaxial  tension  yield  stress  of  the 
Substituting  Aa  =  “r"  -  r  and  equation  (10) 

p  p 

(9)  yields  an  expression  for  A: 


into 


A  = 


K  ^  <  T  ^  -  1) 
s 


(11) 


where  the  overload  ratio  t  Is  defined  by; 


(12) 


Observing  equation  (11)  It  Is  noted  A  is  only  a  function 
of  material  properties  and  test  conditions.  A  therefore 
cannot  be  used  as  an  adjustable  parameter  to  fit 
experimental  data.  This  leaves  the  parameter  ei  given  in 
equation  (8),  to  be  adjusted  to  fit  experimental  data.  The 
value  of  «  chosen  determines  K  ..  immediately  after  an 

*  X  X 

overload  application  and  therefore  can  be  used  to  model  the 
reduced  value  of  sustained-load  crack  growth  rate.  The 
effect  of  varying  a  on  the  ratio  of  K  ../K  is  seen  in 

«  X  X  5 

figure  2.  Also  this  figure  demonstrates  how  the  Overload 


model  uses  a  decreasing  exponential  to  asymptotically 


approach  the  normal  crack  grot#th. 


In  a  prior  Investigation*  Harms  made  experimental 


measurements  of  the  delay  time  before  normal  sustained-load 


crack  growth  resumed  after  an  overload  was  applied.  Haras 


[71  noted  that  delay  times  and*  hence*  the  value  of 
changed  systematically  with  K.  Plots  of  the  delay  tines  At 


versus  K  level  for  the  20\  and  50%  overload  cases  tested  by 


Harms  are  shown  In  figures  3  and  4.  The  boundaries  of  the 


data  were  used  to  approximate  the  average  delay  time 


curve.  Using  an  iteration  process*  values  of  a 


were  chosen  with  a  specified  value  of  K  to  generate  a  ^t 


for  the  model.  Adjustments  In  a  were  made  until  At. 


calculated  and  At  agreed  resonably  well.  Harms  repeated 

avg 


this  process  over  the  full  range  of  K  values  to  define  a 


functional  relationship  for  At  »  K,  and  a. 

MM  r 


The  functionals  developed  were  expressed  in 


^  At  it  it 

non-dimensional  terms  as  «/«  and  K/K  where  a  and  K  are 


threshold  values.  K  ,  a  material  property,  was  substituted 


into  equation  <l)  for  Keff  with  /\a  =  0  and  «  =  «*.  Solving 


for  a*  yields  «*  =  1  -  <  K*  /  )  with  the  limiting  values 


being: 


1  total  crack  arrest 


0  no  retardation. 
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tlELAY  TIME  (SEC) 


Figure  3  Delay  Tines  Resulting  from  20  Percent  Overloads  at 


Various  Stress  Intensities. 


DELAY  TIME  (SEC) 


Figure  4  Delay  Times  Resulting  from  50  Percent  Overloads 


at  Various  Stress  Intensities. 


V 


( 


The  a  versus  K  curves  generated  by  Harms  for  20\  and 
50\  overloads,  normalized  to  the  threshold  values,  are  shown 
In  figure  5.  The  curves  were  fit  with  a  polynomial 
equation.  The  equation  for  the  20\  overload  case  is 

Of  =  «*  [<-.0730791E-01)  <K/K*>^  +  <0. 303086)  <K/K*>^ 

-(0.422108)  <K/K*)  +  <0. 1 17517E-01 )  1  (15) 

while  that  for  the  50%  overload  case  is 

a  =  a*  C(-.121127E-01)  (K/K*)^  +  (0. 239231E-01 )  (K/K*) 

+  (0.987133)  1.  (16) 

With  a  and  A  defined,  the  Overload  model's  retardation 
effect  was  calculated  using  equation  (8).  The  value 

obtained  from  this  equation  is  used  to  account  retardation 
while  within  the  overload  plastic  zone. 

It  is  apparent  all  three  of  the  retardation  models 
discussed  are  numerically  cumbersome.  Therefore,  computer 
programs  were  developed  to  implement  the  retardation 


calculat ions . 


Conputer  Model  PeveloPMent 


I 


In  this  study  two  cosputer  programs  were  used  to 
perform  the  numerical  calculations  required  to  predict  the 
crack  growth.  First,  a  computer  program  called  'Overload' 
was  written  to  integrate  a  creep  crack  growth  rate  (da/dt> 
equation  and  incorporate  the  Overload  retardation  model. 

The  second  program  used  was  the  CRACKS  program  which 
Integrates  a  fatigue  crack  growth  (da/dn>  equation.  To  use 
this  program  for  sustained-load  crack  growth  required 
converting  sustained  load  time  into  equivalent  fatigue 
cycles.  The  crack  growth  due  to  a  simple  spectrum  of 
sustained-load  (represented  by  equivalent  fatigue  cycles) 
with  periodic  overloads  was  analyzed  using  the  CRACKS 
program.  Both  the  Wheeler  and  Willenborg  retardation 
schemes  were  used  in  CRACKS.  After  completing  the  growth 
predictions,  equivalent  cycles  were  converted  back  to 
sustained  load  time.  The  theory  and  assumptions  associated 
with  each  model  will  be  discussed  in  detail  in  the  following 
sect i ons . 


ov<?ri94«?  Pr9flr.Aj 

This  program  was  written  to  carry  out  the  numerical 
calculations  required  to  use  the  Overload  retardation 
model.  A  listing  of  the  program  is  contained  in  appendix  3. 
The  input  data  for  this  program  are  the  initial  crack  size. 


load  history,  sustained-load  crack  growth  rate,  and  coepact 
tension  specleen  diaensions.  At  anytlae,  for  a  given  crack 
length  and  sustained-load  aapLltude,  the  stress  intensity 
factor  is  calculated  using  the  coapact  tension  solution  [91 
below : 


P  (2-t^a/w)  f<a/w) 

=  - 372 - 

b  Vw”  <l-a/w)^ 

where  K  =  Stress  Intensity  factor 

f  (a/w)>:(0. 8864-4. 64(a/w)- 13. 32(a/w)^f  14. 72(a/w> 
a  s  Crack  length 
P  3  Applied  load 
b  »  Speclaen  thickness 
w  s  Speclaen  width 


<17) 


^+5.6<a/w)^) 


Once  K  is  known,  the  crack  growth  rate  is  deterained 
froa  the  crack  growth  rate  equation  relating  da/dt  .:o  K.  A 
Modified  Sigmoidal  Equation  (MSE>  is  used  to  represent  this 
relationship.  The  MSE  model,  developed  by  General  Electric 
[101,  uses  the  various  coefficients  in  the  following 
equation  to  fit  a  sigmoidal  curve  through  the  data. 

da/dt  =  [exp<B)]  tK/K*l^  (ln<K/K*)l®  (ln<K  /K)l°  <18> 

c  — 

where  da/dt  =  Crack  growth  rate 

K  =  Current  stress  intensity  value 

K  =  Threshold  stress  intensity  value 

K  =  Critical  Stress  Intensity  value 
c 

B,P,Q,D  =  Fitting  parameters  of  the  curve 
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In  this  equation,  K  and  K  are  in  units  of  (MPa  ■  )  and 

c 

da/dt  is  calculated  in  units  of  (n/sec).  The  remaining 
constants  are  non-dimensional.  Harms  13]  determined  the 
sigmoidal  coefficients  by  fitting  the  data  from  his  constant 
sustained-load  baseline  test.  The  test  data,  along  with  the 
best-fit  MSE  coefficients,  are  shown  in  figure  6. 

With  the  relationship  between  a,  K  and  da/dt  known,  it 
is  possible  to  find  the  time  it  takes  to  grow  from  an 
initial  crack  size  of  a^  to  final  crack  length  a^  using: 


A  t  = 


da 


(19) 


This  Integration  is  carried  out  numerically  by  dividing  the 
region  of  crack  growth  <a^  "  ^i^  into  a  finite  number  of 
Increments  Aa  and  summing  the  time  to  grow  each  Increment. 
For  each  small  Increment,  Aa,  the  value  of  K  is  calculated 
for  the  end  points  and  for  any  intermediate  crack  lengths 
from  equation  (17).  Similarly,  da/dt  is  determined  for  the 
end  points  and  for  discrete  points  in  between  using  equation 
(18).  The  numerical  Integration  of  each  Increment  Aa  is 
performed  via  Simpson's  rule  fill.  The  total  time  to  grow 
from  a^  to  a^  is  obtained  by  summing  the  tine  to  grow  each 
Aa  Interval.  Since  Simpson's  is  rule  only  a  numerical 
approximation  of  the  exact  integral  of  the  function,  care 
was  taken  to  ensure  that  proper  accuracy  was  carried  through 


the  calculations. 


The  accuracy  of  the  Integration  was  controlled  In  two 


ways.  First,  the  interval  of  total  crack  length  Ca^  “ 
was  divided  into  a  nuaber  of  increments  (  Aa) .  The  size  of 
each  increaent  was  based  on  the  Initial  stress  intensity 
level  at  the  beginning  of  the  Interval  and  was  determined  as 
follows.  Values  of  (da/dt)^^^  were  calculated  at  various 
K-levels.  These  values  were  multiplied  by  100  seconds  to 
obtain  the  crack  length  interval  ( Aa^^^)  which  produces  100 
seconds  of  sustained  load  growth.  Linear  functions  were 
used  to  represent  Aa^^^  as  a  function  of  K,  in  a  manner 
that  insured  the  growth  tine  interval  was  always  less  than 
100  seconds.  The  effect  of  this  first  step  was  to  use  small 
Aa  Intervals  at  lower  K  values  and  larger  Aa  intervals  as 
K  increased,  thus  limiting  the  total  tine  required  for  the 
crack  to  grow  through  any  Interval  to  less  than  100  seconds, 
for  any  value  of  crack  length  or  K. 

Each  Aa  interval  was  then  numerically  integrated  using 
Simpson’s  rule.  The  numerical  integration  errors  were 
minimized  by  testing  the  Simpson's  integration  subroutine 
for  convergence.  This  was  accomplished  by  integrating  each 
crack  length  interval  at  least  twice.  The  program  started 
by  using  two  subintervals  within  each  Aa  Interval.  This 
interval  was  then  integrated  again  doubling  the  number  of 
subintervals.  The  difference  of  the  values  of  the  time 


calculated  to  grow  through  the  Interval  using  these  two 
integrations  was  compared  with  a  convergence  value;  this 


value  was  calculated  by  multiplying  a  tolerance*  set  at 
0.0001,  by  the  time  calculated  to  grow  through  the  crack 
interval.  Since  this  time  Increment  was  always  near  100 
seconds,  the  convergence  value  corresponded  to  approximately 
0.01  seconds  difference  between  the  calculated  times  using 
the  two  successive  Integrations.  If  the  difference  exceeded 
0.01  seconds,  another  integration  was  performed,  with  the 
number  of  subintervals  being  doubled  again.  This  procedure 
was  repeated  until  the  convergence  criterion  was  satisfied. 
Assuming  that  each  Aa  Interval  resulted  in  an  error  of  0.01 
seconds,  the  cumulative  error  for  the  total  time  predictions 
is  no  greater  than  10  seconds  based  on  a  maximum  number  of 
1000  intervals  corresponding  to  a  total  test  time  of 

5 

approximately  10  seconds. 

The  numerical  integration  scheme  was  applied  to  a 
loading  history  which  involved  a  constant  sustained-load 
with  periodic  overloads.  The  input  load  history  for  each 
specimen  specified  the  times  at  which  overloads  were  applied 
and  the  percentage  of  overload.  An  example  of  how  the  load 
history  was  entered  into  the  program  is  shown  in  figure  7. 

The  program  Integrates  repeated  Aa  intervals,  equal  to 
approximately  100  seconds  of  sustained  load  time,  until  the 
total  time  exceeded  the  time  when  the  next  overload  was 
applied.  At  this  point,  the  time  when  the  last  ^a  interval 
started  and  the  time  it  took  to  grow  the  last  Aa  Interval 
were  known.  Using  a  linear  interpolation,  the  length  of  the 


last  Aa  Interval  was  reduced  so  that  the  total  growth  tl*e 
equaled  the  time  when  the  overload  was  applied.  At  this 
point*  the  retardation  effect  due  to  the  overload  was  added 
in  the  model.  This  was  accomplished  by  calculating  a 
reduced  stress  intensity  factor*  defined  by  equation 

(8).  In  this  equation*  the  modeling  parameter  4  was  related 
to  the  overload  plastic  zone  size*  while  the  parameter  a  was 
defined  by  fitting  experimental  data.  During  the  retarded 
growth  was  substituted  for  K  in  equation  (18>  until  the 

crack  grew  through  the  overload  plastic  zone. 


Time 

%  Overload 

0 

0  \ 

T1 

20  \ 

T2 

20  % 

T3 

20  % 

Figure  7  Overload  Program  Load  History  Input  Example. 
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Harms  [71  noted  that  each  time  an  overload  cycle  was 

applied  an  apparent  jump  in  crack  length  occurred.  This 

same  jump  phenomenon  was  also  noted  by  Larsen  and  Nicholas 

(121  in  their  study  of  crack-growth  transients  at  elevated 

temperature.  The  amount  of  crack  jump  seemed  to  correlate 

with  the  K  level  at  overload  application.  However,  Haras's 

data  for  jump  versus  K  level,  shown  in  figure  8,  contained  a 

large  amount  of  scatter  making  it  difficult  to  fit  a 

functional  relationship.  Haras  therefore  assumed  a  constant 

value  of  0.381  mm  jump  in  crack  length  at  each  overload.  In 

an  attempt  to  Improve  the  estimation  of  the  jump  function, 

several  other  functional  relationships  were  Investigated. 

First,  a  function  relating  the  amount  of  jump  to  the  K 

level  at  overload  application,  normalized  to  the  threshold 
*  1/2 

K  value  of  24  MPa  m  ,  was  tried.  This  function  takes  the 
form  of 

Jump  =  0.2  *  <K/K*)  mm  (20) 

and  is  labeled  curve  1  in  figure  8.  The  second  function  used 
the  Log  of  the  K  level  at  overload  application,  normalized 
to  the  threshold  K*  value  of  24  MPa  m*^^.  The  resulting 
function  takes  the  form  of 

Jump  =1.25  *  Log(K/K*)  mm  (21) 

and  is  labeled  curve  2  in  figure  8.  Finally,  a  linear 
function  with  an  initial  jump  of  0.381  mm  at  threshold 


lacreaalng  to  .508  mm  Jump  near  the  critical  stress 
intensity  was  tried.  This  function  is  labeled  curve  3  in 
figure  8.  Also  shown  in  this  figure  is  the  constant  0.381  mm 
Jump  used  by  Harms,  labeled  curve  4. 

The  proof  test,  conducted  by  Harms,  was  analyzed  using 
each  jump  function  to  predict  the  increment  in  crack  length 
caused  by  the  overload  cycle.  The  resulting  predictions  are 
shown  in  figure  9.  Curves  1  and  2  correspond  to  jump 
functions  in  equations  20  and  21*  respectively  these 
functions  underestimated  the  jump  at  low  K  values.  The 
resulting  predictions  had  significant  delay  times.  Curve  3 
corresponded  to  the  linearly  Increasing  Jump  function,  and 
was  the  most  exact  prediction  of  the  total  time  to  failure 
for  this  test.  The  constant  jump  of  0.381  mm.  labeled  curve 
4.  predicted  the  total  time  to  failure  to  within  4  percent. 
Since  the  constant  jump  used  by  Harms  gave  predictions 
within  normal  test  scatter,  it  was  decided  to  use  a  constant 
jump  of  0.381mm  for  all  further  calculations.  This  also 
allowed  a  direct  comparison  of  Harms's  work  with  the  other 
retardation  models. 


CRACK  JUMP  (MM) 


Figure  8  Crack  Juaps  froa  Overload  Cycles  at 


Various  K  Levels 


TIME  (SEC) 


Figure  9  Overload  Model  Predictions  to  Proof  Test 
64-507  Using  Various  Jump  Functions. 
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CRACKS  Program 

This  prograM  was  developed  by  R.  M.  Engle  [41  to 
predict  crack  growth  In  airframe  applications  where  high 
frequency,  low  temperature  spectrum  loading  occurs.  The 
objective  here  was  to  modify  the  original  program  so  that  It 
could  be  used  to  analyze  high  temperature  sustained  loading 
with  periodic  overloads.  A  technique  was  developed  for 
converting  sustained  loading  to  equivalent  fatigue  cycles. 
This  technique  and  all  required  programing  changes  made 
within  CRACKS  are  described  next. 

The  process  of  converting  sustained  load  creep  crack 

growth  to  equivalent  fatigue  cycles  was  achieved  by  setting 

da/dn  equal  to  da/dt  at  equivalent  AK  and  K  stress 

max 

intensity  levels.  A  one-to-one  correspondence  would  equate 
one  fatigue  cycle  to  one  second  of  sustained  load.  Since 
CRACKS  can  only  use  one  crack  growth  rate  equation,  a  method 
of  representing  both  the  rate  of  growth  due  to  sustained 
loading  and  overload  fatigue  cycles  with  the  same  equation 
was  needed.  The  growth  rate  due  to  an  overload  cycle  was 
approximated  using  previously  generated  test  data  [131  for 
Inconel  718  at  650  C,  with  an  R  ratio  of  0.1  and  frequency 
of  0,01  Hz.  This  frequency  is  approximately  that  of  the 
single  overload  cycle  in  the  experimental  part  of  the 
investigation.  This  data  is  labeled  83266G  test  data  and 
shown  in  figure  10.  Also  shown  in  the  figure  is  the  baseline 


Figure  10  MSB  Curve  Shift  Using  the  B  Parameter  and  R  Ratio 


sustained-load  signoldal  curve,  developed  by  Hams,  and 
labeled  curve  1.  The  sustained-load  growth  rate  da/dt 
(«/sec>  given  by  curve  1  equals  the  cyclic  growth  rate  da/dn 
(a/cycle).  Thus,  one  second  of  sustained-load  Is  equal  to 
one  fatigue  cycle.  In  order  to  have  the  saae  curve 
represent  both  the  sustained  loading  and  overload  cycles  a 
vertical  and  horizontal  shift  of  curve  1  was  needed. 

Referlng  to  the  sustained  load  sigaoidal  equation  below 

da/dt=[exp<B) l*(K/K*l^*tln<K/K*) )®*(ln<K  /K)l°  (22> 

c 

the  parameter  B  can  be  used  to  absorb  a  rate  multiplication 
constant  Into  the  equation.  Using  a  trial  and  error 
procedure,  a  factor  of  20  was  found  to  vertically  shift  the 
curve  to  a  position  where  an  additional  horizontal  shift 
placed  the  curve  on  the  crack  growth  rate  data  for  overload 
cycles.  The  new  B  value  was  found  by  solving  the  following 
equation  for  B 

exp  B  =  20  exp  (23) 

The  new  value  for  B  and  the  associated  new  sigmoidal  curve 
is  labeled  curve  2  in  figure  10.  The  cyclic  growth  da/dn 
(m/cycle)  now  equals  20  times  the  sustained-load  growth  rate 
da/dt  (m/sec).  Thus,  20  seconds  of  sustained-load  is  equal 
to  one  fatigue  cycle.  The  horizontal  shift  was  accomplished 
using  an  R  ratio  shift.  The  cyclic  sigmoidal  equation  shown 


da/dn=texp<B)  l*(^K/^K*l*^*t  lii<^K//lK*)  ]®*t  ItiC^K^/^K)  <24) 

can  be  Modified  for  R  ratio  effects  by  replacing  AK  with 

Ak*<1  -  R).  Ak*  with  AK**<1  -  R>»  and  A^  with  A*  *<  I  “ 

c  o 

R).  Again*  using  trial  and  error, It  was  found  that  an  R 
ratio  of  0.4  produced  the  desired  horizontal  shift.  The 
resulting  slgaoldal  equation  with  fixed  a  value  for  Ak*  = 
Ak**<1  -  0.4)  and  AK  *  A*^  *<  1  ”  0.4)  Is  labeled  curve  3 
in  figure  10.  Both  the  sustained  load  and  overload  cycle  can 
be  modeled  by  varying  the  R  ratio  to  represent  the  different 
cycles.  The  sustained  loading  Is  modeled  using  I  cycle 
equals  20  seconds  of  sustained  loading  at  an  R  ratio  of  0.4. 

The  overload  cycle  Is  modeled  as  1  cycle  equals  I  overload 
cycle  at  an  R  ratio  of  0.0.  A  graphical  representation  of 
these  cycles  is  shown  in  figure  11.  Since  AK  *  AK*(  1  -  R) 
and  R  -  0.4  for  the  sustained  loading,  it  is  easily  seen 
that  the  R  ratio  dependence  cancels  out  of  equation  (24)  for 
sustained  loading.  The  sustained  load  crack  growth  rate 
then  In  essence  is  defined  by  curve  2  In  figure  10. 

Remembering  that  the  threshold  and  critical  stress  intensity 
factors  are  fixed  at  values  corresponding  to  an  R  ratio  of 
0.4  and  the  overload  cycle  is  modeled  with  an  R  ratio  of 
0.0,  the  overload  cycle  crack  growth  rate  Is  calculated 
using  curve  3  in  figure  10  with  Ak  =  Ak»<  1  -  0.0)  or  just 


iMplementation  of  the  method  developed  to  change 
sustained  loading  into  equivalent  fatigue  cycles  required 
several  changes  to  the  CRACKS  program.  The  source  listing 
of  the  program  used  in  this  study  is  contained  in  appendix 
3.  Before  describing  the  detailed  changes  within  the 
program,  a  brief  description  of  the  overall  program  will  be 
presented^ 

The  CRACKS  program  consists  of  twenty  two  routines  of 
which  sixteen  are  basic  to  crack  growth  calculations  and  six 
others  are  used  to  implement  the  retardation  models. 

Detailed  descriptions  of  each  routine  is  contained  in  the 
CRACKS  manual  [41.  The  overall  supervisory  routine  CRACKS4 
calls  each  of  the  other  subroutines  as  needed  during  the 
calculations.  Only  the  major  changes  made  within  each 
subroutine  will  be  described,  although  each  change  usually 
required  changing  other  subroutines  that  used  the  same 
common  blocks. 

The  first  change  was  to  include  the  sigmoidal  crack 
growth  rate  equation  in  the  RATE  subroutine.  This  required 
changing  the  input  and  output  subroutines  to  read  and  write 
the  sigmoidal  coefficients.  Also,  equation  (24)  was 
programed  in  the  RATE  subroutine  for  use  whenever  the  crack 
growth  rate  was  required  by  the  main  CRACKS4  routine. 

The  second  change  was  to  include  the  ASTH  compact 
tension  stress  intensity  solution,  given  by  equation  (17). 
in  the  BETA  subroutine.  The  CRACKS  program  calculates  all  K 


values  using  the  equation 


K  =  r  Y  Via  (25) 

where  K  -  Stress  Intensity  level 
r  -  Stress  or  load 
Y  s  Variable  defining  case  solution 
a  =  Current  crack  length 

The  variable  Y  Is  used  to  define  the  case  solution  for  the 
type  geometry  being  analyzed.  Y  was  set  equal  to  equation 
(17)  divided  by  rVw~a'  .  Whenever  the  stress  Intensity  value 
for  a  given  crack  length  was  needed*  the  subroutine  K  was 
called.  This  subroutine  called  the  BETA  subroutine  where  a 
new  value  for  Y  was  calculated  for  substitution  Into 
equation  (25). 

The  final  change  made  In  CRACKS  affected  the  numerical 
integration  method  used  In  the  program.  CRACKS  was  written 
to  handle  a  very  large  number  of  fatigue  cycles.  Crack 
growth  due  to  large  spectra  of  cycles  are  calculated  using  a 
linear  approximation  technique  In  order  to  save  computation 
time.  The  basis  for  the  approximation  is  the  assumption 
that  the  damage  parameters  remain  constant  over  some  small 
increment  of  crack  growth  Aa .  Engle  1141  found  the  linear 
approximation  in  CRACKS  to  be  an  excellent  balance  between 
accuracy  and  computational  efficiency  for  very  large 
spectra.  However,  when  smaller  constant  amplitude  spectra 
were  analyzed,  the  program  was  more  efficient  using  a 


cycle-by-cycle  Runge-Kutta  nuaerical  integration  nethod 
[141.  During  the  check-out  phase  of  the  CRACKS  prograa#  it 
was  found  that  the  linear  approxlaat ion  did  not  provide  the 
saae  accuracy  as  the  Overload  program  for  constant  sustained 
loading.  Therefore#  the  program  was  changed  to  eliminate 
the  linear  approximation  and  use  the  Runge-Kutta 
cycle-by-cycle  integration  method.  With  this  change  in 
place  both  the  Overload  and  CRACKS  programs  predicted 
exactly  the  same  results  for  constant  sustained-load  growth 
containing  no  overloads. 

Results  from  the  CRACKS  program  were  obtained  in  terns 
of  equivalent  fatigue  cycles  of  growth.  This  required  a 
separate  Fortran  program  to  convert  cycles  back  into 
sustained  load  time  by  multiplying  the  equivalent  cycles  by 
20  seconds  per  cycle. 

Both  the  Overload  and  CRACKS  programs  were  now  capable 
of  predicting  sustained  load  crack  growth  with  periodic 
overloads.  The  retardation  effect  overload  cycles  produced 
was  estimated  using  the  Overload  model  in  the  Overload 
program  and  the  Wheeler  and  Willenborg  models  in  the  CRACKS 
program.  The  next  step  in  developing  the  programs  was  to 
compare  the  analytical  predictions  of  each  model  to 
experimental  test  data. 


Each  retardation  aodel  was  applied  to  typical  test 

segments  using  the  computer  programs  described  in  section 

III  to  predict  how  sustained  load  crack  growth  was  affected 

by  overload  cycles.  Six  segments  of  crack  growth  data  with 

average  delay  tines  were  selected  from  Harms's  work  for 

analysis.  The  segments  included  both  20  %  and  50  % 

overloads  applied  at  low  (30  MPa  medium  (40  MPa  m^^^> 
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and  high  (50  MPa  n  )  stress  Intensity  levels.  The  test 
data  for  each  segment,  shown  in  figures  12  through  17,  start 
at  an  initial  crack  length  which  already  includes  the  crack 
length  Jump  produced  by  the  overload  cycle.  Crack  growth 
predictions  using  the  Overload,  Wheeler,  and  Willenborg 
retardation  models  were  computed  for  each  segment  to  provide 
a  comparison  between  the  predictive  capabilities  of  each 
model.  Several  observations  were  made  on  the  flexibility  of 
each  retardation  model  to  predict  the  test  data.  Discussion 
of  each  of  the  observations  and  their  effects  follows. 

The  Overload  model's  flexibility  to  match  test  data  is 
contained  in  the  a  parameter  in  equation  (8).  This  parameter 
controls  the  delay  time  before  normal  crack  growth  resumes 
after  an  overload.  Harms  used  his  test  data  to  develop 
expressions  for  a  as  a  function  of  stress  intensity  level 
for  20  %  and  50  %  overload  ratios.  These  expressions,  given 
by  equations  (15)  and  (16),  were  used  in  the  Overload 


model.  Thus  all  the  model  variables  were  predefined  before 
the  test  segment  were  analyzed.  The  predictions  for  each  of 
the  six  segments  are  shown  in  figures  12  through  17. 

The  Wheeler  model  had  flexibility  to  fit  test  data  by 
changing  the  shaping  exponent  a.  The  shaping  exponent  m  is 
an  empirical  parameter  dependent  upon  material  and  stress 
history  [151.  The  value  of  m  for  fatigue  cycling  generally 
ranges  from  1.0  to  3.S  .  During  analysis  of  the  six  test 
segments  the  shaping  exponent  m  was  treated  as  a  variable  to 
fit  the  Wheeler  model  to  the  test  data.  Figure  18  shows  how 
the  best  fit  value  of  m  was  found  by  trial  and  error  for  a 
typical  test  segment.  A  constant  value  of  m  equal  to  6.0 
accurately  predicted  the  20  \  ov»‘rload  segments.  The  best 
fit  value  of  m  for  the  50  %  overload  segments  was  related  to 
the  stress  intensity  level  at  overload  application  and 
varied  between  6.0  and  3.5  .  This  relationship  between  m  and 
K  is  shown  in  figure  19.  In  general,  the  shaping  exponent 
decreases  as  K  increases  at  higher  overload  ratios.  A 
similar  trend  was  seen  in  the  a  parameter  which  is  used  to 
fit  experimental  data  in  the  Overload  model.  Accurate 
predictions  of  retardation  at  higher  overload  ratios  depends 
upon  relating  a  or  m  to  the  stress  intensity  at  overload 
application.  At  lower  overload  ratios  a  was  still  related 
to  the  stress  intensity  level  but  a  constant  value  of  m  was 
found  to  adequately  predict  the  retardation  affect.  The 
resulting  Wheeler  predictions  using  the  best  fit  values  of  m 


are  shown  In  figures  12  through  17. 


The  Willenborg  model  does  not  have  a  parameter,  like 
the  Overload  and  Wheeler  models,  for  use  In  fitting  test 
data.  Instead  the  model  accounts  for  retardation  by 
reducing  the  equivalent  sustained  load  fatigue  cycle  of 


Ak 

with  R  =  K 


K  -  K  . 
max  min 


.  /K  =  .4 
min  max 


(26) 


to  an  effective  value  calculated  by  substituting  IK  , 

Adt  X 

K  ,  ]  from  equation  (7)  Into  equation  (26)  to  get 
min  eff  ^  -a 


Ak  =  cAk  1  -  cAk  ,  1 

^  eff  max  eff  “  min  e 


ff 


®eff  ^  ^*^mlnUff  '  ^’^max^eff 


(27) 


Application  of  the  Willenborg  model  to  the  test  segments 

showed  that  <lld  not  equal  the  equivalent  fatigue  cycle 

R  ratio  of  0.4.  This  Is  due  to  the  truncation  of  the  minimum 

K  ,,  value  at  zero  and  the  reduction  of  both  K  and  K  , 
eff  max  min 

by  3is  defined  in  equation  (7).  Therefore,  the  model 

would  not  predict  the  same  retardation  effect  If  the 
equivalent  sustained  load  cycles  were  modeled  at  different  R 
ratios.  The  resulting  predictions  using  the  Willenborg 
model  with  the  equivalent  fatigue  cycle  modeled  with  an  R 
ratio  of  0.4  are  shown  in  figures  12  through  17.  Since  the 
Willenborg  retardation  model  was  dependent  on  the  modeling 
technique  used  to  represent  the  sustained  load  time  It  was 
eliminated  from  further  consideration  in  this  study. 
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13  Crack  Length  versus  Time  for 
20  Percent  Overload  at  Medium  K. 
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Figure  14  CracK  Length  versus  Time  for 
20  Percent  Overload  at  High  K. 
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Figure  18  Best  Fit  Shaping  Exponent  (•> 
for  20  Percent  Overload  at  Low  K. 
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Figure  19  Best  Fit  Shaping  Exponent  (■>  versus  K. 
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Several  proof  tests  were  conducted  to  provide 


additional  test  data  to  verify  the  prediction  capability  of 
the  retardation  models.  A  description  of  the  test  apparatus 
and  crack  measurement  procedures  used  during  testing 
follows. 


Ifta.t-AppAra,tu8 

The  experimental  data  for  the  proof  tests  was  gathered 
using  a  semi-automated  creep  test  system  employing  electric 
potential  drop  and  optical  readings  to  monitor  crack 
length.  The  Air  Force  Wright  Aeronautical  Laboratories* 
Materials  Laboratory*  Wr ight-Patterson  Air  Force  Base*  Oh. 
provided  the  facilities  and  equipment  to  conduct  the 
testing.  The  test  system  schematic  is  shown  in  figure  20. 
The  test  setup  included  the  following  components: 

1.  12*000  lb  Swedish  creep  test  frame 

2.  Resistance  heated  furnace  with  power  controlers 

3.  Tektronix  4051  microcomputer 

4.  Daytronic  9000  signal  processor 

5.  Hewlett-Packard  3478A  IEEE-488  digital  voltmeter 

6.  Two  Gaertner  traveling  microscopes 

7.  Current  supply  source 


The  12*000  lb. -capacity  Swedish  creep  frame  was  used  to 


load  the  specimens.  The  frame  is  constructed  using  a  lever 
and  fulcrum  principle.  The  weights  were  suspended  at  the 
end  of  a  20  to  1  lever  arm.  The  other  end  reacts  to  this 
mechanical  advantage  with  a  load  line  containing  the 
specimen.  Dynamic  loading  of  the  specimen  is  avoided  by  the 
use  of  a  hydraulic  ram  which  can  support  a  fraction  or  all 
of  the  suspended  weight.  As  pressure  is  added  or  removed 
from  the  ram,  it  removes  or  adds  the  load  to  the  specimen  in 
a  smooth  manner.  Periodic  overloads  were  applied  by 
unloading  the  specimen,  adding  the  calculated  overload 
weight  to  the  end  of  the  lever  arm.  and  reloading  the 
specimen.  Removing  overloads  was  done  using  the  same 
procedure  except  for  removing  the  weights. 

The  compact  tension  specimens  were  mounted  in  the  load 
line  using  load  rods  with  Inconel  718  clevises  and  holding 
pins.  An  electric  potential  technique,  to  be  discussed 
below,  was  used  to  measure  crack  length.  To  accommodate 
this  system,  the  specimen  was  electrically  isolated  from  the 
creep  frame  by  using  insulated  sleeves  and  pins  to  connect 
the  load  transfer  rods.  A  load  cell  to  measure  the  applied 
load  to  the  specimen  was  included  in  the  load  line.  The 
load  cell  readings  were  used  to  ensure  the  hydraulic  ram 
released  the  entire  load  to  the  specimen.  The  specimen's 
configuration  and  nominal  dimensions  are  shown  in  figure  21. 
In  addition,  the  individual  test  specimen's  dimensions  and 
the  applied  loads  are  listed  in  table  1. 


Speciaen 

Width  (aa) 

Thickness  (aa) 

Load  (KN) 

84-502 

10.381 

39.969 

10.008 

12.188 

84-503 

10.665 

40 . 008 

10.033 

12.023 

84-504 

10.267 

40.018 

10.008 

12.677 

Table  1  Speciaen  Diaensions  and  Applied  Load. 

The  two  piece  resistance  furnace  was  aounted  on  the 
creep  fraae  as  shown  in  figure  22.  The  oven  was  closed 
around  the  speciaen  and  sealed  with  a  flaaeproof  wadding. 

The  wadding  aaterial  also  served  to  Insulate  the  leads  for 
the  electric  potential  systea  froa  the  oven  and  fraae.  The 
furnace  was  constructed  with  four  independently  controlled 
power  zones.  The  power  to  each  zone  is  controlled  by  a 
t iae-proport loning  West  controller.  The  controller  used  a 
theraostat  feedback  loop  to  hold  the  oven's  teaperature 
constant.  Two  K-type,  chroael-aluael »  theraocouples  were 
spot  welded  to  the  speciaens  and  are  shown  in  figure  23.  One 
was  used  by  the  controller  and  the  other  served  as  a  backup 
(connected  to  the  Daytronic  controller  for  periodic 
aonitoring) . 

The  Tektronix  4051  aicrocoaputer  was  used  as  the 
data-acqu is  1 1 ion  systea  for  the  tests.  The  Tektronix  was 
prograaed  to  take  electric  potential  readings  at 
predeterained  tiae  intervals.  The  tiae  input  signal  sent  to 
the  Tektronix  was  provided  by  the  Daytronic  controller. 


Figure  22  Creep  Frame  Photograph. 
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Similarly,  the  Hewlett  Packard  3478A  IEEE-488  prograaable 
voltmeter  provided  the  voltage  reading  when  requested  by  the 
Tektronix' s  program.  The  two  received  signals  were  then 
recorded  on  magnetic  tape.  This  information  was  later 
transferred  to  the  host  POP  11/24  computer  for  use  in  the 
data  reduction  procedures  to  obtain  crack  length  from 
voltage.  In  addition,  optical  crack  length  measurements 
were  obtained  through  the  viewing  ports  on  each  side  of  the 
resistance  furnace  using  two  Gaertner  traveling 
microscopes.  The  crack  length  measurement  was  displayed  on 
a  digital  readout  to  an  accuracy  of  0.0254  mm.  However,  due 
to  optical  problems  in  determining  exactly  where  the  crack 
ended,  measurements  were  reproducible  to  an  accuracy  on  the 
order  of  0.05  mm.  The  optical  crack  length  measurements 
were  manually  recorded  periodically  along  with  time  and  the 
corresponding  voltage  readings  from  the  electric  potential 
system  in  order  to  provide  data  for  later  correlation  of 
crack  length  to  voltage. 

Crack  Measurement  Procedure 

The  proof  experiments  utilized  a  D.C.  Electric 
Potential  (EP)  measurement  system  augmented  with  optical 
readings  to  monitor  crack  length.  The  EP  system  is  based  on 
the  fact  that  when  a  current  is  passed  through  a  conducting 
body  an  electric  field  is  generated.  The  field  shape  and 
intensity  depend  upon  such  factors  as  applied  current. 
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geometry  shape  and  material  properties.  The  EP  system 
relies  on  relating  changes  in  the  EP  field  at  the  output 
leads  to  changes  in  geometry  due  to  crack  growth. 


One  of  the  instrumented  specimens  is  shown  in  figure 
23.  The  input  and  output  leads  were  welded  to  the  specimen 
at  the  locations  shown  in  figure  24.  The  input  leads,  made 
of  Inconel  718,  were  connected  to  a  constant  10.0  amp 
current  source.  The  output  voltage  was  measured  from 
nichrome  wire  leads  spot  welded  on  the  specimen  front 
surface.  Nichrome  wire  was  used  for  the  output  leads  due  to 
its  superior  oxidation  resistance  at  elevated  temperature. 
However,  joining  two  dissimilar  materials  produces  a 
thermocouple  effect  as  the  joint  temperature  is  changed. 

This  thermal  voltage  adds  algebraically  to  the  voltage 
generated  due  to  the  resistance  in  the  specimen.  It  should 
be  noted  that  the  thermal  voltage  is  present  even  when  the 
applied  input  current  is  removed  from  the  specimen.  It  was 
therefore  possible  to  periodically  measure  the  thermal 
voltage  by  shutting  the  current  supply  off.  Plots  of  the 
output  voltage  versus  thermal  voltage  (Vth)  were  generated 
for  each  test  specimen  and  are  shown  in  figures  25  through 
27.  The  data  were  fitted  with  a  linear  equation  using  a 
least  squares  regression  to  yield  thermal  voltage  as  a 
function  of  output  voltage  reading.  This  thermal  voltage 


was  calculated  at  each  data  point  and  subtracted  from  the 
output  voltage  reading. 


Speclaen  84-*504~EE5  was  tested  to  provide  a  calibration 
curve  relating  crack  length  to  output  voltage  under 
conditions  of  sustained-load  crack  growth  with  no  applied 
overloads.  The  initial  voltage  reading*  minus  thermal 
voltage*  was  called  the  initial  reference  voltage  Vo.  The 
initial  crack  size  corresponding  to  Vo  was  measured  using  a 
method  similar  to  the  one  recommended  in  ASTM  E399  (161.  An 
optical  microscope*  equipped  with  a  measurement  table*  was 
used  to  make  a  five  point  throagh-the-thlckness  measurement 
of  the  Initial  crack  tip  profile  on  the  ruptured  specimen. 
The  weighted  average*  defined  in  figure  28*  was  used  to 
calculate  an  average  crack  size  and  a  crack  tunneling 
correction  factor  to  account  for  the  increased  crack  depth 
at  the  center  of  the  crack  front.  The  tunneling  correction 
factor  was  used  to  bias  the  surface  optical  measurements* 
made  at  locations  a^  and  a^  in  figure  28*  to  obtain  a 
average  crack  depth.  Calibration  specimen  84-504-EE5  had  an 
average  tunneling  value  of  0.445  am.  This  value  was  added 
to  each  opt  leal . reading  to  get  a  corrected  value. 

A  functional  relationship  between  a^_.  and  voltage 

opt 

readings  (minus  thermal  voltage*  and  normalized  to  the 
referenced  voltage  Vo)  was  developed  using  a  least  squares 
polynomial  fitting  program.  H.  H.  Johnson  (171  developed 
complex  functions  for  calibrating  crack  length  measurements 
to  voltage.  However*  Johnson's  results  showed  that  the 
crack-starter  geometry  strongly  influences  the 


V 


calibration  between  crack  length  and  Measured  potential. 
Therefore  for  simplicity,  an  eighth  order  polynomial  was 
used  which  was  found  to  adequately  fit  the  experimental  data 
and  is  shown  in  figure  29.  The  polynomial  equation  for  crack 
length  as  a  function  of  normalized  voltage  Is, 

a  =  -332.1731  ♦  1503.861  <V/Vo)^  -  2932.647  <V/Vo>^ 

+  3220.104  <V/Vo)^  -  2176.264  <V/Vo>^ 

5  6  ” 

+  927.0209  <V/Vo)®  -  243.1114  (V/Vo>® 

+  35.90090  <V/Vo)^  -  2.286718  <V/Vo)*  (Inch) 


The  Inverse  of  this  function,  normalized  voltage  as  a 
function  of  wais  also  generated  and  Is  shown  In  figure 

30.  The  corresponding  polynomial  equation  Is, 

V/Vo  ■  -29.66971  +  364.3393  (a)*  -  1859.148  <a>^ 

♦  5319.197  (a)®  -  9330.943  (a>^ 

5  6  ^29) 

+  10274,08  <a)^  -  6926.514  (a)® 

+  2612.860  <a>^  -  422.2600  (a>^ 

where  (a)  is  measured  in  inches. 


•  •  •  »  •  • 
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This  equation  was  used  to  calculate  the  Initial  voltage 
(Vo)  for  specimens  with  different  initial  crack  sizes. 
Knowing  the  Initial  crack  size  and  corresponding  voltage 
reading  (V),  an  effective  Vo  was  calculated  using  equation 
29.  Substituting  the  ratio  of  V/Vo^^^  into  equation  28 
yields  the  correct  Initial  flaw  size  for  the  specimen. 
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In  the  tests  Involving  overloads.  It  was  found  that 
each  overload  caused  the  versus  voltage  to  shift  fro* 

the  calibration  curve.  This  effect  is  shown  in  figure  3i. 
After  each  overload,  there  is  a  period  where  voltage 
increases  without  an  appreciable  Increase  in  the  optical 
crack  size.  One  explanation  is  that  the  crack  changes  shape 
with  the  edge  breaking  through  during  the  overload  cycle. 
This  would  account  for  the  juap  in  optical  aeasureaent  after 
each  overload.  As  the  crack  starts  to  grow  again,  evidenced 
by  increasing  voltage  readings,  the  surface  growth  nay  lag 
since  the  crack  nay  tunnel  to  restore  its  preferred  flaw 
shape.  Although  the  exact  nechanisn  that  cause  the 
deviations  fron  the  calibration  curve  are  not  fully 
understood,  their  effect  nay  be  negated  by  recalibrating  the 
equation  after  each  overload.  The  procedure  is  the  sane  one 
used  to  calibrate  the  equation  to  different  initial  flaw 
sizes.  Using  equation  29  and  the  a^^^  and  voltage  reading 
taken  after  each  overload,  a  new  Vo  was  calculated.  This  Vo 
was  used  to  nornallze  voltage  readings  until  the  next 
overload. 

The  experinental  voltage  readings  recorded  by  the 
Tektronix  4051  conputer  were  reduced  using  the  procedures 
described  in  this  section.  The  resulting  crack  growth 
histories  were  conpared  with  the  analytical  predictions  and 
are  discussed  in  the  next  section. 


Eiperiaents  were  performed  in  the  course  of  this  study 
in  order  to  obtain  additional  crack  growth  data  for 
evaluating  the  prediction  capability  of  the  retardation 
models.  An  electric  potential  crack*measureaent  system  was 
used  for  taking  crack  length  measurements  during  sustained 
loading  with  periodic  overloads.  The  applicability  of  this 
system  to  crack  growth  following  overload  was  investigated. 

Three  specimens  were  used  in  the  experimental  work. 

One  was  used  as  a  calibration  specimen  for  the  electric 
potential  system.  The  other  two  specimens  were  used  for  the 
proof  tests.  In  addition*  a  third  proof  test*  conducted  by 
Harms*  was  also  used  to  evaluate  the  retardation  models. 

The  crack  growth  prediction  for  the  calibration 
specimen  84*-504  is  shown  in  figure  32.  This  specimen  was 
subjected  to  sustained  loading  at  650  C  with  no  overloads. 

As  expected*  both  the  Overload  and  Wheeler  models  predicted 
the  same  growth*  since  no  retardation  cycles  were  applied. 
The  analytical  prediction  for  time  to  failure  was  within  20% 
of  the  actual  test  data.  This  is  well  within  the  normal  2X 
scatter  in  crack  growth  data  associated  with  variations  in 
material  properties. 

Specimen  84-503  was  tested  with  20  %  overloads  applied 
each  hour.  The  test  data  and  retardation  model  predictions 
are  shown  in  figure  33.  The  test  was  divided  into  two 
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segments*  due  to  an  overnight  shutdown  In  the  testing.  In 
the  first  segment*  the  crack  growth  appears  to  be  dominated 
by  the  crack  jump  due  to  the  overload  cycle.  The  Overload 
model  overestimated  this  jump*  while  the  calculated  growth 
Increment  from  the  0.01  Hz  overload  cycle  before  the  Wheeler 
model  is  applied  underestimated  it.  During  the  overnight 
shutdown*  the  load  was  removed  from  the  specimen  while  the 
temperature  was  held  constant  at  650  C.  Since  the  specimen 
experienced  thermal  soaking  during  this  shutdown*  reloading 
of  the  specimen  was  followed  by  at  least  an  hour  of 
sustained  load  growth  to  re-establish  the  previous  growth 
rate.  The  overnight  shutdown  provided  a  heat  tinted  region 
visible  on  the  fracture  surface  after  the  specimen 
fractured.  Optical  readings  were  taken  of  the  crack  tip 
profile  to  adjust  the  tunneling  correction  factor  which  was 
then  added  to  the  optical  surface  measurements.  The  Wheeler 
model  overestimated  the  total  time  to  failure  for  the  second 
segment  by  1  %.  The  Overload  model  underestimated  the  growth 
In  the  same  segment  by  approximately  16  %  . 

Specimen  84-502  was  tested  with  20  \  overload  cycles 

applied  at  low  30  <MPa  m^^^)*  medium  (40  MPa  m^^^>  and  high 
1/2 

50  (MPa  m  >  stress  Intensity  levels.  The  effects  of  the 
jumps  on  the  crack  length*  due  to  overload  cycles*  was 
minimized  by  applying  just  three  overload  cycles  during  the 
test.  The  test  data  and  retardation  predictions  are  shown 
In  figure  34.  Both  models  predicted  the  total  time  of 
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Figure  34  Crack  Length  versus  Tlee  Specleen  84-502. 
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Figure  35  Crack  Length  versus  Tlse  Specisen  84-507 


growth  within  10  %  of  the  test  data.  It  was  also  noted  that 
the  retardation  predictions  were  similar  In  shape  to  the 
test  data.  In  general »  the  models  gave  better  predictions 
of  total  tlme-to-fallure  when  the  effect  of  the  crack  jumps 
was  minimized. 

In  addition  to  the  tests  conducted  as  part  of  this 
study,  the  proof  test  conducted  by  Harms  was  analyzed.  The 
test  data  and  retardation  predictions  are  shown  In  figure 
35.  The  test  spectrum  consisted  of  20  %  overload  applied 
each  hour.  The  prediction  of  the  Overload  model  was  within 
4  %  of  the  total  tlme-to-fallure,  but  was  dependent  on  the 
jump  In  crack  length  to  account  for  crack  growth  at  lower 
stress  intensity  levels.  The  Wheeler  model  again 
underestimated  the  growth  due  to  the  overload  cycles  and 
therefore  predicted  a  growth  time  longer  than  the  actual 
data. 

Accurate  prediction  of  the  retardation  caused  by 
overload  cycles  requires  empirically  defining  the  a 
parameter  for  the  Overload  model  and  the  shaping  exponent  a 
for  the  Wheeler  model.  These  variables  were  defined  for  a 
specific  overload  ratio,  temperature  and  test  geometry. 

Both  01  and  m  were  found  to  depend  upon  the  stress  Intensity 
level  at  overload  application.  Using  the  predefined  values 
of  01  and  m  from  section  IV,  both  models  predicted  the 
retardation  effect  of  overload  cycles,  when  the  cycles  were 
applied  far  enough  apart  so  that  no  interations  occured 


between  cycles.  When  the  overload  cycles  were  applied  at 
closer  intervals,  the  crack  growth  was  doninated  by  the 
Juaps  in  the  crack  length  caused  by  the  overload  cycles.  As 
seen  in  the  predictions  in  figures  33  and  35,  the  Overload 
model  did  not  accurately  in  predict  the  Jump  in  crack  length 
when  overload  cycles  were  applied.  Defining  a  function 
relating  the  jump  in  crack  length  to  the  current  stress 
intensity  level  was  difficult  due  to  the  large  scatter  in 
crack  length  juap  aeasureaents .  An  attempt  was  made  to 
reduce  the  scatter  in  crack  length  jump  measurements  by 
using  an  electric  potential  system  to  measure  the  crack 
length.  This  attempt  was  unsuccessful  mainly  due  to  the 
electric  potential  system  requiring  recalibration  via  an 
optical  crack  measurement  after  each  overload.  Optical 
crack  measurements  on  a  specimen  under  sustained  loading  are 
difficult  because  the  crack  tip  is  not  sharply  defined. 
Instead,  the  crack  tip  looks  like  a  deformation  zone,  with 
the  exact  tip  location  unknown.  The  accuracy  of  both  models 
for  predicting  crack  growth  when  the  overload  cycles 
interact  could  be  improved  by  further  refinement  of  the  jump 
function. 

Overall  both  retardation  programs  predicted  the  total 
time  to  failure  within  the  normal  2X  scatter  associated  with 
variations  in  material  properties.  Using  the  CRACKS  program 
it  is  possible  to  analyze  spectra  of  engine  fatigue  cycles 
with  the  sustained-load  growth  between  cycles  included. 


VII. 


During  the  course  of  this  study  several  observations 
were  make  on  the  applicability  of  using  existing  retardation 
Models*  developed  for  airframe  application*  to  predict 
sustained  load  crack  growth  retardation.  It  was  found  that: 

1. )  The  crack  growth  rate  for  sustained  loading  (da/dt) 
and  overload  fatigue  cycles  <da/dn>  can  be  represented  by 
one  crack  growth  rate  equation.  This  is  accoaplished  by 
modeling  the  sustained  load  time  as  equivalent  fatigue 
cycles.  The  equivalent  cycle's  period  and  R  ratio  are 
adjusted  to  obtain  the  overload  fatigue  cycle's  crack  growth 
rate.  The  equivalent  sustained  load  fatigue  cycles  and 
overload  cycles  were  analyzed  using  the  CRACKS  crack  growth 
prediction  program  developed  for  airframe  cyclic  loading. 

The  CRACKS  program  is  capable  of  predicting  the  total  time 
to  failure  within  20  %  of  experimental  data. 

2. )  The  Overload  and  Wheeler  retardation  models  depend 
on  empirical  parameters  that  are  related  to  the  stress 
intensity  level  at  overload  application. 

3. )  The  jump  function  dominates  the  crack  growth  in  the 
Overload  model  when  the  overload  cycles  are  spaced  close 
enough  to  interact. 

4. )  The  Willenborg  retardation  model  modifies  the 
fatigue  cycle's  R  ratio  when  accounting  for  retardation. 
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Thus>  the  retardation  model  was  dependent  on  the  R  ratio 
chosen  for  the  equivalent  sustained  load  fatigue  cycles  and 
was  deemed  unacceptable  for  use. 

5.)  The  electric  potential  crack  measurement  system  is 
affected  by  overload  cycles  and  requires  recalibration  after 
an  overload  cycle  is  applied. 

This  study  verified  that  the  Overload  retardation  model 
predicts  sustained-load  crack  growth  with  periodic  overload 
within  normal  test  data  scatter.  In  addition*  procedures 
were  developed  to  convert  sustained  loading  into  equivalent 
fatigue  cycles  and  analyze  crack  growth  using  the  CRACKS 
program.  The  CRACKS  program*  with  the  Wheeler  retardation 
model  selected  to  account  for  retardation  effects*  is 
capable  of  predicting  the  total  time  to  failure  within  20  \ 
for  tests  consisting  of  sustained-load  with  periodic 
overloads. 

The  modified  CRACKS  program  offers  the  unique 
capability  to  analyze  sustained-loading  and  fatigue  cycle 
loading  together..  This  capability  can  be  readily  applied 
to  complex  engine  spectra*  consisting  of  fatigue  cycles  with 
hold  times*  to  predict  the  crack  growth  in  engine 


components . 


In  this  study  it  was  found  that  the  electric  potential 
systea  for  Measuring  crack  length  was  affected  by  periodic 
overloads.  Additional  investigation  is  required  to 
understand  how  electric  potential  crack  neasurenents  made  on 
fatigue  loaded  speciaens  with  periodic  overloads  coapare 
with  sustained-loaded  speciaens  with  periodic  overloads. 

This  study  used  Inconel  718  material  exclusively  for 
experiaental  testing.  Additional  testing  should  be 
perforaed  to  determine  if  the  models  apply  to  other 
materials.  Another  nickel-base  superalloy  such  as  Rene  95 
is  recoaaended. 

Finally,  the  retardation  parameters  for  the  Overload 
and  Wheeler  models  should  be  developed  for  a  wider  range  of 
overload  ratios  and  used  to  analyze  a  more  complex  spectrum. 


V 


Blbl lograohy 


1.  Hlll>  R.  J.«  Reiaann,  W.  H.,  and  Ogg,  J.  S.,  ‘A 
Ret ireaent-For-Cause  Study  of  an  Engine  Turbine  Disk:, 
AFWAL-TR-81-2094.  Wr Ight-Patterson  APB,  OH  1981. 

• 

2.  Harris.  J.  A..  Sins.  0.  L.,  and  Ann  is.  C.  G..  "Concept 
Definition:  Retlreaent  for  Cause  of  FIDO  Rotor  Coaponents”. 
AFWAL-TR-80-41 18.  Ur ight-Patterson  AFB.  OH  1980. 

3.  Military  Standard:  Engine  Structural  Integrity  Prograa, 

^  MIL-STD-1783  (USAF)  30  Noveaber  1984.) 

4.  Engle,  R.M.,  "CRACKS,  A  FORTRAN  IV  Digital  Coaputer 
Prograa  for  Crack  Propagation  Analysis,"  AFFDL-TR-70-107, 

Air  Force  Flight  Dynaalcs  Laboratory,  Ur Ight-Patterson  Air 
Force  Base,  Oh,  Oct.  1970 

5.  Uheeler,  O.E.,  "Spectrua  Loading  and  Crack  Growth", 
Journal  of  Basle  Engineering.  March  1972,  pp  181-186. 

6.  Ulllenborg,  J.,  Engle,  R.M.,  and  Uood,  H.A.,  "A  Crack 
Growth  Retardation  Model  Using  an  Effective  Stress  Concept", 

^  AFFDL-TM-71-1-FBR,  Ur Ight-Patterson  AFB,  Oh,  Jan  1971 

9 

7.  Haras  K.  E.,  "Overload  Effects  on  Sustained  Load  Crack 
Growth  at  Elevated  Teaperatures* ,  AFIT  Thesis, 
AFIT/GAE/AA/84D-8,  Ur Ight-Patterson  AFB,  Oh,  1984. 

8.  Ueerasoor lya,  T.,  Nicholas,  T.,  "Overload  Effects  in 

9  Sustained  Load  Crack  Growth  in  Inconel  718", 

AFUAL-TR-85-4121 .  Ur Ight-Patterson  AFB,  OH  1985. 

9.  Annual  Book  of  ASTM  Standards,  Part  10; 

Metals-Mechanlcal ,  Fracture,  and  Corrosion  Testing;  Erosion 
and  Uear;  Effect  of  Teaperature,  E  647-83,  page  719,  1983. 

10.  Utah,  0.  A.  General  Electric  Coapany,  "Crack  Growth 
Modeling  In  an  Advanced  Powder  Metallurgy  Alloy", 
AFUAL-TR-80-4098.  Ur ight-Patterson  AFB,  OH  1980. 

11.  Kreyszig,  E.,  EMiMtrtna  Fl^th 

Edition,  New  York,  John  Uiley  and  Sons,  1983,  pp  788-789. 

12.  Larsen  J.  M.  and  Nicholas,  T.  "Load  Sequence  Crack 
Growth  Transients  in  a  Superalloy  at  Elevated  Teaperature", 
Fracture  Mechanics:  Fourteenth  Svaposlua  Voluae  II;  Testing 
and  AddI icat Ions .  ASTM  Society  for  Testing  and  Materials, 
1983,  pp  II-536-II-552. 


-79- 


13.  Harltos,  G.K.*  Nicholas*  T.  and  Painter*  G.O.* 
"Evaluation  of  Crack  Growth  Models  for  Elevated  Temperature 
Fatigue*"  Eighteenth  National  Symposium  on  Fracture 
Mechanics*  ASTM  STP  XXX*  American  Society  for  Testing  and 
Materials,  1985*  pp  x*-yy 

14.  Engle*  R.  M.  "Damage  Accumulation  Techniques  in  Damage 
Tolerance  Analysis"*  ASXtf,^TP_842*  American  Society  for 
Testing  and  Materials*  1984*  pp.  25*35. 

15.  Collins*  J.  A.*  Failure  of  Materials  In  Mechanical 
Design.  New  York*  John  Wiley  and  Sons*  1981*  Ch.  8. 

16.  Annual  Book  of  ASTM  Standards*  Part  10; 
Metals-Mechanical *  Fracture*  and  Corrosion  Testing*  Erosion 
and  Wear;  Effect  of  Temperature*  1981. 


17.  Johnson*  H.  H.*  "Calibrating  the  Electric  Potential 
Method  for  Studying  Slow  Crack  Growth"*  Materials  Research 
and  Standards,  Vol.  5,  No.  9,  1965 


Appendix  A 


Heat  Treatment  History  of  Test  Soecleens 


Anneal  at  968  C  for  1  hour  -  air  cool 


Age  harden  at  718  C  for  8  hours  **  furnace  cool  to  621  C 


Age  harden  at  621  C  for  a  total  of  an  additional  10  hours 
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Append!  x  g., 
OVERLOAD  PROGRAM 


C  PROGRAM  OVERLD 

DOUBLE  PRECISION  V 

REAL«4  LOWER, KUP, KLOU, N, LOUM, KINIT, LOUERR, KLOUU, 

1KUPP,L0WR,L0AD,K 
INTEQER*4  PIECES, D,YNCR,NAME<7) 

DIMENSION  TITLE(60),SBTITL(60),TIME(20),0L(20) 

COMMON/LYI A1 /PI , Q, TOL, U, B, T, R, 

IQUE, PEA, BEE, DEE, DKSTAR, TOUGH, DKCRIT, 

2TR1G, LOMU, DAGER, NSKIP 
C  DEFAULTS  FOR  DIAGNOSTIC  OPTION 

DATA  V/0/ , PI /3. 1 4 1 59/ , TOL/ . 0001 /, R/ . 0/ , ZERO/0. 0/ , T/ .  25/ ,  TR I G/ 1 .  / 

C  *****  OPEN  FILES  2: INPUT  IsOUTPUT  3;PL0T  ************************* 
WRITE  (5,40) 

40  FORMAT  (»  WHAT  IS  YOUR  INPUT  FILE?  »,♦) 

READ  (5,90)  NAME 

OPEN  <UN1T=2,NAME*NAME,TYPE=*0LD') 

WRITE  (5,50) 

50  FORMAT  (/,'  WHAT  SHOULD  YOUR  OUTPUT  FILE  BE  CALLED?  ',») 

READ  (5,90)  NAME 

OPEN  (UNIT=1,NAME*NAME,TYPE=»NEW') 

WRITE  (5,60) 

60  FORMAT  (/,’  WHAT  SHOULD  YOUR  PLOT  FILE  BE  CALLED?  ',*) 

READ  (5,90)  NAME 

OPEN  (UNIT®3, NAME*NAME, TYPE=»  NEW’ ) 

90  FORMAT  (11A2) 

*****  all  files  ATTACHED  *************************************** 

*****  READ  INPUT  FILE  ****************************************** 

READ  (2,1120)  TITLE 
READ  (2,1120,END=1590)  SBTITL 
READ  (2,*)  TOUGH, LOWW,R, YIELD 
READ  (2,*)  BEE, DKSTAR, PEA, QUE, DEE 
READ  (2,*)  B,LOAD,W 
WRITE  (5,105) 

WRITE  (5,110) 

105  FORMAT  (/,’  DO  YOU  WANT  TO  BREAK  OUT  ON  A  CRACK  LENGTH’) 

110  FORMAT  (’  OR  AT  THE  KIC  STRESS  INTENSITY  C1=A,  0=KIC]— >',*) 

READ  (5,*)  AFIN 

IF  (AFIN.EQ.l)  GO  TO  150 

AFIN*10000. 

GO  TO  160 
150  WRITE  (5,155) 

155  FORMAT  (/,’  INPUT  YOUR  FINAL  CRACK  LENGTH,  INCHES— >’,♦) 

READ  (5,*)  AFIN 

*****  END  OF  INPUT  ******************************************** 

*****  PRINT  OUTPUT  FORMATS  ON  FILES  *************************** 


-"•^rvT^lL^K^  *4*.  »■  ’ 


ISO 

WRITE 

(1, 

1130)  TITLE 

WRITE 

(3, 

1130)  TITLE 

WRITE 

(5, 

1130)  TITLE 

WRITE 

(1, 

1130)  SBTITL 

WRITE 

(3, 

1130)  SBTITL 

WRITE 

(5, 

1130)  SBTITL 

WRITE 

a, 

740)  TOUGH, R 

WRITE 

<1, 

850)  YIELD 

WRITE 

(1. 

830)  B,W 

WRITE 

(1, 

930)  BEE,DKSTAR,PEA 

WRITE 

(1, 

940)  DKSTAR,QUE 

WRITE 

<1, 

950)  TOUGH, DEE 

WRITE 

(1, 

1270) 

WRITE 

(1, 

1280) 

740 

FORMAT 

(' 

FRACTURE  TOUGHNESS: ' ,F6. 2, /, '  STRESS  RATIO; ' ,F5.3) 

830 

FORMAT 

(' 

THICKNESS:  ',F7.5,/,'  WIDTH: ' ,F7.5) 

850 

FORMAT 

(' 

YIELD  STRENGTH; ',F6. 2) 

930 

FORMAT 

(' 

da/dN=EXPC ' , F6. 2, ' ]*C  CDELTA-K/' , F6. 2, ' ]*»' , F6. 2, ' ]*' ) 

940 

FORMAT 

(' 

[ [ 1 nIDELTA-K/' , F6. 2, ' 3 1**' , F6. 2, ' 3*' ) 

950 

FORMAT 

(' 

CClnC',F6.2,'/DELTA-K33*»',F6.2,'3') 

1120  FORMAT  <60A2) 

1130  FORMAT  (1X,60A2) 

1270  FORMAT  (35X,'  d*lta») 

1280  FORMAT  (IX,'  N(xlOOO) ' , lOX, 'a' ,9X, 'LOAD' ,2X, 
1'  K  ',1X,'  da/dN  ', IX, 'PIECES* ) 

C  END  OF  HEADERS 


C 

C  ttttt  READ  IN  SPECTRUM  LOAD 


NA*1 


95  READ  (2,*,END=100)  TIME(NA) , OLCNA) 
NA=NA+1 


GO  TO  95 
100  NA=NA-1 

C  tttt**tt*t*t**ttttil^**t**t*tt*ttttt*tttttttttttt***t*******t*t**tt 
T0UGH=(1.0-R)*T0UGH 
DKCRIT=TOUGH 
T0UQH=T0U6H*.95 

1210  STAR=1.0 
KUP=0 
V=TIME(1) 

SIGMA»LOAD 

C  *****  START  OF  SPECTRUM  INTEGRATION  ***************************** 
DO  1570  NS=1,NA 
DAGER=OL(NS) 

NSKIP=0 


D0L=0 
U=LOUU 
1290  LOWER=U 

C  ****************************************************************** 

c 

1320  CALL  SIMPCSUM, LOWER, UPPER, SIGMA, KLOH,KUP, PIECES, YIELD, STAR, PYTCH, 
1CALPHA,BETA,D0L) 


4  * 
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% 


c 

C  ttttt  Check  if  initial  K  is  greater  than  Klc 
IF  (KLOW.LT. TOUGH )t30  TO  1360 
WRITE  (5,1330) 

WRITE  (1,1330) 

1330  FORMAT  dOX,’  **«*SPECIMEN  FAILED  ON  LOADINGm**' ) 

GO  TO  1490 
C 

1360  VN=V 

V=V+ (SUM/1000) 

C  *****  CHECK  IF  INTEGRATED  TIME  IS  MORE  THAN  THE  NEXT  OVERLOAD  * 

C 

IF  (NS.EQ.NA)  GO  TO  1370 

IF  (V.LT.TIME(NS+1))  GO  TO  1370 

IF  (NSKIP.EQ.l)  GO  TO  1370 

PYT I T=PYTCH* ( T I ME (NS+ 1 ) -VN) / ( SUM/ 1000. ) 

UPPER=LOWER+PYTIT 

C  ************  EXTRA  PRINT  TO  CHECK  LIhEAR  PYTCH  INCREMENT  ***** 

C  WRITE  (1,1331)  V,VN, NS, PYTIT, LOWER, UPPER, SUM 

C1331  FORMAT  (IX, 'V=  »,F15.6,'  VN  =  »,F15.6,»  NS=  ’,12,/,'  PYTIT=  ’, 
C  1F15. 10,/,’  LOWER=  ’ ,F16. 13, 'UPPER=  ’,F16. 13,’  SUM  =  ’,F9.3) 

C  *************************************************************** 

V=VN 
NSKIP=1 
130  TO  1320 
C 

1370  CALL  FINDF2  (KL0W,F2) 

DADNL=1/F2 

CALL  FINDF2  (KUP,F2) 

DADNU=1/F2 

IF  (U.EQ.LOWW)  WRITE  (1, 1421) VN, LOWER, SIGMA, KLOW,DADNL, PIECES 

C 

C  *****  CONVERSION  FROM  ENGLISH  TO  METRIC  TO  WRITE  PLOT  FILE  ******* 

C0ETA=BETA 
CBETA=CBETA/25.4 
C  1/INCHES— >1/MILLIMETERS 

L0WERR=L0WER*25. 4 
C  INCHES— >MILLIMETERS 

DADNLL=DADNL*0. 0254 

C  INCHES  PER  CYCLE— >METERS  PER  CYCLE 

C 

KLOWWsKLOU* 1.0989 

C  KSI  ROOT  INCHES— >MPA  ROOT  METERS 

C 

UPPERR=UPPER*25.4 
DADI«JU=DADNU*0. 0254 
VV*V*1000 
KUPP=KUP*1 . 0989 
VVN=VN*1000. 

C  IF  (U.EQ.LOWW)  WRITE(3,1375)  CALPHA,CBETA 
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Vi.-.T.V. 


% 


•  •  •  .  . 


v.v 


^  / 


o  o  o 


IF  (U.EQ.LOUU)  URITE  (3,«>  VVN,LOUERR,DADNLL,KLOUU 
IF  (U.EQ.LOWW)  WRITE  (5, 1420) VN, LOWER, SIGMA, KLOW,DADNL, PIECES 
WRITE  (1,1420)  ,V, UPPER, SIGMA, KLff>,DADNU, PIECES 
WRITE  (3,*)  ,VV,UPPERR,DADNUU,KljPP 
WRITE  (5,1420)  ,V, UPPER, SIGMA, KUP,DADNU, PIECES 
C1375  FORMAT(»  ALPHA  =  ',Fa.5,/,'  BETA  =  ',F8.5) 

1420  FORMAT  (1X,F15.3,F9.6,F9.2,F7.2,E10.3, 17) 

1421  FORMAT  ( IX, F15. 4, F9. 6, F9. 2,F7.2,E10. 3, 17) 

1430  IF  ( (KUP.6T. TOUGH). OR. (UPPER. GE.AFIN))  GO  TO  1490 
IF  (NSKIP.EQ.DQO  TO  1560 
U=U+PYTCH 

CALL  CPYTCH  (KUP,PYTCH) 

GO  TO  1290 
1490  WRITE  (1,1500) 

1500  FORMAT  (6X,»  C-T  SPECIMEN,  SANS  PLASTIC- ZONE  CORRECTION') 

GO  TO  1590 

*****  ENTER  JUMP  FUNCTION  ************************* 

1560  XJUMP=.015 
LOWW=U-t-XJUMP 

C 

WRITE  (1,1561)  XJUMP,KUP 

1561  FORMATdX, 'XJUMP  =  ’,F15.10,'  KUP=  ',F6.2) 

1570  CONTINUE 

C  ********************  END  DO  LOOP  *************** 

1590  CALL  EXIT 
END 

C  SUBROUTINES  IN  DESCENDING  ORDER  OF  USE, FIRST, SIMPSON'S  RULE  APPROX 
SUBROUTINE  SIMP(SUM, LOWER, UPPER,SIGMA,KLOW,KUP, PIECES, YIELD, STAR, 
IPYTCH, CALPHA, BETA, DOL) 

COMMON/LY lAl /P I , Q, TOL, W, B,  T, R, 

IQUE, PEA, BEE, DEE, DKSTAR, TOUGH, DKCRIT, 

2TRIG, LOWW, DAGER, NSKIP 
REAL*4  K,N, LOWER, KLOW,t<UP, LOWW 
INTEGER*4  PIECES 
PIECES=2 
X=LOWER/T 

CALL  CT (SI GMA, X, K, F2, YIELD, STAR, PYTCH, CALPHA, BETA, DOL) 

IF  (X.EQ.LOWW/T) 

ICALL  CALCSIQMA, X, K, F2, YIELD, STAR, PYTCH, CALPHA, BETA, DOL ) 

KLOW=K 

IF  (X.EQ.LOWW/T)  CALL  CPYTCH(K, PYTCH) 

IF  (NSKIP. NE.l)  UPPER=LOWER+PYTCH 
ESUM=F2 

DELTA= ( UPPER-LOWER ) /P I ECES 
EVSUM=0 

X=(LOWER+DELTA)/T 

CALL  CT (SIGMA, X, K, F2, YIELD, STAR, PYTCH, CALPHA, BETA, DOL) 

0DSUM=F2 

X=UPPER/T 

CALL  CT (SIGMA, X, K, F2, YIELD, STAR, PYTCH, CALPHA, BETA, DOL) 


KUP=K 

ESUM=ESUM+F2 

SUM=(ESUM+4*0DSUM)*DELTA/3 
1600  PIECES=PIECES*2 
SUM1=SUM 

DELTA= ( UPPER-LOWER) /P I ECES 

EVSUM=EVSUM+QDSUM 

0DSUM=O 

L=IFrx (FLOAT (PIECES) /2) 

DO  1610  1=1,10000 
Z  =LOWER+ DELTA* ( 2* I - 1 ) 

X=Z/T 

CALL  CT (SIGMA, X, K, F2, YIELD, STAR, pyT(>l, CALPHA, BETA, DOL) 
0DSUM=0DSUM+F2 
1610  IF  (I.EQ.L)  GO  TO  1620 
1 620  SUM= (ESUM+4*0DSUM+2*EVSUM) *DELTA/3 

IF  (ABS( SUM-SUM 1).GT.ABS(T0L*SUM))  GO  TO  1600 

RETURN 

END 

C  SUBROUTINE  FOR  CALCULATING  CALPHA  FOR  THE  RETARDATION 
C  MODEL  (20y.  AND  SOX  OVERLOAD  CASES) 

SUBROUTINE  CAL (SIGMA, X , K, F2, YIELD, STAR, PYTCH, CALPHA, BETA, DOL) 
COMMON/LYIAl/PI , Q, TOL, W,B, T,R, 

1 QUE , PE A , BEE , DEE , DKSTAR , TOUGH , DKCR I T , 

2TR IG, LQWW, DAGER, NSKIP 
REAL*4  K,N,LOUU 
IF  (DAGER. £0.0.)  CALPHA=0. 

IF  (DAGER. EQ.O.)  GO  TO  1679 

SK=K/DKSTAR 

CAS=1-1/SK 

IF  (DAGER. EQ. 50)  GO  TO  1675 

BET A= ( SORT ( 2 . 0 ) *P I **2*Y I ELD**2) / ( .  44*K**2) 

CACAS= (-. 73079 1E-01)*SK**3  +  (. 303086) *SK**2  - 
1  (.422108)*SK  +  0.117517E01 
GO  TO  1677 

1 675  BETA= ( SORT ( 2 . 0 ) *P I **2* Y I ELD**2 ) / ( 1 , 25*K**2 ) 

CACAS=(-. 121127E-01)*SK«*2  +  (. 23923 1 E-0 1 )*SK  +  0.987133 
1677  CALPHA=CACAS*CAS 
1679  DaL=l 

CALL  CT (SIGMA, X, K, F2, YIELD, STAR, PYTCH, CALPHA, BETA, DOL) 

RETURN 

END 

C  SUBROUTINE  FOR  COMPACT  TENSION  SPECIMENS 
C  IN  THIS  ROUTINE,  "SIGMA"  IS  A  LOAD,  NOT  A  STRESS!!! 

SUBROUTINE  CT(SIGMA,X,K,F2, YIELD, STAR, PYTCH, CALPHA, BETA, DOL) 
COMMON/LYIAl/PI,a,TOL,W,B,T,R, 

1 QUE , PE A , BEE , DEE , DKSTAR , TOUGH , DKCR I T , 

2TR I G , LOUU , DAGER , NSK I P 
REAL*4  K, Ml, M2, M3, N, LOUU 
ALPHA=X*T/W 
M1=SIGMA/(B*(SQRT(W))) 

M2= ( 2+ALPHA ) / ( ( 1 -ALPHA ) ** 1 . 5 ) 
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M3=0. 806+ (4. 64CALPHA) - ( 13.32* (ALPHA**2) ) + 

1 (14. 72*(ALPHA**3) ) -(5. 6*  <ALPHA»*4) ) 

K=Ml*M2*M3*(l.0-R) 

IF  ((X.EQ.LOWU/T)  .AMO.  (DOL.EQ.O))  GO  TO  1760 
DELA=X»T-LOUW 

STAR*CALPHA«EXP<-BETA*DELA) 

ErFK=K*(l-STAR) 

K=EFFK 

1760  CALL  FINDF2(K,F2) 

RETURN 

END 

SUBROUTINE  F1NDF2(K,F2) 

C  SUBROUTINE  SELECTING  F2  BASED  MODIFIED  SIGMODAL  EON  ONLY 
REAL*4  K,N,LOVM 
INTEGER*4  CINH 

COMMON/LYI A1 /PI , Q, TOL, W, B, T, R, 

IQUE, PEA, BEE, DEE, DKSTAR, TOUGH, DKCRIT, 

2TRI6, LOUU, OAGER, NSKIP 

CINH=2 

C  **  USE  MODIFIED  SIGMODAL  EON  OM.Y  CINH=2  ******* 

IF  <CINH.NE.2)  GO  TO  1840 
IF  (K.QT. DKSTAR)  GO  TO  1830 
URITE  (5,1800)  K 
WRITE  (5,1810)  DKSTAR 
WRITE  (5,1820) 

1800  FORMAT  ('  YOUR  INITIAL  DELTA-K  IS  »,F5.2,»,AND  THIS  IS  SMALL-') 
1810  FORMAT  ('  ER  THAN  YOUR  DELTA-K  n«ESHOLD  OF  »,F5.2,'.  THIS') 
1020  FORMAT  (’  WILL  GIVE  THE  MSE  INDIGESTION.  TRY  AGAIN.') 

CALL  EXIT 
STOP 

1830  F2= 1 . 0/ ( EXP ( BEE ) » ( (K/DKSTAR ) **PEA) * ( ( ALOG 1 0 ( 

IK/DKSTAR) )**QUE)*( (AL0G10(DKCRIT/K) )**DEE) ) 

1040  IF  (CINH.EQ. 1)  F2=l .0/( 10. 0**(S1*(SINH(S2*( (AL0610(K) )+S3) ) ) 
1+S4)) 

IF  (CINH.EQ.O)  F2=(K**(-N))/C 

RETURN 

END 

C  *****  SUBROUTINE  TO  CALCULATE  PYTCH  INCREMENT  FOR  SIMP  ***** 
SUBROUTINE  CPYTCHXK, PYTCH) 

REAL*4  K 

COMMON/LYI A1 /PI , Q, TOL, W, B, T, R, 

IQUE, PEA, BEE, DEE, DKSTAR, TOUGH, DKCRIT, 

2TRI6, LOUU, DAGER, NSKIP 
IF  (K  .LE.  22.77)  PYTCH=. 00001 
IF  (22.77  .LT.  K  .AND.  K  .LE.  30.0) 

1 PYTCH= ( . 002*K- . 045) /53. 5 
IF  (30.0  .LT.  K) 

1 PYTCH= ( . 0004*K- . 0 1 ) / 7 . 0 
RETURN 
END 
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Sample  Input 


TYPE  I502T.DAT 

30  SEPT  85  SPECIMEN  84-502,  EE3  RLH 
SPECIMEN  84-502,  EES  AT  1200  DEG  P 
272.73, .4087,0., 120. 
-8.69,21.00,-1.1,1.8,-1.8 
.394,2.740,1.5736 
0.0, 0.0 
1.426,20.0 
32.888,20.0 
44.707,20.0 


Appendix  C 


Modified  CRACKS  Program 

PROGRAM  CRACKS4( INPUT, OUTPUT, TAPe5=INPUT, TAPE6=0UTPUT , TAPE7, TAPEl >000100 


COMMON/DATA/  EQN, NASA,  JIPR, J2PR, J3PR, J4PR, J5PR, AZERO, AMAX, NZERO  000110 
INTEGER  EQN  000120 

REAL  NZERO  000130 

COMMON/RDATA/  MODEL, RETARD, PLSTRN,OVLD, SIGMAX, SIGMIN, ASUBP, SMALLM  000140 
INTEGER  RETARD, PLSTRN  000150 

COMMON/MDATA/  MATID( 18) ,C,SMALLN,CARRAY( 100) ,SNARAY( 100) ,KSUBC,  000160 
/  KSUBQ,SIGMAY,DELKTH,RMULT,RCUT,OLMAX  000170 

REAL  KSUBC,KSUBQ  000180 

COMMON/LDATA/SMAX<  20, 10),SMIN<  20, 10) , CYCLES!  20, 10) ,NLYRS(10),  000190 

/  NBLKS, IBLKS(  SO  ),ISEGS(  SO  ),NSEGS  000200 

COMMON/STEPS/ISEQ, J1 , J2, J3, J4, J5, ISTOP, NORTRD  000210 

COMMON/CORFAC/  ISURF, RATIO, PHI, THICK, 1BETA( 10) , BETA! 10) ,NPTS,  000220 

/  AOVERB!100),BTABLE!100),NPTS2,AOVRB2!100),  000230 

/  BTABL2!100),ASTART!10),AST0P!10)  000240 

COMMON/DUTPOT/KMAX,KMAXA,DELTAK, 1FLT,DADNPR  000250 

REAL  KMAX,KMAXA  000260 

WRITEf6,3000)  000270 

IRSTRT  s  0  000280 

MODI  =  0  000290 

RTARDl  =  0  000300 

I CASE  »  0  000310 

1  I CASE  =  I CASE  +  1  000320 

ISPEC=0  000330 

5  IFLT=0  000340 

RETARD  =  0  000350 

JSTART=1  000360 

OVLD  =0.0  000370 

IFfRTARDl.EQ.O)  GO  TO  10  000380 

NORTRD  =  0  000390 

MODEL  =  MODI  000400 

10  I  SPEC  =  I  SPEC  1  000410 

CALL  INPUT! ICASE, ISPEC, IRSTRT)  000420 

IF! ISTOP. NE.O)  GO  TO  1  000430 

IF! IRSTRT. GE.4)  GO  TO  1150  000440 

20  A  =  AZERO  000450 

CYC  =  NZERO  000460 

ABLK  =  AZERO  000470 

DO  1000  J1=JSTART, NBLKS  000480 

DO  900  J2=1,NSEGS  000490 

ISEG  =  ISEGS!J2)  000500 

IBLK  =  IBLKS!J2)  000510 

NLYR  =  NLYRS!ISEG)  000520 

DO  800  J3=1,IBLK  000530 

IFLT  =  IFLT  +  1  000540 

DO  700  J4=1,NLYR  000550 


DN  *  CYCLES(J4,ISEG)  000560 

CYCSVE  =  CYC  000570 

500  CALL  6RWCRK(CYC,A,DN)  000580 

IF(A.LT.AMAX)  GO  TO  550  000590 

URITE(6,3600)  A,AMAX,irLT  000600 

60  TO  1100  000610 

550  IF  (IST0P.NE.2)  GO  TO  600  000620 

DN  =  CYCSVE  ♦  DN  -  CYC  000630 

I STOP  =  0  000640 

I SURF  =  0  000650 

IFCDN.GT.O.O)  60  TO  500  000660 

600  IF(ISTOP.NE.O)  GO  TO  1100  000670 

IFdFLT.EQ.  1)  CALL  OUTPUT(CYC, A)  000680 

IF(J4PR.EQ.O)  GO  TO  700  000690 

IF(M0D(1FLT,J4PR).EQ.0)  CALL  QUTPUT(CYC,A)  000700 

700  CONTINUE  000710 

IF(J3PR.EQ.O)  GO  TO  800  000720 

IF(MOD(IFLT,J3PR).NE.O)  GO  TO  GOO  000730 

WRITE(6,3050)  IFLT, ISEG.A  000740 

800  CONTINUE  000750 

IF<J2PR.EQ.O)  60  TO  900  000760 

IF(MOD(J2,J2PR).NE.O)  GO  TO  900  000770 

URITE(6,3100)  J2,J1,A  000780 

900  CONTINUE  000790 

IF(IRSTRT.EQ.3)  CALL  RESTRTICYC, A, IRSTRT)  000800 

IF(JIPR.EQ.O)  GO  TO  1000  000810 

IF<M0D<J1,J1PR).NE.0)  GO  TO  lOOO  000820 

DELA  =  A  -  ABLK  000830 

GROWTH  *  A  -  A2ER0  000840 

ABLK  *  A  000850 

WRITE<6,3200)  Jl, A, DEL A, GROWTH  000860 

1000  CONTINUE  000870 

GROWTH  =  A  -  AZERO  000880 

WR I TE ( 6 , 3500 )  C YC , A , GROWTH  000890 

1100  IF(NORTRD.EQ.O)  GO  TO  1300  000900 

1150  ICHECK  =  IRSTRT  ♦  1  000910 

GO  TO  <1,5, 10, 1,1200, 1250), ICHECK  000920 

1200  CALL  RESTRTCCYC, A, IRSTRT)  000930 

AZERO  »  A  000940 

NZERO  *  CYC  000950 

JSTART  =  Jl  -  1  000960 

GO  TO  20  000970 

1250  CALL  RESTRTCCYC, A, IRSTRT)  000980 

AZERO  »  A  000990 

NZERO  =  CYC  001000 

JSTART  =  Jl  +  1  001010 

NBLKS  *  NBLKS  +  Jl  001020 

60  TO  20  001030 

1300  ISPEC  =  ISPEC  +  1  001040 

WRITE(6,2000)  001050 

WRITE(6,2900)  I CASE, ISPEC  001060 

WRITE(6,3300)  001070 


RTARDl  »  1  001080 

MODI  =  MODEL  001090 

MODEL  =  0  001100 

ISTOP  =*  0  001110 

RETARD  =0  001120 

NORTRD  *  1  001130 

OVLD  =0.0  001140 

lELT  =  0  001  ISO 

GO  TO  20  001160 

2000  FORMAT(lHl)  001170 

2900  F0RMAT<1H1,70<1H»)/26X,5HCASE  ,I2,4HRUN  , I2/lX,70flH»))  001180 

3000  FORMATCIHI,  70<IH«)/1X,  70(1H*)//IX,29(1H»), 12H  CRACKS  IV  29( 1H*)001190 
1//  20X,9HVERS10N  6, 12X,9H09/21/79  /30X, 13HR.M. ENGLE  JR./IBX,  001200 

2  36HAIR  FORCE  FLIGHT  DYNAMICS  LABORATORY/ /26X, lOHAFFDLCFBE)/  001210 

3  26X, 19HATTN(R.M. ENGLE  JR. )/26X,17HU-PAFB, OHIO  45433//  001220 

4  26X,12H513-2S5-6104//1X,  70(1H*)/1X,  70(1H*)  )  001230 

3050  F0RMATC15H  END  OF  FLIGHT  I5,9H  MISSION, 12, 7X, 14HCRACK  LENGTH  =,F  001240 

/10.5)  001250 

3100  FORMATdSH  END  OF  SEGMENT,  16,9H  OF  BLOCK,  IS, 5X,  14HCRACK  LENGTH  =,  001260 
/F10.5)  001270 

3200  F0RMAT(14H0  END  OF  BLOCK, 15, 20X, i4HCRACK  LENGTH  =,F10.5/  001280 

/5X,19HGR0WTH  THIS  BLOCK  =,F10.5,5X, 14HT0TAL  GROWTH  *,F10.5//)  001290 

3300  F0RMAT(///1X,  70( 1H»)/24X,33HRERUN  OF  CASE  WITH  NO  RETARDATION/  001300 
/  IX,  70(1H»))  001310 

3500  FORMATdHO,  70dH*)/  5X,14HT0TAL  CYCLES  =,F12.2,5X,20HF1NAL  CRACK  001320 
/LENGTH  =,F10.5/22X,20HTOTAL  CRACK  GROWTH  =,F10.5/1X,  70dH*)/lHl)  001330 


3600  FORMATdHO,  70dHC)/  1X,20HCURRENT  CRACK  LENGTHF10.5/31H  EXCEEDS  A001340 
/LLOWABLE  CRACK  LENGTH, F10.5/9X17H  IN  FLIGHT  NUMBER, I6/1X,  70dH*) >001350 


END  001360 

SUBROUTINE  OUTPUT (CYC, A)  001370 

C  001380 

COMMON/RDATA/  MODEL, RETARD, PLSTRN, OVLD, SIGMAX , SIGMIN, ASUBP, SMALLM  001390 
INTEGER  RETARD, PLSTRN  001400 

COMMON/LDATA/SMAX(  20, 10),SMIN(  20, 10), CYCLES (  20, 10) ,h«.YRSdO) ,  001410 

/  NBLKS,IBLKS(  50  ),ISE6S(  50  ),NSEGS  001420 

COMMON/MDATA/  MATIDd8),C,SMALLN,CARRAYdOO>,SNARAYdOO),KSUBC,  001430 
/  KSUBa,SIGMAY,DELKTH,RMULT,RCUT,OLMAX  001440 

CQMMON/DATA/EQN, NASA, JIPR, J2PR, J3PR, J4PR, J5PR, AZERO, AMAX , NZERO 
REAL  KSUBC,KSUBQ  001450 

COMMON/STEPS/ ISEG, J1 , J2, J3, J4, J5, ISTOP, NORTRD  001460 

COMMON/OUTPOT/KMAX.KMAXA, DELTAK, IFLT, DADNPR  001470 

REAL  KMAX,KMAXA,KMXEFF  001480 

DATA  NOP/0/ 

IF (RETARD. EQ.O)  GO  TO  100  001490 

C  001500 

C  OBTAIN  OUTPUT  QUANTITIES  FOR  RETARDATION  001510 

C  001520 

GO  TO  (10,20,30,40,50),  MODEL  001530 

10  CALL  WHELER(CYC,A,DDNRET)  001540 

GO  TO  90  001550 

20  CALL  WLNBRG(CYC,A,DDNRET)  001560 

GO  TO  90  001570 


o  o  o  o 


30  DDNRET  =  DADNPR  001580 

GO  TO  90  001590 

40  CONTINUE  001600 

50  CONTINUE  001610 

90  KMXEFF  =  KMAX  001620 

DLKEFF  =  OELTAK  001630 

001640 

DETERMINE  UNRETARDED  QUANTITIES  CORRESPONDING  TO  CURRENT  CRACK  001650 

LENGTH  001660 

001670 

100  CONTINUE  001680 

SIGMAX  =  SMAX(J4, ISEG)  001690 

SIGMIN  =  SMIN<J4,ISEG)  001700 

CALL  RATE(CYC,A,DADN)  001710 

IFCMODEL.NE.O. AND. RETARD. NE.O)  GO  TO  200  001720 

KMXEFF  =  KMAX  001730 

DLKEFF  =  DELTAK  001740 

DDNRET  =  DADN  00 1750 

C  001760 

C  PRINT  HEADER  AT  TOP  OF  EACH  PAGE  001770 

C  001780 

200  CONTINUE  001790 

IF(J4.EQ. 1)  URITE(6,2000)  001800 

IF(J4.EQ.l  .AND.  NQP.EQ.O)  URITE(1,2000> 

IF(J4.E0. 1  .AND.  NOP.EQ.O)  WRITE(1,2001)N2ERO,AZERO 
2001  F0RMAT(16X,F10.1,1X,F9.6) 

N0P=N0P+1 

RETFAC  *  DDNRET /DADN  001810 

IF(DADN  .EQ.  0.0)  RETFAC  =  0.0  001820 

WRITE(6,2100)  IFLT,ISEG,J4, CYC, A, IM.KEFF, KMXEFF, DDNRET, RETFAC  001830 

WRITE(1,2100)  IFLT, ISEG,J4, CYC, A, DLKEFF, KMXEFF, DDNRET, RETFAC 
RETURN  001840 

2000  F0RMAT(2X*FLT  MSN  LYR  CYCLES* 7X*A»4X  001850 

/*DELTA  K  K  MAX*5X»DA/DN  RETARD*  ) 

2100  F0RMAT(I5,2X,13,2X,I3, 1X,F10.1,1X,F9.6, IX,  001870 

/F7.2, 1X,F7.2, 1X,E10.3, 1X,F7.3)  001880 

END  001890 

SUBROUTINE  INPUTCICASE, ISPEC, IRSTRT)  001900 

C  001910 

C  READS  LABELED  SECTIONS  OF  DATA  DECK  IN  ANY  ORDER.  001920 

C  PRINTS  OUT  PROBLEM  AND  SOLUTION  DESCRIPTION.  001930 

C  001940 

COMMON/DATA/  EQN, NASA, J IPR, J2PR, J3PR, J4PR, J5PR, AZERO, AMAX, NZERO  001950 

INTEGER  EQN  001960 

REAL  NZERO  001970 

COMMON/RDATA/  MODEL, RETARD, PLSTRN, OVLD,SIGNAX, SIGMIN, ASUBP, SMALLM  001980 
INTEGER  RETARD, PLSTRN  001990 

COMMON/MDATA/  MATID< 18) , C, SMALLN,CARRAY(100) , SNARAY(IOO) , KSUBC,  002000 

/  KSUBQ,SIGMAY,DELKTH,RMULT,RCUT,OLMAX  002010 

REAL  KSUBC, KSUBQ  002020 

COMMON/LDATA/SMAX <  20,10),SMIN(  20, 10) , CYCLES!  20, 10) ,NLYRS(10) ,  002030 

/  NBLKS,IBLKS<  50  ),ISEGS(  50  ),NSE6S  002040 


COMMON/STEPS/ I SEG, J1 , J2, J3, J4, J5, ISTOP, NORTRD  002050 

COMMON/CORFAC/  ISURF, RATIO, PHI, THICK, IBETA(10),BETA(10),NPTS,  002060 

/  AOVERB(100),BTABLE(100),NPTS2,AOVRB2<100),  002070 

/  BTABL2(100),ASTART(10),ASTOP(10)  002080 

DIMENSION  CARD(1B),TITLE(18),SEGTTL(1B,50)  002090 

REAL  LODLABO)  002100 

REAL  LABEL (3)  002110 

REAL  SPCTRM(3),ENDLDS<3), IEND(3),LTITLE(3)  002120 

REAL  TA6(2)  002130 

REAL  EQNS(3),MATL(3),LIMITS(3),ANAL(3),L0ADS(3),END(3)  002140 

REAL  DADN(3),UALK(3)  002150 

REAL  EQUAT(3),DELk(3)  002160 

REAL  F0RMAN(3),PARIS(3),NASAL(3),BLAhtt<(3),SIGM0ID(3} 

REAL  SURF(3),RET(3),BETAL(3),BETAE(3),SIGMAS(3),DELTAS(3)  002180 

REAL  AMEANS(3),ENI>SPC(3),TRSHLD(3),LRSTRT(3),KPRTC3>  002190 

THE  NEXT  3  CARDS  ADDED  FOR  CL0S1«  002200 

CQMMQN/CLOS/CF, CCOEF , CFEXP, B, BOL, NSAT  002210 

REAL  NSAT  002220 

COMMON/CLOSIC/SC, SPEAK, PRVMX, APEAK, SCI , SC2, SC3, PRVMN  002230 

DATA  SPCTRM  /4HSPEC,4HTRUM, IH  /.ENDLDS  /3HEND,4HL0AD, IHS/  002240 

+  ,IEND  /3HEND,4HDATA,1H  /.LTITLE  /4HTITL, IHE, IH  /  002250 

DATA  EONS  /4HEQUA,4HTI0N, IH  /,MATL  /4HHATE,4HRIAL, IH  /  002260 

+  , LIMITS  /4HLIMI,2HTS, IH  /.ANAL  /4HANAL,4HYSIS, IH  /  002270 

*  , LOADS  /4HL0AD, IHS, IH  /,END  /3HEND,2*1H  /  002280 

DATA  DADN  /4HDA/D, IHN, IH  /,UALK  /4HUALK,2HER, IH  /  002290 

DATA  FORMAN  /4HF0RM, 2HAN, IH  /.PARIS  /4HPARI , IHS, IH  /  002300 

+  .NASAL  /4HNASA,2»1H  /.BLANK  /3*1H  /,SIGM0ID/4HSI6M,3H0ID, 1 
+H  / 

DATA  SURF  /4HSURF,3HACE, IH  /,RET  /4HRETA, 2HRD, IH  /  002320 

+  .BETAL  /4HBETA,2»1H  /.BETAE  /3HEND,2*1H  /  002330 

+  .SIGMAS  /4HMAX-,3HMIN, IH  /.DELTAS  /4HR-DE,4HLTA  , IH  /  002340 

DATA  AMEANS  /4HMEAN, 4H-ALT, IH  /.ENDSPC  /3HEND, 4HSPEC, IHR/  002350 

-►  .TRSHLD  /4HTHRE,4HSH0L,1HD  /.LRSTRT  /4HREST,3HART, IH  /  002360 

+  ,KPRT  /4HPRIN, IHT, IH  /  002370 

ISTOP  =  0  002380 

IF(ISPEC.GT.l)  GO  TO  6  002390 

INITIAL  CONDITIONS  FOR  CLOSUR  002400 

SC  =  0.  002410 

SCI  =0.  002420 

SC2  =  0.  002430 

SC3  =  0.  002440 

SPEAK  =  0.  002450 

PRVMX  =  0.  002460 

PRVMN  =  0.  002470 

APEAK  =  0.  002480 

END  OF  INITIAL  CONDITIONS  FOR  CLOSUR  002490 

IRSTRT  =  0  002500 

NORTRD  =  1  002510 

EON  =  1  002520 

NASA=0  002530 

NPTS=0  002540 

NPTS2*0  002550 


ISURF=0  0025&0 

M0DEL=0  002570 

JIPR  =  1  002580 

J2PR  =  0  002590 

J3PR  =  0  002600 

J4PR  =  0  002610 

J5PR  =  0  002620 

RATIO  =  1.0  002630 

RCUT  =  1.  002640 

KSUBC  =  68000.  002650 

KSUBQ  =  KSUBC  002660 

C  =  5.0E-13  002670 

SMALLN  =3.0  002680 

DELKTH  =  0.  002690 

RMULT  =  0.  002700 

DO  5  J=l,10  002710 

BETA(J)  =0.0  002720 

5  IBETA(J)  =0.0  002730 

6  REAO(S, 1000)  LABEL  002740 

IF (EOF (5))  9999,7  002750 

7  MRITE(6,2300)  I CASE, I SPEC  002760 

IF(IRSTRT.NE.O)  WRITE(6,3400)  002770 

60  TO  10  002780 

1  READ (5, 1000) LABEL  002790 

IF(E0F(5))  9998,10  002800 

10  IF(  LABEL(l)  .EQ.  LTITLEd)  .AND.  LABEL(2)  .EQ.  LTITLE(2)  002810 

+  .AND.  LABEL(3)  .EQ.  LTITLEO)  )  60  TO  100  002820 

IF(  LABELd)  .EQ.  KPRTd)  .AND.  LABEL(2)  .EQ.  KPRT(2)  002830 

+  .AND.  LABEL(3)  .EQ.  KPRT(3)  )  60  TO  150  002840 

1F(  LABELd)  .EQ.  EQNSd)  .AND.  LABEL(2)  .EQ.  EQNS(2)  002850 

+  .AND.  LABEL(3)  .EQ.  EQNS(3)  )  60  TO  200  002860 

IF(  LABELd)  .EQ.  MATLd)  .AND.  LABEL(2)  .EQ.  NATL(2)  002870 

.AND.  LABEL (3)  .EQ.  MATL(3)  )  60  TO  300  002880 

IF(  LABELd)  .EQ.  TRSHLDd)  .AND.  LABEL(2)  .EQ.  TRSHLD(2)  0028-?0 

♦  .AND.  LABEL (3)  .EQ.  TRSHLD(3)  )  60  TO  350  002900 

IF(  LABELd)  .EQ.  LIMITSd)  .AND.  LABEL(2)  .EQ.  LIt1ITS(2)  002910 

+  .AND.  LABEL<3)  .EQ.  LIMITSCS)  )  60  TO  400  002920 

IF(  LABELd)  .EQ.  ANALCl)  .AND.  LABEL(2)  .EQ.  ANAL(2)  002930 

.AND.  LABEL(3)  .EQ.  ANAL(3)  )  GO  TO  500  002940 

IF(  LABELd)  .EQ.  LOADSCl)  .AND.  LABEL(2)  .EQ.  L0ADS(2)  002950 

+  .AND.  LABEL(3)  .EQ.  L0ADS(3)  )  GO  TO  600  002960 

IF(  LABELd)  .EQ.  SPCTRMd)  .AND.  LABEL(2)  .EQ.  SPCTRM(2)  002970 

+  .AND.  LABEL<3)  .EQ.  SPCTRM(3)  )  60  TO  660  002980 

IF(  LABELd)  .EQ.  LRSTRTd)  .AND.  LABEL(2)  .EQ.  LRSTRT(2  )  002990 

+  .AND.  LABEL(3)  .EQ.  LRSTRT(3)  )  GO  TO  680  003000 

IF(  LABELd)  .EQ.  lENDCl)  .AND.  LABEL(2)  .EQ.  IEND(2)  003010 

+  .AND.  LABEL(3)  .EQ.  IEND(3)  )  GO  TO  700  003020 

WR I TE ( 6 , 9020 )  LABEL  003030 

ISTOP  =  1  003040 

GO  TO  1  003050 

100  READC5, 1010)  NTITLE  003060 

DO  no  1  =  1, NTITLE  003070 


-  94 


READ<5, 1000)  CARD  003080 

WRITE (6, 2005)  CARD  003090 

WRITE! 1,2005)  CARD 

no  CONTINUE  003100 

GO  TO  1  003110 

150  READ(5, 1010)  JIPR, J2PR, J3PR, J4PR, J5PR  003120 

GO  TO  1  003130 

200  RE AD (5, 1002)  EQUAT.DELK  003140 

IF(  EQUAT(l)  .EQ.  FORMAN(l)  .AND.  EQUAT(2)  .EQ.  F0RMAN(2)  003150 

+  .AND.  EQUAT<3)  .EQ.  F0RMAN(3)  )EQN  =  1  003160 

IF!  EQUAT!!)  .EQ.  PARIS!!)  .AND.  EQUAT!2)  .EQ.  PARISI2)  003170 

+  .AND.  EQUAT!3)  .EQ.  PARIS!3)  )EQN  =  2  003180 

IF!  EQUAT!1)  .EQ.  DADN!1)  .AND.  EQUAT!2)  .EQ.  DA0N!2)  003190 

.AND.  EQUAT!3)  .EQ.  DADN!3)  )EQN  =  3  003200 

IF!  EQUAT!1)  .EQ.  WALK!!)  .AND.  EQUAT!2)  .EQ.  WALK!2)  003210 

.AND.  EQUAT!3)  .EQ.  UALK!3>  )EQN  =  4  003220 

IF!  EQUAT!!)  .EQ.  SIGMOID!!)  .AND.  EQUAT!2)  .EQ.  SIGM0ID!2} 

+  .AND.  EQUAT!3)  .EQ.  SIGM0ID!3)  )EQN  =  5 

IF!  DELK!1)  .EQ.  NASAL!1)  .AND.  DELK!2)  .EQ.  NASAL!2)  003230 

+  .AND.  DELK!3)  .EQ.  NASAL!3)  )NASA  =  EON  003240 

IF!NASA.NE.O)  GO  TO  260  003250 

GO  TO  !2!0,220,230,240,250),EQN  003260 

210  WRITE!6,2010)  003270 

GO  TO  !  003280 

220  WR I TE ! 6 , 2020 )  003290 

GO  TO  !  003300 

230  WRITE! 6, 3030)  003310 

GO  TO  1  003320 

240  CONTINUE  003330 

WRITE! 6, 3050)  003340 

GO  TO  1  003350 

250  CONTINUE  003360 

WRITE!6,3070) 

GO  TO  1  003370 

260  GO  TO  !270, 280, 285, 290, 295),  NASA  003380 

270  WRITE! 6, 2030)  003390 

GO  TO  1  003400 

280  WRITE (6, 2040)  003410 

GO  TO  1  003420 

285  WRITE!6,3040)  003430 

GO  TO  1  003440 

290  CONTINUE  003450 

WRITE!6,3060)  003460 

GO  TO  1  003470 

295  CONTINUE  003480 

300  READ!5,!000)  MATID  003490 

WRITE! 6, 2050)  MATID  003500 

CALL  CNDIN!EQN)  003510 

GO  TO  1  003520 

350  READ!5, 1020)  DELKTH,RMULT  003530 

WRITE! 6, 2260)  DELKTH,RMULT  003540 

GO  TO  1  003550 


400  READ(5, 1020)  AZER0,AMAX,NZER0,RCUT  003560 

IF(RCUT.EQ.O.)  RCUT  =1.  003570 

IFCAMAX.EQ.O. )  GO  TO  410  003580 

WRITE(6,2070)  AZERO.AMAX.NZERQ  003590 

MRITE(6,2075)  RCUT  003600 

GO  TO  1  003610 

410  AMAX=1.0E+50  003620 

URITG(6,2080)  AZERO.NZERO  003630 

MRITE(6,2075)  RCUT  003640 

GO  TO  1  003650 

500  READ(5, 1070)  LABEL, CON 1,C0N2,C0N3,CQN4, CONS, C0N6  003660 

IF(  LABEL(l)  .EQ.  SURF(l)  .AND.  LABEL<2)  .EQ.  SURF(2)  003670 

+  .AND.  LABEL (3)  -EQ.  SURF(3)  )  GO  TO  510  003680 

1F(  LABEL(l)  .EQ.  RET(l)  .AND.  LABEL(2)  .EQ.  RET(2)  003690 

+  .AND.  LABEL (3)  .EQ.  RET (3)  )  GO  TO  520  003700 

IF(  LABEL(l)  .EQ.  BETAL(l)  .AND.  LABEL(2)  .EQ.  BETAL(2)  003710 

+  .AND.  LABEL (3)  .EQ.  BETAL(3)  )  GO  TO  530  003720 

IF(  LABEL(l)  .EQ.  BETAE(l)  .AND.  LABEL(2)  .EQ.  BETAE(2)  003730 

+  .AND.  LABEL<3)  .EQ.  BETAE(3)  )  GO  TO  1  003740 

WRITE (6, 9020)  LABEL  003750 

I STOP  =  1  003760 

GO  TO  1  003770 

510  ISURF=1  003780 

CZERQ=C0N1  003790 

TH1CK=C0N2  003800 

RATI0=AZER0/(2.*CZER0)  003810 

SMALLK  =  SQRT<(CZER0**2  -  AZER0**2)  /  CZER01:»2)  003820 

CKSQD  *  1.0  -  SMALLK*l:2  003830 

CALL  CELI2(PHI , SMALLK, 1 . 0, CKSQD, lER)  003840 

WRITE(6,2090)  RATIO, THICK, PHI  003850 

IFCIER.EQ.O)  GO  TO  500  003860 

WRITE (6, 9030)  003870 

I STOP  =  1  003880 

GO  TO  500  003890 

520  M0DEL=C0Nl+0.5  003900 

PLSTRN*C0N2+1 . 5  003910 

NORTRD  =  C0N4  +0.5  003920 

OVLD  =  CONS  003930 

ASUBP  =  C0N6  003940 

IF (NORTRD. NE.O)  WRITE(6,2170)  003950 

GO  TO  (521, 522, 523,524, 525), MODEL  003960 

521  SMALLM»C0N3  003970 

WRITE(6,2100)  SMALLM  003980 

GO  TO  526  003990 

522  WRITE(6,2105)  004000 

OLMAX  =  C0N3  004010 

IF(OLMAX.NE.O.)  WRITE(6,2106)  OLMAX  004020 

GO  TO  526  004030 

523  CONTINUE  004040 

WRITE(6,2107)  004050 

READ(5,1020)  CF,CCOEF,CFEXP,B,BOL,NSAT  004060 

WRITE(6,2108)  CCOEF,CF,CCOEF,CFEXP,B,BOL,NSAT  004070 


RETARD  =  1 
GO  TO  526 

524  CONTINUE 

525  WRITE (6, 3000) 

M0DEL=2 

526  GO  TO  (527,528),PLSTRN 

527  WRITEC6,2180) 

PLSTRN  =  PLSTRN  -  1 
GO  TO  500 

528  WRITE(6,2190) 

PLSTRN  =  PLSTRN  -  1 
GO  TO  500 

530  I  =  CONI  +0.5 
IBETA(I)=I 

GO  T0(531 , 532,533,534,535,536, 537, 538, 539) , I 

531  BETA<I)=C0N2 
ASTART(I)=C0N3 
AST0P(I)=C0N4 

IF(ASTART(I).EQ.O.)  ASTART(l)  =  A2ER0 
IF(ASTQP(I).EQ.O.)  ASTOP(I)  =  1.E50 
WRITE<6,2110)  BETA(1),ASTART(I),AST0P(I) 

GO  TO  500 

532  BETA(I)=C0N2 
ASTART(I)=C0N3 
AST0P(I)=C0N4 

IF(ASTART(I).EQ.O.)  ASTART(I)  =  AZERO 
IF<ASTOP(I).EQ.O. )  ASTOP(I)  =  I.E50 
WR I TE ( 6 , 2 1 20)  BETA ( 2 ) , ASTART ( I ) , ASTOP ( I ) 

GO  TO  500 

533  BETA(I)=C0N3 
NPTS=C0N2+0.5 
ASTART(I)=C0N4 
AST0P(I)=C0N5 

IF(ASTART(I).Ea.O.)  ASTART(I)  =  AZERO 
IFCASTOPCn.EQ.O.)  ASTOP(I)  =  1.E50 
IFCNPTS.LE.O.OR.NPTS.GT. 100)  GO  TO  540 
READ(5,1000)  CARD 

READ(5, 1040)  (AOVERB(J) , BTABLE(J) , J=1,NPTS) 

WF. .£(6,2130)  BETA(I),ASTART(I),ASTOP(I),CARD, 
1  (AOVERB(J),BTABLE(J),J=l,NPTS) 

GO  TO  500 

534  BETA(I)=C0N3 
NPTS2=C0N2+0. 5 
ASTART(1)=C0N4 
AST0P<I)=C0N5 

IF(ASTART(I).EQ.O.)  ASTART(I)  =  AZERO 
IF(ASTOP(I).EQ.O. )  ASTOP(I)  =  1.E50 
IF(NPTS2.LE.0,0R.NPTS2.QT. 100)  GO  TO  540 
READ(5,1000)  CARD 

READ (5,1040)  ( A0VRB2 ( J ) , BTABL2 ( J ) , J= 1 , NPTS2 ) 
WRITE(6, 2140)  BETA( I ) , ASTART ( I ) , ASTOP( 1 ) , CARD, 
1  (A0VRB2(J),BTABL2(J),J=1,NPTS2) 
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004130 

004140 
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004200 

004210 

004220 

004230 

004240 

004250 

004260 

004270 

004280 

004290 

004300 

004310 

004320 

004330 

004340 

004350 

004360 

004370 

004380 

004390 

004400 

004410 

004420 

004430 

004440 

004450 
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004470 

004480 

004490 

004500 
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004520 

004530 

004540 

004550 

004560 

004570 

004580 

004590 


GO  TO  500 

535  CONTINUE 
BETA(I)  =  C0N2 
ASTART(I)  =  C0N3 
ASTOP(I)  =  C0N4 

IF(ASTART(I).EQ.O.O)  ASTART(I)  =  AZERO 
IFCASTOPCD.EQ.O.O)  ASTOP(I)  =  1.0E50 
WR I TE ( 6 , 2200 )  BETA ( I ) , ASTART ( I ) , ASTOP ( I ) 

GO  TO  500 

536  CONTINUE 
BETACI)  =  C0N2 
ASTART (I)  =  C0N3 
ASTOP(I)  =  C0N4 

IF < ASTART ( I). EQ. 0.0)  ASTART (I)  =  AZERO 
IFCASTOPCD.EQ.O.O)  ASTOPCI)  =  1.0E50 
WR I TE  C  6 , 2205 )  BETA ( I ) , ASTART ( I ) , ASTOP ( I > 

GO  TO  500 

537  CONTINUE 

BETA(I)  =  C0N2 
THICK  =  C0N3 
ASTART(I)  =  C0N4 
ASTOP(I)  =  CaN5 

IF  (  ASTARTCI)  .EQ.  0.0  )  ASTARTd)  =  AZERO 
IF  (  ASTOPd)  .EQ.  0.0  )  ASTOP(I)  =  1.E50 
WRITE (6, 2206)  BETA ( I ), THICK, ASTART! I ) , ASTOP! I ) 
GO  TO  500 

538  CONTINUE 

BETA! I)  =  C0N2 
THICK  =  C0N3 
ASTART!!)  =  C0N4 
ASTOP! I)  =  CONS 

IF  !  ASTART (I)  .EQ.  0.0  )  ASTART!!)  =  AZERO 
IF  !  ASTOP!I)  .EQ.  0.0  )  ASTOP!I)  =  1.E50 
WR I TE ! 6 , 2207 )  BE TA ! I ) , TH I CK , ASTART ! I ) , ASTOP ! I ) 
GO  TO  500 

533  CONTINUE 

BETA! I)  =  C0N2 
THICK  =  CONS 

ASTART! I)  =  C0N4 
ASTOP! I)  =  CONS 

IF  !  AST ART! I)  .EQ.  0.0  )  ASTART! I)  =  AZERO 
IF  !ASTOP!I)  .EQ.  0.0  )  ASTOPII)  =  1.E50 
WRITE!6, 2208)  BETA! I ) , THICK, ASTART! I ) , ASTOP! I ) 
GO  TO  500 

540  WRITE!6,9000)  IBETA!I) 

IBETA!I)=0 
GO  TO  500 

600  READ!5,1<^0)  NBLKS,LPRT, !TITLE!I) , 1=1, 18) 
WRITE!6,,^00)!TITLE!I),I=l,l8),NBLKS 
IF!LPRT.Cf.O)  GO  TO  1 
ISEG  =  0 
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004760 

004770 

004780 

004790 
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004830 

004840 

004850 

004860 

004870 

004880 

004890 
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004960 

004970 
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• 

£05  ISE6  =  ISE6  +  1 

005040 

READ (5, 1002)  LODLAB,  (SEGTTLd,  ISEG),  1  =  1, 18) 

005050 

• 

LYR=1 

005060 

_J 

IF(  LODLAB(l)  .EQ.  SIGMASd)  .AND.  L0DLAB(2)  .EQ.  SIGMAS(2) 

005070 

+  .AND.  LODLABO)  .EQ.  SIGMASO)  )  60  TO  610005080 

•  .  » 

• 

IF(  LODLABCl)  .EQ.  DELTASd)  .AND.  L0DLAB(2)  .EQ.  DELTAS(2) 

005090 

•’  y..' 

+  .AND.  LODLABO)  .EQ.  DELTASO)  )  GO  TO  620 

005100 

IF(  LODLABCl)  .EQ.  AMEANSCl)  .AND.  L0DLABC2)  .EQ.  AMEANSCZ) 

005110 

• 

+  .AND.  LQDLABO)  .EQ.  AMEANS(3>  )  GO  TO  630 

005120 

i 

IFC  LODLABCl)  .EQ.  ENDLDSCl)  .AND.  LQDLABC2)  .EQ.  ENDLDSC2) 

005130 

•/ 

+  .AND.  LQ0LABC3)  .EQ.  ENDLDS(3)  )  GO  TO  662 

005140 

WR I TEC 6, 9020)  LODLAB 

005150 

I STOP  =  1 

005160 

GO  TO  1 

005170 

1.  -- 

• 

c 

005180 

• 

4 

c 

LOAD  SPECTRUM  INPUT  AS  MAX  AND  MIN  STRESSES 

005190 

c 

005200 

610  READC5, 1030)  TAQ,SMAXCLYR, ISE6) ,SMINCLYR, I SEG ), CYCLES CLYR, ISEG) 

005210 

IFC  TA6C1)  .EQ.  ENDCl)  .AND.  TA6(2)  .EQ.  ENDC2)  )  GO  TO  650 

005220 

c 

005230 

c 

PRESENT  PROGRAM  DOES  NOT  CONSIDER  COMPRESSIVE  LOADS 

005240 

J 

c 

005250 

c 

EXCEPT  FOR  CLOSURE  MODEL  WHICH  DOES 

005260 

c 

005270 

IF  C  MODEL  .EQ.  3  )  GO  TO  618 

005280 

IFCSMAXCLYR, ISEG) .LT.O. )  SMAXCLYR, ISEG)=0. 

005290 

IFCSMINCLYR,ISEe).LT.O. )  SMINCLYR, ISEG)  =  0. 

005300 

• 

( 

w 

618  LYR=LYR+1 

005310 

GO  TO  610 

005320 

620  IFCLPRT.EQ.O) 

005330 

1WRITEC6,3010)  ISEG, CSEGTTLCI, ISEG), 1=1, 18) 

005340 

IFCLPRT.EQ.O)  WRITEC6, 2210) 

005350 

625  READC5, 1030)  TAG, DELSIQ, R, CYCLESCLYR, ISEG) 

005360 

A 

a 

o 

IFC  TAGCl)  .EQ.  ENDCl)  .AND.  TA6C2)  .EQ.  ENDC2)  )  GO  TO  650 

005370 

*  .  • 

IFCLPRT.EQ.O) 

005380 

1WRITEC6,2220)  L YR, T AG, DELS I Q,R, CYCLES CLYR, ISEG) 

005390 

/*y 

c 

005400 

c 

PRESENT  PROGRAM  DOES  NOT  CONSIDER  COMPRESSIVE  LOADS 

005410 

c 

005420 

*■  “  * 

1  \ 

IFCDELSIG.LT.O.)  DELSIG=0. 

005430 

■L— 

IFCR.LT.O. )R=0. 

005440 

SMAXCLYR, ISEQ)=DELSIG/C1.-R) 

005450 

SMINCLYR, ISEG)=SMAXCLYR, ISEG)-DELSI6 

005460 

■  »  * 

LYR=LYR+1 

005470 

GO  TO  625 

005480 

1  t 

630  IFCLPRT.EQ.O) 

005490 

• 

1WRITE(6,3010)  ISEG, CSEGTTLCI, ISEG), 1=1, 18) 

005500 

IFCLPRT.EQ.O)  WRITEC6,2230) 

005510 

•  * , 

*.y 

635  READC5,1030)  TAG, SMEAN, SALT, CYCLESCLYR, ISEG) 

005520 

•  * 

IFC  TAGCl)  .EQ.  ENDCl)  .AND.  ENDC2)  .EQ.  TAGC2))  GO  TO  650 

005530 

.’’y 

IFCLPRT.EQ.O) 

005540 

*'*/•*-' 

1 

1WRITEC6,2240)  LYR, TAG, SMEAN, SALT, CYCLES CLYR, ISEG) 

005550 

IL.^ 

-  99  - 

o  o  o 


C  005560 

C  PRESENT  PROGRAM  DOES  NOT  CONSIDER  COMPRESSIVE  LOADS  005570 

C  005580 

SMAX ( LYR , I SEG ) =SMEAN+SALT  005590 

IF(SMAX(LYR, ISEQ).LT.O.)  SMAX (LYR, ISEG)=0.  005600 

SMIN(LYR,ISEG)=SMEAN-SALT  005610 

IF(SMIN(LYR, ISEG).LT.O.)  SMIN<LYR, ISEG)=0.  005620 

LYR=LYR+1  005630 

GO  TO  635  005640 

650  NLYRS(ISEG)  =  LYR-1  005650 

ijO  to  605  005660 

005670 

READ  IN  COMPOSITION  OF  SPECTRUM  005680 

005690 

660  READ(5, 1010)NSEGS, IPRT  005700 

READ(5, 1060)  (IBLKSd),  ISEGS(I),  I*1,NSEGS)  005710 

IF(IPRT.EQ.O)  WRITE(6,2250)  NSEGS  005720 

IFLTS  =  0  005730 

DO  661  1=1, NSEGS  005740 

IFLTS  =  IFLTS  +  IBLKS(I)  005750 

IF(IPRT.EQ.O)  WRITE(6,2255)  I, IBLKS(I), ISEGS(I), IFLTS  005760 

661  CONTINUE  005770 

GO  TO  1  005780 

662  NNSEG  =  I  SEG  -  1  005790 

IF(LPRT.NE,0)  GO  TO  1  005800 

DO  670  ISEG  =1, NNSEG  005810 

WRITE(6,2270)  ISEG, (SEGTTLCI, ISEG), 1=1, 18)  005820 

LYRS  =  NLYRSdSEQ)  005830 

TOTCYC  =0.0  005840 

DO  665  J2=1,LYRS  005850 

TOTCYC  »  TOTCYC  +  CYCLES( J2, ISEG)  005860 

WRITE<6,2280)  J2, SMAX <J2, ISEG), SMIN<J2, ISEG), CYCLES! J2, ISEG),  005870 

1  TOTCYC  005880 

665  CONTINUE  005890 

670  CONTINUE  005900 

GO  TO  1  005910 

680  READ(5, 1060)  IRSTRT  005920 

I CHECK  =  IRSTRT  +  1  005930 

60  TO  (1,685,685,690,695,695), I CHECK  005940 

685  WRITE(6,2150)  005950 

GO  TO  1  005960 

690  WRITE(6,2160)  005970 

GO  TO  1  005980 

695  WRITE<6,2165)  005990 

GO  TO  1  006000 

700  WRITE(6,2015)  (TITLE! I ), 1=1 , 18)  006010 

RETURN  006020 

9998  IFdRSTRT  .NE.  0)  WRITE(6,9010)  006030 

9999  STOP  006040 

C  006050 

1000  FORMAT (20A4)  006060 

1002  FORMAT (2A4,A2,17A4,A2)  006070 


1010  FORMAT (1615) 

1020  FORMAT (8E 10.0) 

1030  FORMAT <A4,A1,7E10.0) 

1040  FORMAT (2E 10.0) 

1050  F0RMAT(2I5,17A4,A2) 

1060  F0RMAT(2I5) 

1070  F0RMAT<2A4,A2,7E10.0) 

2000  FORMATdHO,  1X,17A4,A2//  5X,I6,19H  BLOCKS  IN  SPECTRUM) 
2005  FORMAT(  2X,17A4,A2) 


006080 

006030 

006100 

006110 

006120 

006130 

006140 

006150 

006160 


2010  FORMATdHO,  1X,50HCRACK  PROPAGATION  ANALYSIS  USING  FORMAN'S  EQUATI006170 
ION  /  4X,44HDA/DN=C* (DELTA  K)**N/(d-R)»KSUBC-DELTA  K)  006180 

2  /  5X,51HWHERE  K  IS  OF  THE  FORM  _  K=SIQMA»SQRT(PITA)*BETA)  006130 

2015  F0RMAT(//lX,29dH*),12HEND  OF  INPUT, 29dH»)/////lHl,  IX, 70dH$)/  006200 

/lX,26dH»),  18HCRACKS  IV  ANALYSrS,26dH»)/lX,  17A4,A2/  006210 

/lX,70dH$)//)  006220 

2020  FORMATdHO,  1X,50HCRACK  PROPAGATION  ANALYSIS  USING  PARIS'  EQUATI0N006230 

1  /16X,21HDA/DN=C»(DELTA  K)**N  /  5X,51HWHERE  K  IS  OF  THE  FORM  .  006240 

2..  K=SIGMA«SQRT(PI«A)*BETA  )  006250 

2030  FORMATdHO,  1X,50HCRACK  PROPAGATION  ANALYSIS  USING  FORMAN'S  EQUATI006260 
ION  /  4X,44HDA/DN=C»( DELTA  K)M:N/(d-R)»KSUBC-DELTA  K)  006270 

2  /  5X,S1HUHERE  K  IS  OF  THE  FORM  ....  K-SI6MA»SQRT(A)4rBETA  )  006280 
2040  FORMATdHO,  1X,50HCRACK  PROPAGATION  ANALYSIS  USING  PARIS'  EQUATI0N006290 

1  /16X,21HDA/DN=C*(DELTA  K)**N  /  5X,51HUHERE  K  IS  OF  THE  FORM  .  006300 

2..  K=SIGMA.l!SQRT(A)*BETA  )  006310 

2050  FORMATdHO,  1X,20A4)  006320 

2070  FORMATdHO,  1X,27HINITIAL  HALF  CRACK  LENGTH  =  E16.a/2X,35HMAXIMUM  00633«5 
IHALF  CRACK  LENGTH  ALLOWED  =  E16.8/1H0,  1X,22HINITIAL  CYCLE  NUMBER  006340 

2  ,F11.2  )  006350 

2075  FORMATdHO,  IX,  11  HR  CUTOFF  =  ,F6.3)  006360 

2080  FORMATdHO,.  IX, 27HINITIAL  HALF  CRACK  LENGTH  *  E16.8  /IHO,  IX, 22HIN006370 

IITIAL  CYCLE  NUMBER  =  , FI 1.2  )  006380 

2030  FORMATdHO,  1 X , 34HSURFACE  FLAW  iWALYSIS  WITH  A/2C  OF  F5.2  /  5X,21H006330 
IMATERIAL  THICKNESS  IS  ,F8.5/  5X, 13HSHAPE  FACTOR  PHI  =  ,F8.5)  006400 

2100  FORMATdHO,  1X,41HHHEELER'S  RETARDATION  MODEL  WITH  SMALLM  =  ,F6. 3)006410 

2105  FORMATdHO,  1X,*W I LLENBORG  RETARDATION  MODEL*)  006420 

2106  FORMATdHO, IX, »GALLAGHER-MODIFIED  WILLENBORG  RETARDATION  MODEL  */  006430 

2X*WHERE. ..*//  5X,27HPHI  =  d-THRESHOLD/K  MAX)/(,F5.2,3H-1) )  006440 

C  FORMATS  FOR  CLOSUR  INPUT/OUTPUT  DATA  006450 

2107  FORMATdHO,  IX,  18HCL0SURE  MODEL  WITH)  006460 

2108  FORMATdHO,  IX, I6HCL0SURE  FACTOR  =,F6.4,3H+  (,F6.4, 1H-,F6.4,  006470 

1  10H)*d.-R)«:,F6.4,/  5X,33HEXP0NENT  FOR  DECREASING  CLOSURE  =,006480 

/  F6.3,/  006430 

2  5X,34HEFFECTIVENESS  AFTER  ONE  OVERLOAD  =,F6.4/  006500 

3  5X,36HNUMBER  OF  OVERLOADS  FOR  SATURATION  =,F6.0)  006510 

2110  FORMATdHO,  IX, 40HC0RRECTI0N  FACTOR  BETAd)  IS  A  CONSTANT  /  5X,9H006520 

IBETACl)  =  ,E16.8/  5X, 16HAPPLIED  FROM  A=  ,E16.8,7H  TO  A  =  ,E16.8)  006530 

2120  FORMATdHO,  1X,54HC0RRECTI0N  FACTOR  BETA(2)  IS  FINITE  WIDTH  C0RREC006540 
ITION  /5X28HBETA(2)  =  SaRT(SEC<PI*A/B) )  /  006550 

25X*WHERE  THE  EFFECTIVE  PLATE  WIDTH  W  =  TE16.8/  006560 

35X*APPLIED  FROM  A  =  *El6.a,»  TO  A  =  *616.8)  006570 

2130  FORMATdHO,  1X,54HC0RRECTI0N  FACTOR  BETA (3)  IS  A  TABULAR  FUNCTION  006580 
lOF  A/L  /  5X,10HWHERE  L  =  ,E16.8  /  5X, 16HAPPLIED  FROM  A  =  ,E16. 8, 7006530 


rj  ro  h  j  M 


2H  TO  A  =  ,E16.8//  lX,iaA4//  aX,3HA/L,  16X,7HBETA(:3)/(  5X,E15.a,5X,  006600 
3E15.8))  006610 

2140  FORMAT (1 HO,  5X, 55HC0RRECTI0N  FACTOR  BETA (4)  IS  A  TABULAR  FUNCTION  006620 
lOF  A/Ll  /  5X,10HWHERE  L  =  E16.8/  5X, 16HAPPLIED  FROM  A  =  ,E16.a,7H006630 
2  TO  A  =  ,E16.8//  1X,18A4//  aX,3HA/L, 16X,7HBETA(4)/C  5X, E15. 8,5X,E1006640 
35.8))  006650 

150  FORMATCIHO,  1X,35HTHIS  CASE  WILL  BE  RESTARTED  ON-LINE)  006660 

160  FORMATCIHO,  1X,36HTHIS  CASE  WILL  BE  RESTARTED  OFF-LINE  /  006670 

1  48H  RESTART  DATA  WILL  BE  WRITTEN  ON  LOGICAL  UNIT  7)  006680 

165  FORMATCIHO, 1X,36HTH IS  CASE  WILL  BE  RESTARTED  OFF-LINE  /  006690 

1  4aH  RESTART  DATA  WILL  BE  READ  FROM  LOGICAL  UNIT  7  )  006700 

170  FORMATCIHO,  IX, 40HAUT0MATIC  UNRETARDED  SOLUTION  SUPPRESSED)  006710 

180  FORMATCIHO,  1X,41HPLANE  STRESS  YIELD  ZONE  CONDITION  ASSUMED  )  006720 

2190  FORMATCIHO,  1X,41HPLANE  STRAIN  YIELD  ZONE  CONDITION  ASSUMED  )  006730 

2200  FORMATCIHO,  IX, 71HC0RRECTI0N  FACTOR  BETAC5)  IS  UNIAXIAL  BOWIE  S0LU006740 

ITION  FOR  A  SINGLE  CRACK/  5X,35HFR0M  A  CIRCULAR  HOLE  OF  RADIUS  R  =  006750 
2E16.8  /  5X, 16HAPPLIED  FROM  A  =  ,E16.8,8H  TO  A  =  ,E16.B)  006760 

2205  FORMATCIHO,  1X,68HC0RRECTI0N  FACTOR  BETAC6)  IS  UNIAXIAL  BOWIE  S0LU0O6770 

ITION  FOR  TWO  CRACKS  /5X,35HFR0M  A  CIRCULAR  HOLE  OF  RADIUS  R  =  ,  006780 

2E16.8  /  5X, 17HAPPLIED  FROM  A  =  ,E16.B,BH  TO  A  =  ,E16.8)  006790 

2206  FORMATCIHO,  1X,S8HC0RRECTI0N  FACTOR  BETAC7)  IS  ASTM  COMPACT  TENSI0006800 

IN  SPECIMEN/  5X,15HWITH  A  WIDTH  0F,E16.B/  006810 

/  5X,16HAND  THICKNESS  0F,E16.5/  006820 

2  5X,17HAPPLIED  FROM  A  =  ,E16.8,BH  TO  A  =  ,E16.B)  006830 

2207  FORMATCIHO,  1X,61HC0RRECTI0N  FACTOR  BETACB)  IS  GRUMMAN  COMPACT  TEN006840 

ISION  SPECIMEN/  5X, 15HWITH  A  WIDTH  0F,E16.8/  006850 

/  5X,16HAND  THICKNESS  0F,E16.5/  006360 

2  5X, 17HAPPLIED  FROM  A  =  ,E16.8,BH  TO  A  =  ,E16.8)  006870 

2208  FORMATCIHO,  IX, 71HC0RRECTI0N  FACTOR  BETAC9)  IS  ASTM  E647-83  FOR  A 


ICOMPACT  TENSION  SPEIMEN/  5X,15HWITH  A  WIDTH  OF, £16.8/ 

/  5X,16HAND  THICKNESS  0F,E16.5/ 

2  5X,17HAPPLIED  FROM  A  =  ,E16.8,8H  TO  A  =  ,E16.8) 

2210  FORMAT  C2X*LAYER.T4X*LAB£L*7X«DELTA*9X.tR:C:4X*CYCLES  PER*/  006880 

/23X*SIGMA*16X*LAYER*//)  006890 

2220  FORMAT  C 3X , 1 2 , 5X , 2A4 , 1 X , E 1 3 . 5 , 1 X , F7. 4 , 2X , F 1 0 . 2 )  006900 

2230  FORMAT C2X*LAYER*4X*LABEL*8X*MEAN*6X*ALTERNAT I NG*3X  006910 

/*CYCLES  PER*/23X*STRESS*8X*STRESS*7X*LAYER*// )  006920 

2240  FaF:MATC3X,I2,5X,2A4,lX,E13.5,lX,E13.5,2X,F10.2)  006930 

2250  FORMATCIHI, IX, I6,52H  SEGMENT  SPECTRUM  APPLIED  IN  THE  FOLLOWING  3EQ006940 
/UENCE  //1X*SEGMENT»11X»FLIGHTS  PER*6X*MISSI0N*11X*CUMULATIVE*/  006950 
/21X*MISSiaN*28X*FLIGHTS*//)  006960 

2255  F0RMATC3X,I2, 17X, I4,12X,I2,17X,I5)  006970 

2260  FORMATCIHO,  IX, 19HTHRESH0LD  DELTA  K  =,E17.8,7H  Cl.0-C,F6.3,4HJ*R) >006980 
2270  FORMATCIHI,  70C1H*)/  5X,34HSTRESS  SPECTRUM  FOR  MISSION  NUMBER, 13/  006990 
/2X, 17A4,A2/1X,70C1H*)//2X*LAYER*5X*MAXIMUM*7X*MINIMUM*7X  007000 

/♦:CYCLES*6X*CUMULATIVE*/13X*STRESS*8X*STRESS*BX*PER*8X  007010 

/tCYCLES  PER*/40X*LAYER*9X*MISSI0N*//)  007020 

2280  F0P.MATC3X,  I2,4X,4CF12.3,2X))  007030 

2900  FORMATCIHI, 70C1H*)/26X,5HCASE  ,I2,5X,4HRUN  , I2/1X, 70C IH*) )  007040 

3000  FORMATCIHO,  70C1H*)/  5X, 66HINACTIVE  RETARDATION  MODEL  CHOSEN. EFFEC007050 
ITIVE  STRESS  MODEL  ASSUMED.  /IX,  70C1H*))  007060 

3010  FORMATCIHI,  5X,33HINPUT  SPECTRUM  FOR  MISSION  NUMBER, 14/  5X, 18A4)  007070 


3030  FORMAT (IHO,  1X,65HCRACK  PROPAGATION  ANALYSIS  USING  DIRECT  INPUT  0F007080 

1  DA/DN  VS. DELTA  K  /  5X,51H  WHERE  K  IS  OF  THE  FORM  _ SIGMA*SQRT(F007080 

2I*A)*BETA  )  007100 

3040  FORMAT! IHO,  1X,65HCRACK  PROPAGATION  ANALYSIS  USING  DIRECT  INPUT  OF007110 

1  DA/DN  VS. DELTA  K  /  5X,51HWHERE  K  IS  OF  THE  FORM  -  K=SIGMA*SQR007120 

2T(A)*BETA  )  007130 

3050  FORMATCIHO,  1X,50HCRACK  PROPAGATION  ANALYSIS  USING  WALKER'S  EQUATI007140 
ION/  7X,35HDA/DN=C*( DELTA  K/((1-R)»»(l-M)) )»*N  /  5X,52HWHERE  K  IS  007150 

20F  THE  FORM _  K=SIGMA*SQRT(PI*A).1:BETA  )  007160 

3060  FORMAT! IHO,  1X,50HCRACK  PROPAGATION  ANALYSIS  USING  WALKER'S  EQUATI007170 
ION/  7X,35HDA/DN=C*!DELTA  K/! ! 1-R)**! 1-M) ) )*»N  /  5X,51HWHERE  K  IS  0007180 

2F  THE  FORM _ K=SIGMA*SQRT!A)*rcTA  )  0071'30 

3070  FORMATCIHO,  1X,50HCRACK  PROPAGATION  ANALYSIS  USING  SIGMOIDAL  EON. 

1  /  7X,60HDA/DN=EXP!B)  CDELTA  K/DELTA  K*)*«P  !LN!DELTA  K/DELTA  K*) 

2)**Q/12X,26H  CLNCDELTA  KC/DELTA  K))**D  /  5X,51HWHERE  K  IS  OF  THE  F... 
30RM  ....  K=SIQMA»SQRT!PI:l:A);tBETA  ) 

3400  FORMATCIHO, 70!1H*)/21X,40HRERUN  OF  CASE  WITH  THE  FOLLOWING  CHANGES007200 
/  /  1X,70!1H*))  007210 

9000  FORMATCIHO,  70C1H*)/  5X,51HNUMBER  OF  POINTS  IN  TABULAR  CORRECTION  007220 
IFACTOR  BETACIl,t)  EXCEEDS  100.*/*  CORRECTION  FACTOR  WILL  BE  IGN0RE007230 
2D.*/1X,  70! IH*)  )  007240 

9010  FORMATCIHO,  70C1H*)/  5X,38HERR0R  IN  DECK  SETUP. E-O-F  ENCOUNTERED.  007250 
1  /IX,  70C1H*))  007260 

9020  FORMATCIHO,  70C1H*)/  5X,48H INCORRECT  LABEL  CARD  ENCOUNTERED. LABEL  007270 
♦READ  WAS  ,2A4,A2,1H*/  1X,64HEXECUTI0N  SUPPRESSED. PROGRAM  WILL  C0MP007280 
2LETE  INPUT  DATA  PROCESSING/ IX,  70C1H*))  007290 

9030  FORMATCIHO,  70C1H*)/  1X,71HERR0R  IN  CALCULATING  PHI. PROGRAM  REQUIR007300 
♦ES  (AZER0/2CCZER0))SQD  .LE.  0.5/  1X,65HEXECUTI0N  SUPPRESSED.  PR0G007310 
♦RAM  WILL  COMPLETE  INPUT  DATA  PROCESSING/ IX,  70(1H*))  007320 

end  007330 

SUBROUTINE  DELTACA,DELTAK,KMAX,R)  007340 

COMMON/RDATA/  MODEL, RETARD, PLSTRN,OVLD,SIGMAX, SIGMIN, ASUBP, SMALLM  007350 


INTEGER  RETARD, PLSTRN 

COMMON/MDATA/  MATIDC 18) ,C,SMALLN,CARRAYC 100) , SNARAYC 100) ,KSUBC, 
/  KSUBQ, SIGMAY, DELKTH, RMULT, RCUT, OLMAX 

REAL  KSUBC, KSUBQ 

COMMON/CORF AC/  ISURF, RATIO, PHI , THICK, IBETAC 10) , BETA! 10) , NPTS, 

/  AOVERB  C 1 00 ) , STABLE  C 1 00 ) , NPTS2 , A0VRB2  C 1 00 ) , 

/  BTABL2(100),ASTART<10),AST0PC10) 

REAL  KMAX,KMIN 
CALL  KCSIGM1N,A,KMIN) 

CALL  KCSIGMAX,A,KMAX) 

R  =  0.0 

IFC  KMAX  .NE.  0.0  )  R  =  KMIN/KMAX 

DELTAK=KMAX-KMIN 

RETURN 

END 

REAL  FUNCTION  TRP2CT, X, Y,M) 

DIMENSION  TC100),ZC4),D(6) 


007360 

007370 

007380 

0073'30 

007400 

007410 

007420 

007430 

007440 

007450 

007460 

007470 

007480 

007490 

1X17500 

007510 

007520 

007530 

007540 

007550 


- v.v 


-  103  - 


I 


o  o  o 


X1=X  007560 

Y1=Y  007570 

I  a  T(l)  /  1000.  +  1.  007580 

J  =  AM0D(T(1),1000.)  +  1.  007590 

L=J*M  007600 

I1=J*3+1  007610 

I2=I:l:J  007620 

M1=M  007620 

DO  10  K=I1, I2,L  007640 

IF(X1-T(K))  20,20,10  007650 

10  CONTINUE  007660 

K=I2+1-J  007670 

20  DO  30  L=4,J,M1  007680 

IF  (Y1-T<L))  40,40,30  007690 

30  CONTINUE  007700 

L=J  007710 

40  L1=L1+1  007720 

DO  50  MN=1,3  007730 

N=L+MN-3  007740 

Nl=K+(J1:(Ll-3))+N-l  007750 

D(MN)=T(N)  007760 

50  D<MN+3)=T<N1)  007770 

60  Z(Ll)=D(4)+(Yl-D(l))1:(<D(5)-D<4))/(D(2)-D(l))+(  007780 

1Y1-D(2))/(D(3)-D(1))*((D(6)-D(5))/(D(3)-D(2))  007790 

2-<D(5)-D(4))/(D(2)-D(l))))  007800 

IF  (Ll-3)40,70,90  007810 

70  DO  80  MN=1,3  007820 

D(MN+3)=Z(:MN)  007830 

Nl=K+(J:|:(MN-3))  007840 

80  D<MN)=T(N1)  007850 

LI  =4  tXi7860 

Y1=X  007870 

GO  TO  60  007880 

90  TRP2=Z(4)  007890 

RETURN  007900 

END  007910 

SUBROUTINE  RESTRTCCYC, A, IRSTRT)  007920 

COMMON/DATA/ II (7), R1 (3)  007930 

COMMON/RDATA  /I2(3),R2C5)  007940 

COMMON/MDATA  /I3<18)  ,R3<209)  007950 

C0MM0N/LDATA/R4(  600  ),I4(  112)  007960 

COMMON/STEPS  /I5<8)  007970 

COMMON/CORFAC  /I6,R5(3) , I7( 10) ,R6(10) , 18, R7(200) , I9,R8(220)  007980 

C0MM0N/MKCRVE/R9<10O)  007990 

COMMON/OUTPOT/RIO(3), 110(2)  008000 

REWIND  7  008010 

INOUT  =  IRSTRT  -  2  008020 

GO  TO  (100,200,200), INOUT  008030 

008040 

WRITE  RESTART  TAPE  008050 

008060 

100  WRITE(7)  I1,R1, I2,R2, I3,R3  008070 


WRITE(7)  R4,I4,  no  008080 

WRITE(7)  15, 16, R5, I7,R6, I8,R7, I9,R8,R8,R10,CYC,A  008090 

RETURN  008100 

008110 

READ  RESTART  TAPE  008120 

008130 

200  READ  (7)  I1,R1, I2,R2, I3,R3  008140 

READ  (7)  R4, 14,  no  008150 

READ  (7)  I5,I6,R5,I7,R6,I8,R7,I9,RB,R9,R10,CYC,A  008160 

WRITE(6,2000)  CYC, A  008170 

RETURN  008180 

2000  FORMATdHO,  70(1H»)/  5X , 43HRESTART  TAPE  READ.  THIS  RUN  BEGINS  AT  C008190 
lYCLE  ,  E16.8/5X,24H  WITH  A  CRACK  LENGTH  OF  ,F9.5/1X,  70(1H*))  008200 

END  008210 

SUBROUTINE  CNDIN(IEQN)  008220 

COMMON/MDATA/  MAT ID ( 18) , C, SMALLN, CARRAY ( 100) , SNARAY ( 100) , KSUBC,  008230 
/  KSUBQ,SIGMAY,DELKTH,RMULT,RCUT,OLMAX  008240 

C0MM0N/PARIS/C1,SN1,DKC0M,C2,SN2  008250 

COMMON  /DIRECT/  NDADN  008260 

COMMON  /WALKER/  CWALK,EXPM,EXPN  008270 

COMMON  /SIGMOID/DKSTAR, TOUGH, BEE, PEA, QUE, DEE 

REAL  KSUBC, KSUBQ  008280 

GO  TO  (100,300,200,400,500), lEQN  008290 

100  READ (5, 1200)  C, SMALLN, KSUBC  008300 

WRITE(6,2700)  C, SMALLN, KSUBC  008310 

GO  TO  600  008320 

200  READ(5, 1100)  PTS  008330 

NDADN  =  PTS  +  0.5  008340 

READ(5,1100)  (CARRAY(I),SNARAY(I),I=  1, NDADN)  008350 

WRITE (6, 2000)  008360 

WRITE(6,2100)(CARRAY(I),SNARAY(I),I  =  1, NDADN)  008370 

DO  250  J=l, NDADN  008380 

250  SNARAY (J)  =  ALDG10(SNARAY( J) )  008390 

GO  TO  600  008400 

300  READC5,1200)  C1,SN1,DKC0M,C2,SN2  008410 

IFfDKCOM.GT.O.)  GO  TO  350  008420 

C2  =  Cl  008430 

SN2  =  SNl  008440 

WRITE<6,2550)  Cl, SNl  008450 

GO  TO  600  008460 

350  WRITE(6,2600)  C1,SN1,DKC0M,C2,SN2,DKCQM  008470 

GO  TO  600  008480 

400  READ(5,1200)  CWALK,EXPM,EXPN  008490 

WRITE(6,2500)  CWALK,EXPM,EXPN  008500 

GO  TO  600  008510 

500  READ (5, 1200) DKSTAR, TOUGH, BEE, PEA, QUE, DEE 

WRITE (6, 2520) BEE, DKSTAR, PEA, DKSTAR, QUE, TOUGH,  DEE 


2520  FORMAT ( 1 X , 1 OHDA / DN*EXP ( , F 1 2 . 7 , 1 2H ) » ( ( DELT A-K/ , F 1 2 . 7 , 3H ) ** , F6 . 2 , 1 H ) 

1,1H*  /  1X,13H((LN<DELTA-K/,F12.7,4H))«*,F6.2,2H)* 

2/,  1X,5H((LN(,F12.7, 12H/DELTA-K)  )«,F6.2,  IH) ) 

600  READCS, 1100)  KSUBQ, SIGMAY  008530 

WRITE(6,2400)  KSUBQ,SIQMAY  008540 


FORMAT ( 2E 10.0)  008550 
FORMAT (8E10.0)  0085E0 
FORMATCIHO,  1X,32HDIRECT  INPUT  OF  DA/DN  VS  DELTA  K//34X, 7HDELTA  K  008570 


1  ,14X,5HDA/DN  )  008580 

)  FORMATC  5X,2E20.8)  008590 

)  FORMATCIHO,  1X,7HKSUBQ  =  , E15. 8, 5X, 15HYIELD  STRESS  =  ,E15.8  )  008600 

)  FORMATCIHO,  1X,3HC  =,E16.8,5X,3HM  =, F6. 4, 5X, 3HN  =,F6.3)  00S610 

)  FORMATCIHO, 1X,30HL1 NEAR  PARIS  FIT  AS  FOLLOWS  /  008620 

/  5X,E12.4, IIHCDELTA  K)**  ,F5.3)  008630 

)  FORMATCIHO,  1X,29HBILINEAR  PARIS  FIT  AS  FOLLOWS/  SX,E12.4, 11HCDELT008640 
lA  K)**,F5.3,22H  FOR  DELTA  K  LESS  THAN, E16.8/5X,E12. 4, IIHCDELTA  K)*008650 
2t  ,F5.3,22H  FOR  DELTA  K  MORE  THAN  ,E16.S  )  008660 

D  FORMATCIHO,  1X,3HC  =,E16.8,5X,aHSMALLN  =,F6.3,5X, 7HKSUBC  =,E16.8)  008670 
RETURN  008680 

END  008690 

BLOCX  DATA  008700 

COMMON/MKCRVE/  MKClOO)  008710 

REAL  MK  008720 

INITIALIZATION  VALUES  FOR  FIRST  84  ELEMENTS  OF  MK  008730 

DATA  MK/  008740 

1  11006. ,0.05,0. 10,0.20,0.30,0.40,0.50,  008750 

2  0.0,1.00,1.00,1.00,1.00,1.00,1.00,  008760 

3  0.1,1.01,1.01,1.01,1.01,1.01,1.00,  008770 

4  0.2,1.03,1.03,1.02,1.02,1.01,1.00,  008780 

5  0.3,1.06,1,06,1.04,1.03,1.02,1.00,  008790 

6  0.4,1.12,1.12,1.08,1.05,1.02,1.00,  008800 

7  0.5,1.22,1.18,1.14,1.08,1.03,1.00,  008810 

3  0.6,1.34,1.30,1.22,1.13,1.06,1.01,  008820 

9  0.7,1.48,1.42,1.31,1.20,1.08,1.02,  008830 

A  0.8,1.64,1.57,1.41,1.26,1.13,1.04,  008840 

B  0.9,1.77,1.68,1.50,1.32,1.18,1.08,  008850 

C  1.0,1.84,1.75,1.59,1.38,1.22,1.10/  008860 

END  008870 

SUBROUTINE  CELI2CRES, AK, A,B, lER)  008880 

008890 

. 0089i)0 

0089 1 0 

SUBROUTINE  CEL I 2  008920 

008930 

PURPOSE  008940 

COMPUTES  THE  GENERALIZED  COMPLETE  ELLIPTIC  INTEGRAL  OF  008950 

SECOND  KIND.  008960 

008970 

USAGE  008980 

CALL  CELI2CRES,AK,A,B, lER)  008990 

009000 

DESCRIPTION  OF  PARAMETERS  009010 

RES  -  RESULT  VALUE  009020 

AK  -  MODULUS  C INPUT)  009030 

A  -  CONSTANT  TERM  IN  NUMERATOR  009040 

B  -  FACTOR  OF  QUADRATIC  TERM  IN  NUMERATOR  009050 

lER  -  RESULTANT  ERROR  CODE  WHERE  009060 


RETURN 

END 

BLOCK  DATA 

COMMON/MKCRVE/  MK  C 100) 

REAL  MK 

INITIALIZATION  VALUES  FOR  FIRST  84  ELEMENTS  OF  MK 
DATA  MK/ 

I  11006. ,0.05,0. 10,0.20,0.30,0.40,0.50, 

Z  0.0,1.00,1.00,1.00,1.00,1.00,1.00, 

3  0. 1, 1.01,1.01, 1.01,1.01,1.01, 1.00, 

\  0.2,1.03,1.03,1.02,1. 02, 1.01,1.00, 

3  0.3,1.06,1,06,1.04,1.03,1.02,1.00, 

i  0.4,1.12,1.12,1.08,1.05,1.02,1.00, 

7  0.5,1.22,1.18,1.14,1.08,1.03,1.00, 

3  0.6, 1.34,1.30,1.22,1.13,1.06,1.01, 

3  0.7,1.48,1.42,1.31,1.20,1.08,1.02, 

^  0.8,1.64,1.57,1.41,1.26,1.13,1.04, 

B  0.9,1.77,1.68,1.50,1.32,1.18,1.08, 

;  1.0,1.84,1.75,1.59,1.38,1.22,1.10/ 

END 

SUBROUTINE  CELI2CRES, AK, A,B, lER) 


SUBROUTINE  CEL I 2 
PURPOSE 

COMPUTES  THE  GENERALIZED  COMPLETE  ELLIPTIC  INTEGRAL  OF 
SECOND  KIND. 

USAGE 

CALL  CELI2CRES,AK,A,B, lER) 

DESCRIPTION  OF  PARAMETERS 
RES  -  RESULT  VALUE 

AK  -  MODULUS  C INPUT) 

A  -  CONSTANT  TERM  IN  NUMERATOR 

B  -  FACTOR  OF  QUADRATIC  TERM  IN  NUMERATOR 
lER  -  RESULTANT  ERROR  CODE  WHERE 


o  u  o  o  o  o  o  o  a  (J  o  o  o  o  o  u  o  o  o  o  o  o  o  o 


IER=0  NO  ERROR 

IER=l  AK  NOT  IN  RANGE  -1  TO  +1 

REMARKS 

FOR  AK  =  +1,-1  THE  RESULT  VALUE  IS  SET  TO  1.E75  IF  B  IS 
POSITIVE,  TO  -1.E75  IF  B  IS  NEGATIVE. 

SPECIAL  CASES  ARE 

k(k:)  obtained  with  a  =  i,  b  =  i 

E(K)  OBTAINED  WITH  A  =  1,  B  =  CK*CK  WHERE  CK  IS 

COMPLEMENTARY  MODULUS. 

B(K)  OBTAINED  WITH  A  -^1,  B  =  0 

D(K)  OBTAINED  WITH  A  -  0,  B  =  1 

WHERE  K,  E,  B,  D  DEFINE  SPECIAL  CASES  OF  THE  GENERALIZED 
COMPLETE  ELLIPTIC  INTEGRAL  OF  SECOND  KIND  IN  THE  USUAL 
NOTATION,  AND  THE  ARGUMENT  K  OF  THESE  FUNCTIONS  MEANS 
THE  MODULUS. 

SUBROUTINES  AND  FUNCTION  SUBPROGRAMS  REQUIRED 
NONE 

METHOD 

DEFINITION 

RES= I NTEQRAL ( ( A+B*T*T ) / (SORT  <  < 1 +TCT ) * ( 1 + ( CK*T ) ) ) * ( 1 +T*T) 
SUMMED  OVER  T  FROM  0  TO  INFINITY). 

EVALUATION 

LANDENS  TRANSFORMATION  IS  USED  FOR  CALCULATION. 

REFERENCE 

R.BULIRSCH,  NUMERICAL  CALCULATION  OF  ELLIPTIC  INTEGRALS 
AND  ELLIPTIC  FUNCTIONS,  HANDBOOK  SERIES  SPECIAL  FUNCTIONS, 
NUMERISCHE  MATHEMATIK  VOL.  7,  1965,  PP.  78-90. 


IER=0 

TEST  RANGE 
CK=AK*AK 

IFi;CK-l.)  20,20,10 
10  IER=1 
RETURN 

COMPUTE  COMPLEMENTARY  MODULUS 

20  GE0=SQRT(1.0-CK) 

IF (GEO)  70,30,70 

SET  RESULT  VALUE  =  OVERFLOW 

30  IF(B)  40,60,50 
40  RES=-1.E38 
RETURN 


>009290 
009300 
009310 
009320 
009330 
009340 
009350 
009360 
009370 
. 009380 
009390 
009400 
009410 
009420 
009430 
009440 
009450 
009460 
009470 
009480 
009490 
009500 
009510 
009520 
009530 
009540 
009550 
009560 
009570 
009580 
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•" 

009120 
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009140 

009150 
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009220 
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009250 

’  -i" 
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i  I 

009280 
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50  RES=1.E38  00'35'30 

RETURN  ..'09E00 

60  RES=A  009610 

RETURN  009620 

C  009630 

C  COMPUTE  INTEGRAL  009640 

C  009650 

70  ARI=1.  009660 

AA=A  009670 

AN=A+B  009680 

W=B  009690 

80  M=W+AA»6E0  009700 

W=W+W  009710 

AA=AN  009720 

AARI=ARI  009730 

ARI=GEO+ARI  009740 

AN=W/ARI+AN  009750 

009760 

TEST  OF  ACCURACY  009770 

C  009780 

I F  ( AAR  I  -GEO- 1 .  E-4  liAAR  1)  100,100,90  009790 

90  GE0=SQRT«:QE0*AARI)  009800 

GEO=GEO+GEO  009810 

GO  TO  80  009820 

100  RES=, 785398 16*AN/AR I  009830 

RETURN  009840 

END  009850 

SUBROUTINE  GRWCRKi:CYC,A,DN)  009860 

COMflON/PDATA/  MODEL, RETARD, PLSTRN,OVLD, SIGMAX , SIGMIN, ASUBP, SMALLM  009870 
INTEGER  RETARD, PLSTRN  009880 

iIOMMON/LDATA/SMAXf  20, 10) , SMIN(  20, 10) , CYCLES!  20, 10) , NLYRS( 10) ,  009890 

/  NBLKS,IBLKS(  50  ),ISE6S<  50  ),NSEGS  009900 

COMMON/STEPS/ ISEG,J1,J2,J3,J4,J5, ISTOP,NORTRD  009910 

EXTERNAL  RATE, WHELER, WLNBRG  009920 

SIGMAX  =  SMAXf J4, ISEG)  009930 

SIGMIN  =  SMIN(J4, ISEG)  009940 

IF f MODEL. 6T.0)  GO  TO  100  009950 

CALL  RK1DES(CYC,A,DN,RATE)  009960 

RETURN  009970 

100  GO  TO  (110, 120, 130, 140, 150), MODEL  009980 

110  CALL  YLD ZNE( CYC, A, DN, WHELER)  009'390 

GO  TO  200  010000 

120  CALL  YLD ZNE( CYC, A, DN, WLNBRG)  010010 

GO  TO  200  010020 

130  CONTINUE  010030 

CALL  CLOSUR(CYC,A,DN)  010040 

GO  TO  200  010050 

140  CONTINUE  010060 

150  CONTINUE  010070 

200  RETURN  010080 

END  010090 


SUBROUTINE  TRANS(A, ATRANS, CYCTR) 

COMMON/LDATA/SMAX C  20,10),SMIN(  20, 10) , CYCLESC  20, 10) , NLYRSC 10) , 

/  NBLKS, IBLKS(  50  ),r3EGS(  50  ) , NSEGS 

COMMON/STEPS/ ISEG, J1 , J2, J3,  J4, J5, ISTOP, NDRTRD 
COMMON/CORFAC/  I SURF , RAT I 0 , PH I , TH I CK , I BETA ( 1 0 ) , BETA ( 1 0 ) , NPTS , 

/  ADVERB ( 1 00 ) , BT ABLE ( 1 00 ) , NPTS2 , A0VRB2 ( 1 00 ) , 

/  BTABL2 ( 1 00 ) , AST ART ( 1 0 ) , ASTOP ( 1 0 ) 

ATRANS  =  THICK 
I F ( A . LT . ATR ANS ) RETURN 
100  AEFF  =  ATRANS/(2.;l:RATI0) 

A  =  AEFF 
ISTOP  =  2 

WRITE (6, 1000)  A,CYCTR 
1000  FORMATCIHO,  70(1H*)/  5X, 55HTRANSITIQN  TO  A  THRU  CRACK  OF  EFFECT I VEO 10230 


1  LENGTH,  AEFF  =  ,FG.5,4H  AT  ,F12.2,7H  CYCLES/IX,  70(1H*))  010240 

RETURN  010250 

END  010260 

SUBROUTINE  RATE  (CYCLE, A, DADN)  010270 

COMMON/DATA/  EQN, NASA, J IPR, J2PR, J3PR, J4PR, J5PR, AZERO, AMAX , NZERO  010280 
INTEGER  EQN  010230 

REAL  NZERO  010300 

COMMON/MDAT A/  MAT I D ( 1 8 ) , C , SMALLN, CARRAY (100), SNARAY (100), KSUBC ,  010310 

J  KSUBQ,SIGMAY,DELKTH,RMULT,RCUT,OLMAX  010320 

REAL  KSUBC,  KSUBQ  010330 

COMMON/LDATA/SMAX (  20,10),SMIN(  20, 10) , CYCLESC  20, 10) , NLYRS( 10) ,  010340 

/  NBLKS, IBLKSC  50  ),ISE6S(  50  ), NSEGS  010350 

CQMMON/STEPS/ISEG, J1,J2, J3, J4, J5, ISTOP,NORTRD  010360 

COMMON/CORFAC/  ISURF, RATIO, PHI, THICK, IBETAC 10) , BETAC 10) , NPTS,  010370 

/  AOVERB ( 1 00 ), STABLE ( 100 ),NPTS2,A0VRB2( 100),  010380 

/  BTABL2(100),ASTART(10),ASTDP(10)  010330 

COMMON/OUTPOT/  KMAX,KMAXA,DELTAK, IFLT,DADNPR  010400 

COMMON/PAR I S/C 1 , SN 1 , DKCOM, C2, SN2  0 1 04 1 0 

COMMON/DIRECT/  NDADN  010420 

COMMON  /WALKER/  CWALK, EXPM,  EXPN  010430 

COMMON  /S IGMO ID/DKSTAR, TOUGH, BEE, PEA, QUE, DEE 

REAL  KMAX  010440 

REAL  KMAX A  010450 

CALL  DELTA(A,DELTAK,KMAX,R)  010460 

IF (ISTOP  .NE.  0)Ga  TO  575  010470 

IF(  R.  GE.  RCUT)  R  =  RCUT  010480 

CALL  K(SMAX(J4, ISEQ),A,KMAXA)  010430 

IF(DELTAK.LE.O.O)  GO  TO  300  010500 

IF ( ISURF. EQ.O)  GO  TO  50  010510 

CALL  TRANS (A, ATRANS, CYCLE)  010520 

IF (A. LT. ATRANS)  GO  TO  50  010530 

RETURN  010540 

50  THRSLD  =  DELKTH:l:(  1 . 0-RMULT*R)  010550 

IF  (DELTAK.LE. THRSLD)  GO  TO  300  010560 

GO  TO  (100, 200, 230, 240, 250), EQN  010570 

100  DENOM  =  (1.0  -  R)  X  KSUBC  -  DELTAK  010580 

IF(KMAXA.GE.KSUBQ)  GO  TO  400  010530 


010100 

010110 

010120 

010130 

010140 

010150 

010160 

010170 

010180 

010130 

010200 

010210 

010220 


IF(DENOM.LE.O.O)  GO  TO  525 


010600 


DADN= ( C* ( DELTAK ) **SMALLN ) / DENOM  010610 

RETURN  010620 

200  DELKC  =  KSUBQ  -  KMAXA  010630 

IF(DELKC.LE.O.O)  GO  TO  400  010640 

C  =  Cl  010650 

SMALLN  =  SNl  010660 

IF (DELTAK. GE.DKCOM)  C  =  C2  010670 

IF ( DELTAK. GE.DKCOM)  SMALLN  =  SN2  010680 

DADN=C*( DELTAK )**SMALLN  010690 

RETURN  010700 

230  DADN  =  TBLKUP(CARRAY,SNARAV,NDADN, 100, DELTAK)  010710 

DADN  =  lO.iWDADN  010720 

IF (KMAXA. GE. KSUBQ)  GO  TO  400  010730 

RETURN  010740 

240  IF (KMAXA. QE. KSUBQ)  GO  TO  400  010750 

DADN=CWALK* (DELTAK/ ( ( 1 . -R)** ( 1 . -EXPM) ) );C*EXPN  010760 

RETURN  010770 

250  IF  (KMAXA  .QE.  KSUBQ)  GO  TO  400 
T0UG=T0U6H»(l.-.4) 

IF  (DELTAK  .QE.  TOUQ)  GO  TO  260 
XKSTAR=DKSTAR» ( 1 . - . 4 ) 

IF  (DELTAK  .LE.  XKSTAR)  GO  TO  300 

DADN=E  XP (BEE ) * ( ( DELTAK/ XKSTAR ) **PEA) * ( ( ALOG ( DELTAK/ XKSTAR ) ) t*DUE ) 
l:t(  (ALOG(TOUQ/DELTAK)  )»*DEE) 

RETURN 
260  IST0P=1 

WRITE(6,270) 

270  FORMATCIHO,  70(1HC)/  5X, 50HDELTA-K  EXCEEDS  THE  TOUGHNESS  PROB.  IS 
1  TERMINATED/1X,70(1H»)/1H0,  IX,26HLAST  CALCULATED  VALUES  ARE/// 

2) 

GO  TO  575 

300  DADN=0.0  010780 

RETURN  010790 

400  I STOP  =  1  010800 

WRITE(6,500)  010810 

500  FORMATdHO,  70(1H*)/  5X,4&HKMAX  APPLIED  EXCEEDS  KSUBQ.  PROBLEM  TER010820 

IMINATED/IX,  70(1H*)/1H0,  1X,26HLAST  CALCULATED  VALUES  ARE///)  010830 

GO  TO  575  010840 

525  WRITE(6,550)  010350 

550  FORMATdHO,  70(1H»)/  5X,49HDELTA  K  EXCEEDS  (l-R)KSUBC.  PROBLEM  IS  010860 

1TERMINATED/1X,70(1H;|:)/1H0,  1X,26HLAST  CALCULATED  VALUES  ARE///)  010870 

I STOP  =  1  010880 

575  CONTINUE  010890 

URITE(6,600)  J1,J2,ISE6, IFLT,J4, CYCLE, A,KMAXA,KMAX, DELTAK, DADN  010900 

600  FORMAT(  5X, 18HBL0CK  IN  SPECTRUM  ,14/  010910 

/  5X, 18HSEGMENT  NUMBER  ,14/  010920 

/  5X,1QHMISSI0N  NUMBER  ,14/  010930 

/  5X,18HFLIGHT  NUMBER  ,16/  010940 

/  5X, 18HLAYER  IN  MISSION  ,14/  010950 

5  5X, 18HACCUMULATED  CYCLES,E16.8/  010960 

6  5X, ISHCRACK  LENGTH  ,E16.8/  010970 

7  5X,18HKMAX  APPLIED  ,E16.8/  010980 
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a  5X,18HKMAX  EFFECTIVE  ,E1£.8/ 

9  5X,18HDELTA  K  ,816.3/ 

/  5X,18HDA/DN  ,£16. 8) 

RETURN 

END 

SUBROUTINE  KISIGMA, A,KM) 

i;OMMON/CORFAC/DUMMY(  14) ,  BETA ( 10)  .DUMMY  1  (422) 

COMMON/DATA/  EQN, NASA, JIPR, J2PR, J2PR, J4PR, J5PR, AZERO, AMAX, NZERO 
INTEGER  EQN 
REAL  NZERO 
REAL  KM 

DATA  PI/3.14159265/ 

CALL  BETAS  (SIGMA, A, BET AT, Q) 

IF(NASA.NE.O)  GO  TO  100 
KM=S  I  GMA*SQRT  (  P I  »A )  1:BET  AT 
RETURN 

100  KM=SIGMA4:SQRT(A)1:BETAT 
RETURN 
END 

SUBROUTINE  BETAS(SIGMA,A,BETAT,Q) 

COMMON/MDAT A/  MATID(18),C, SMALLN , CARRAY (100), SNARAY (100), KSUBC , 
/  KSUBQ , S I GMAY , DELKTH, RMULT , RCUT , OLMAX 

REAL  KSUBC, KSUBQ 

COMMON/STEPS/ ISEG, J1 , J2, J3, J4, J5, ISTOP, NORTRD 
COMMON/CORFAC/  ISURF, RATIO, PHI, THICK, IBETA(IO) , BETA( 10) , NPTS, 

/  AQVERB( 100) , BTABLE( 100) , NPTS2, A0VRB2 ( 100) , 

/  BTABL2 ( 1 00 ) , ASTART ( 1 0) , ASTOP ( 1 0 ) 

CQMMON/MKCRVE/  MK(IOO) 

REAL  MK,MSUBK,M1 
BETAT=1.0 
MSUBK=1 . 0 
0=1.0 

PI  =  3.14159265 
5  DO  100  1=1,10 
J=IBETA(I) 

IF(J.EQ.O)  GO  TO  100 

IF(A.LT.ASTART(J).  OR  .A.GT.ASTOP(J))  GO  TO  100 
GO  TO ( 10,20,30,40, 50,60,70,80,90) , J 

CONSTANT  MULTIPLIER 

10  BETAT  =  BETAT  *  BETA(J) 

GO  TO  100 

FINITE  WIDTH  SECANT  CORRECTION 

20  HOLE  =  BETA(5)+BETA(6) 

SMALLC=A 

IF ( ISURF. NE.O)  SMALLC=A/(2.*RATI0) 

SMALLK  =  (SMALLC  +  HOLE) /BETA( J) 

IF(SMALLK.GT.0.5)  GO  TO  200 

BETAT  =  BETAT*SQRT(1./ COS (PI ♦SMALLK)) 
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GO  TO  100 


TABULAR  CORRECTION  FACTOR 
30  SMALLK=A/BETA(J) 

BET AT=BET ATCTBLKUP ( AOVERB , BT ABLE , NPTS ,100, SMALLK ) 

GO  TO  100 

SECOND  TABULAR  CORRECTION  FACTOR 
40  SMALLK=A/BETA(J) 

BETAT=BETAT*TBLKUP(A0VRB2, BTABL2,NPTS2, 100, SMALLK) 

GO  TO  100 

BOWIE  SOLUTION  FOR  SINGLE  CRACK  FROM  CIRCULAR  HOLE 

50  BETAT  =  BETAT  *  (0.6762062  +  (0.8733015/ (0.3245442+A/BETA( J) )) ) 
GO  TO  100 

BOWIE  SOLUTION  FOR  DOUBLE  CRACK  FROM  CIRCULAR  HOLE 

60  BETAT  =  BETAT  t  (0.3438510  +  (0. 6805078/(0. 277 1 965+A/BETA( J) )) ) 
GO  TO  100 

SOLUTION  FOR  ASTM  COMPACT  TENSION  SPECIMEN  USING  J.  C.  NEWMAN 
EQUATION  (12)  FROM  'STRESS  ANALYSIS  OF  THE  CIDMPACT  SPECIMEN 
INCLUDING  THE  EFFECTS  OF  PIN  LOADING' 

70  AW  =  A/BETAd) 

BETAT  =  BETAT:»:(4. 55-40. 32:l:AW+414.7«AW.»:1:2. -1698.  tAW.t1:3. 

1  -•■3781 . 1:AW**4. -4287.  YAW KtS. -^2017.  YAWYY6. )  /SORT (BETA ( I )  YAYPI )  / 

2  THICK 
GO  TO  100 

SOLUTION  FOR  tSRUMMAN  COMPACT  TENSION  SPECIMEN 
H/W  =  .35  D/W  =  .25 

30  AW  =  A/BETA(I) 

POLY  =  .  1 229-^  1 6 . 4098YAW-37 . 395YAWYY2 .  -^54 . 7667 YAWY Y3 . 

IF  i:  AW  .LE.  0.5  )  GO  TO  81 

POLY  =  114.054-830. 132YAW-^2327.177YAWYY2.-2890.B11YAWYY3. 

1  •■13a2.306YAWYY4. 

81  BETAT  =  BET ATYPOLY/SQRTCPIYA) /THICK 

GO  TO  100 

SOLUTION  FOR  COMPACT  TENSION  SPECIMEN  FROM  ASTM  E  747-83 
CONSTANT-LOAD-AMPLITUDE  FATIGUE  CiRACK  GROWTH  RATES  ABOVE 
10-08  M/CYCLE 
30  AW=A/BETA(I) 

BET AT=BET ATY  ( 2 .  AW )  Y  ( .  886-^4 . 64YAW- 1 3 . 32Y AWYY2  .+14.  72YAWY Y3 . 

1  -5 . 6 YAWY Y4 . ) / ( ( TH I CK Y ( 1 - AW ) YY 1 . 5 ) YSQRT ( BET A ( 1) YP I YA ) ) 

100  iCONTINUE 

IF(ISURF.EQ.O)  RETURN 

SURFACE  FLAW  CORRECTION 
FROM  NEWMAN  -  NASA  TN  D-8244 
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011510 
011520 
011530 
011540 
011550 
011560 
011570 
011580 
011590 
011600 
011610 
011620 
011630 
011640 
011650 
011660 
011670 
011680 
011630 
011700 
011710 
011720 
011730 
011740 
011750 
011760 
011770 
011730 
011730 
011800 
011810 
01 1820 
011830 
011840 
011850 
011S60 
011870 
011880 
011330 
011900 
01 1310 


011330 

011340 

011350 

011360 

011370 


.'v 


'  - '  • 

>  -‘-V- 

I . . 


•V.-, 


c 

Q=PH  I U-Z .  0-0 . 2 1 2*  C S I QMA/S I GMAY )  **2 . 0 
A0VERT=A/THICK 
AaVERC=2.:t:RATI0 
P=2.+8.*A0VERC**3. 

Ml  =  1.13-0.1:|:AaVERC 
MSUBK=  (SORT  (Q/AOVERC)  -Ml )  .1:A0VERT**P 
BETAT=BETAT*(M1+MSUBK) /SQRT(Q) 

RETURN 

200  WRITE (6, 1000) 

1000  FORMATCIHO,  TOCIH*)/  5X,52HCRACK  LENGTH  EXCEEDS  PLATE  WIDTH.  TERMI012080 


INATE  PROBLEM.  /  TOCIH*))  012090 

ISTOP  =  1  012100 

RETURN  012110 

END  012120 

SUBROUTINE  YLDZNECCYC, A, DN, FR)  012130 

C  012140 

C  THIS  ROUTINE  CONTROLS  APPLICATION  OF  THE  WHEELER  MODEL  AND  THE  012150 

C  EFFECTIVE  STRESS (WILLENBORG) MODEL  BASED  ON  THE  PROGRESS  THRU  A  012160 

C  YIELD  ZONE (ASUBP) DUE  TO  AN  OVERLOAD.  012170 

C  012180 

COMMON/RDATA/  MODEL , RETARD , PLSTRN, OVLD , S I GMAX , S I GM I N , ASUBP , SMALLM  012190 
INTEGER  RETARD, PLSTRN  012200 

COMMON/LDATA/SMAX(  20,10),SMIN(  20, 10) , CYCLES!  20, 10) , NLYRSC 10) ,  012210 

/  NBLKS,IBLKS(  50  ),ISEGS(  50  ),NSEGS  012220 

COMMON/MDATA/  MAT I D  < 1 8 ) , C , SMALLN, CARRAY (100), SNARAY (100), KSUBC ,  012230 

/  KSUBQ,SIGMAY,DELKTH,RMULT,RCUT,OLMAX  012240 

REAL  KSUBC, KSUBQ  012250 

COMMON/STEPS/ I SEG, J1 , J2, J3, J4, J5, ISTOP, NORTRD  012260 

REAL  KMAXjKOVLD  012270 

EXTERNAL  FR,RATE,DNDA  012280 

CYCF  =  CYC  +  DN  012290 

IF (RETARD. EQ. 1)  GO  TO  100  012300 

C  012310 

C  IS  THIS  LAYER  SUBJECT  TO  RETARDATION  S.  012320 

C  012330 

IF(SMAX(J4,ISEG).GE.0VLD)  GO  TO  800  012340 

C  012350 

C  RETARDATION  IS  APPLIED.  012360 

C  012370 

RETARD  =  1  012380 

C  012390 

C  DETERMINE  EXTENT  OF  ELASTIC-PLASTIC  INTERFACE  012400 

C  012410 

CALL  K(OVLD,A,KOVLD)  012420 

RSUBY  =  RY(KOVLD, PLSTRN)  012430 

ASUBP  =  A  +  RSUBY  012440 

C  012450 

C  WILL  FIRST  CYCLE  OF  THIS  LAYER  CAUSE  ASUBP  TO  BE  EXCEEDED  &  012460 

012470 
012480 
012490 


011980 

011990 

012000 

012010 

012020 

012030 

012040 

012060 

012070 


100  CALL  SIGEFF(A,ASUBP,SIGMAY) 
IF (RETARD. EQ.O)  GO  TO  800 


o  o  o  o  o  o  o  o  o  o  o  o  o  o  o 


CALL  K(SIGMAX,A,KMAX)  012500 

RSUBYl  =  RY(KMAX,PLSTRN)  012510 

IF(A+RSUBY1.LT.ASUBP)  GO  TO  200  012520 

RETARD  =  0  012530 

GO  TO  800  012540 

012550 

THIS  LAYER  IS  SUBJECT  TO  RETARDATION.  012560 

ASSUME  THAT  RETARDATION  APPLIES  OVER  THE  ENTIRE  LAYER.  012570 

012580 

200  CYCl  =  CYC  012590 

AA  =  A  012600 

CALL  RK1DES(CYC1,AA,DN,FR)  012610 

IF(ISTOP.NE.O)  RETURN  012620 

IFCAA.GT.ASUBP)  GO  TO  400  012630- 

CALL  K<SIGMAX,AA,KMAX)  012640 

RSUBY  =  RY(KMAX,PLSTRN)  012650 

012660 

CHECK  ASSUMPTION.  012670 

012680 

IF<AA+RSUBY.GE.ASUBP)  GO  TO  400  012690 

CYC  =  CYCl  012700 

A  =  AA  012710 

RETURN  012720 

012730 

ENTIRE  LAYER  IS  NOT  RETARDED  012740 

CALCULATE  DELTA  A  AND  ITS  ASSOCIATED  YIELD  ZONE(RSUBY)  SUCH  THAT  012750 
(A+DA+RSUBY  =  ASUBP)  TO  A  GIVEN  TOLERANCE (TOL)  012760 

012770 

400  TOL  =  l.OE-5  'M2780 

RED  =  1.0  012790 

SUBT  =  0.1  012800 

CYCl  =  CYC  012810 

AA  =  A  012820 

500  DA  =  (ASUBP-CAA+RSUBYD)  t  RED  012830 

CALL  SIGEFF<AA+DA,ASUBP,SIGMAY)  012840 

CALL  K(SIGMAX,AA+DA,KMAX)  012850 

RSUBY  =  RYCKMAX.PLSTRN)  012860 

IF CAA+DA+RSUBY.LT. ASUBP)  GO  TO  700  012870 

RED  =  RED  -  SUBT  012880 

GO  TO  500  012890 

600  RED  =  RED  +  SUBT  012900 

SUBT  =  SUBT/10.  012910 

GO  TO  500  012920 

700  IFCSUBT.GT.TOL)  GO  TO  600  012930 

012940 

KNOWING  DN,  DETERMINE  NUMBER  OF  CYCLES  IN  THIS  LAYER  REQUIRED  012950 

TO  PRODUCE  DA.  012960 

C  012970 

CALL  RUNKUT(AA,CYC1,DA,DNDA)  012980 

IF(ISTOP.NE.O)  RETURN  012990 

CYC  =  CYCl  013000 

A  =  AA  013010 


C 

c 

c 

c 


REMAINDER  OF  LAYER  IS  NOT  SUBJECT  TO  RETARDATION. 
USE  UNRETARDED  DADN  FUNCTION (RATE) 

RETARD  =  0 
DN  =  CYCF  -  CYC 
IF(DN.LE.O.O)  RETURN 


C 

C 


c 

c 


c 

c 

c 

c 


UNRETARDED  DADN  FUNCTION 

800  SISMAX  =  SMAX(J4,ISEG) 

SIGMIN  =  SMIN(J4, ISE6) 

OVLD  =  SI6MAX 

CALL  RK1DES(CYC,A,DN,RATE) 

RETURN 

END 

REAL  FUNCTION  RY<K,PLSTRN) 

COMMON/MDATA/  MATID ( 18) , C, SMALLN, CARRAYI 100) , SNARAY ( 100) ,KSUBC, 

/  KSUBQ, SIGMAY, DELKTH, RMULT, RCUT, OLMAX 

REAL  KSUBC, KSUBQ 
REAL  K 

INTEGER  PLSTRN 

DATA  PI,RQQT2  /3. 1415926, 2.828428/ 

IF(K.LE.O.)  GO  TO  999 

RY  =  <<K/SIGMAY)*»2.)/(2.»P1) 

IF(PLSTRN.NE.O)  RY  =  RY/R00T2 
RETURN 
999  RY=  0.0 
RETURN 
END 

SUBROUTINE  DNDA(A,CYC,DDNINV) 

THIS  ROUTINE  CALCULATES  THE  NUMBER  OF  CYCLES  IN  A  LAYER  OVER 
WHICH  RETARDATION  IS  APPLIED 

COMMON/RDATA/  MODEL , RETARD , PLSTRN, OVLD , S I GMAX , S I GM I N , ASUBP , SMALLM 
INTEGER  RETARD, PLSTRN 
GO  TO  (10, 20), MODEL 
10  CALL  WHELER<CYC,A,DADN) 

GO  TO  30 

20  CALL  WLNBRG(CYC,A,DADN) 

30  DDNINV  =  1.0/DADN 
RETURN 
END 

SUBROUTINE  ULNBRGCCYC, A,DADN) 

THIS  ROUTINE  BRINGS  THE  EFFECTIVE  STRESSES  GENERATED  BY  THE 
WILLENBORG  MODEL  INTO  THE  UNRETARDED  GROWTH  RATE  EQUATIONS. 

COMMON/RDATA/  MODEL, RETARD, PLSTRN, OVLD, SI GMAX, SIGMIN, ASUBP, SMALLM 
INTEGER  RETARD, PLSTRN 

COMMON/MDATA/  MATIDC 18) ,C,SMALLN,CARRAY( 100) ,SNARAY( 100) , KSUBC, 


013020 

013030 

013040 

013050 

013060 

013070 

013080 

013090 

013100 

013110 

013120 

013130 

013140 

013150 

013160 

013170 

013180 

013190 

013200 

013210 

013220 

013230 

013240 

013250 

013260 

013270 

013280 

013290 

013300 

013310 

013320 

013330 

013340 

013350 

013360 

013370 

013380 

013390 

013400 

013410 

013420 

013430 

013440 

013450 

013460 

013470 

013480 

013490 

013500 

013510 

013520 

013530 
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c 

c 

c 


c 

c 


/  KSUBQ, S IGMAY, DELKTH,RMULT, RCUT, OLMAX 

013540 

*  .  'v'v 

REAL  KSUBC.KSUBQ 

013550 

CALL  SIGEFF(A,ASUBP,3IGMAY) 

013560 

N.  . 

CALL  RATE(CYC,A,DADN) 

013570 

RETURN 

013580 

END 

013590 

SUBROUTINE  SIGEFF(A, ASUBP.SIGMAY) 

013600 

013610 

THIS  ROUTINE  COMPUTES  THE  EFFECTIVE  STRESSES  FOR  USE  IN  THE 

013620 

WILLENBORG  MODEL.  THE  EFFECTIVE  STRESSES  ARE  STORED  IN  LOCATIONS 

013630 

’SIGMAX*  AND  ’SIGHIN’  IN  COMMON  BLOCK  /RDATA/. 

013640 

013650 

013660 

COMMON/DATA/ I DUMl , NASA, IDUM2 (5) , DUM3(3) 

013670 

COMMON/RDATA/  MODEL, RETARD, PLSTRN, OVLD, SIGMAX , SIGMIN, DUMMY, SMALLM 

013680 

INTEGER  RETARD, PLSTRN 

013690 

COMMON/LDATA/SMAX (  20,10),SMIN(  20, 10) , CYCLES!  20, 10) ,NLYRS(10) , 

013700 

/  NBLKS,IBLKS<  SO  ),ISEGS(  50  ),NSEGS 

013710 

COMMON/STEPS/ ISEG, J1 , J2, J3, J4, J5, ISTOP, NORTRD 

013720 

'  ■  V 

COMMON/MDATA/MATIDC 18) , C, SMALLN, CARRAY! 100) , SNARAY! 100) , KSUBC, 

013730 

1  KSUBQ, DUMMYY, DELKTH, RMULT, RCUT, OLMAX 

013740 

*  • ..  • .  • 

REAL  KSUBC, KSUBQ, KMAX 

013750 

013760 

PUT  APPLIED  STRESSES  IN  SIGMAX  AND  SIGMIN 

013770 

013780 

SIGMAX  =  SMAX(J4,ISEG) 

013790 

SIGMIN  =  SMIN(J4, ISEG) 

013800 

v- 

IF(M0DEL.EQ.2)  GO  TO  100 

013810 

1^-  . 

RETURN 

013820 

;• 

100  CALL  BETAS<SIGMAX,A,BETAT,QMAX) 

013830 

IF(A.GT.ASUBP)  GO  TO  200 

013840 

'  * 

SIGREF  =  (SIQMAY*SQRT(2.0*<ASUBP-A)/A))/BETAT 

013850 

IF(NASA.NE.O)  SIGREF=SIGREF.1:SQRT(3. 1415926) 

013860 

*  »■ 

IF (PLSTRN. NE.O)  SIGREF  =  SIGREF  *  SORT <2. 828428) 

013870 

1.  ,  . 

SIGRED  =  SIGREF  -  SIGMAX 

013880 

•  *  -/'■/ 

CALL  DELTA (A, DELTAK, KMAX, R) 

013890 

THRSLD  =  DELKTH* (l.-RMULTTR) 

013900 

PHI  =  (;i.-THRSLD/KMAX)/(OLMAX-l.) 

013910 

.-''1 

IF (OLMAX. EQ.O.)  PHI  =  1. 

013920 

SIGRED  =  PHI  *  SIGRED 

013930 

I 

IF (SIGRED. LE. 0.0)  GO  TO  200 

013940 

SIGMAX  =  SIGMAX  -  SIGRED 

013950 

IF (SIGMAX.lt. 0.0)  SIGMAX  =  0.0 

013960 

■  ■- 

SIGMIN  =  SIGMIN  -  SIGRED 

013970 

IF(SIGMIN.LT.O.O)  SIGMIN  =  0.0 

013980 

•  - 

RETURN 

013990 

200  RETARD  =  0 

014000 

RETURN 

014010 

•.**.**• 

END 

014020 

SUBROUTINE  WHELER(CYC, A,DADN) 

014030 

014040 

THIS  ROUTINE  APPLIES  THE  WHEELER  CORRECTION  TO  THE  UNRETARDED 

014050 

1 _ 

1.  • 


116  - 


O  O  O  K~>  o 


C  GROWTH  RATE  014060 

C  014070 

COMMON/RDATA/  MODEL, RETARD, PLSTRN, OVLD, SI6MAX, SIGMIN, ASUBP, 3MALLM  014080 


INTEGER  RETARD, PLSTRN 

0140G0 

CSUBP  =1.0 

014100 

c 

014110 

c 

DETERMINE  EXTENT  OF  CURRENT  YIELD  ZONE 

014120 

c 

014130 

CALL  K(SIGMAX,A,KMAX) 

014140 

RSUBY  =  RYCKMAX, PLSTRN) 

014150 

IF(A  +  RSUBY  .GE.  ASUBP)  GO  TO  20 

014160 

c 

014170 

c 

CALCULATE  WHEELER'S  RETARDATION  PARAMETER 

014180 

c 

014190 

CSUBP  =  CRSUBY/CASUBP-A))**  SMALLM 

014200 

20 

CALL  RATE(CYC,A,DADN) 

014210 

DADN  =  CSUBP  t  DADN 

014220 

RETURN 

014230 

END 

014240 

SUBROUTINE  RKIDESICYC, A,DCYC,F) 

014250 

c 

014260 

EXTERNAL  F 

014270 

COMMON/STEPS/ I SEQ, J1 , J2, J3, J4, J5, ISTOP, NORTRD 

014280 

IFIDCYC  .GE.  1.)  GO  TO  300 

CYCF=CYC+DCYC 

014300 

A0=A 

014310 

50 

H=0.005*A0 

014320 

A1=A0+H 

014330 

AGR0W=A0+2.0*H 

014340 

CALL  F ( CYC, Al, DADN) 

014350 

IFCDADN.LE.O.)  GO  TO  200 

014360 

IF C ISTOP. NE.O)  GO  TO  250 

014370 

DCYCR=2.TH/DADN 

014380 

IF(DCYCR.GT.DCYC)  GO  TO  100 

014390 

CYC=CYC+DCYCR 

014400 

DCYC=DCYC-OCYCR 

014410 

AO=AGROW 

014420 

GO  TO  50 

014430 

100 

DA=(CYCF-CYC)*DADN 

014440 

A0=A0+DA 

014450 

200 

CYC=CYCF 

014460 

A=A0 

014470 

250 

RETURN 

014480 

300 

CALL  RUNKUT(CYC,A,DCYC,F) 

014490 

RETURN 

014500 

END 

014510 

FUNCTION  TBLKUP<X,Y,N,NMAX,ARG) 

014520 

c 

014530 

c 

014540 

DIMENSION  X(NMAX),Y(NMAX) 

014550 

DO  10  1  =  1, N 

014560 

IF(X(I)-ARG)  10,20,20 

014570 

o  o  o  o 


p.  n.  n .■  I  ■  u  •  I.’  •  ■ 


10  CONTINUE 
I=N 

20  IF(I-1)30,30,40 
30  1*2 

40  SLaPE=(Y(r)-Y(r-i))/(X(i)-x(r-i)) 
TBLKUP=SLOPEt.  ( ARQ-X  <  I  - 1) )  +Y  ( I  - 1 ) 
RETURN 
END 

SUBROUTINE  CLOSURICYC, A1 , ON) 


COMMON/OUTPOT/XK, XKA, XKEFFN, IFLT,DADNPR 
COMMON/CLOS/CF , CCOEF , CFE  XP , B , BOL , NSAT 
COMMON/CLOSIC/SC, SPEAK, PRVMX , APEAK, SCI , SC2, SC3, PRVMN 
REAL  ISUM,KTHS6,KCOEF,KEXP,NOL,NSAT,NPR,NPREV,KTH 

DIMENSION  Q(6),QQ(2) 

DATA  Pr/3. 14159265/ 


CLOSURE  FUNCTION 

CL0SE(61,a2)  =  G1*(CC0EF  +  (CF  -  CCOEDUCl.  +  G2)  t^JCFEXP) 
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014500 

014590 

014600 

014610 

014620 

014630 

014640 

014650 

014660 

014670 


c 

EQUIVALENCE 

BETWEEN  CRACKS  2  AND  CLOSUR 

014680 

c 

CRACKS  2 

CLOSUR 

014690 

c 

ASUBP 

AP 

014700 

c 

SMAX 

S 

014710 

c 

CYCLES 

CYCLS 

014720 

p 

ISEG 

I 

014730 

c 

J4 

K 

014740 

c 

SMALLN 

AN  (NOT  USED  AS  SMALLN  ANYWAY) 

014750 

c 

KSUBQ 

AKC 

014760 

c 

DELKTH 

KTH 

014770 

c 

A 

A1 

014780 

c 

AMAX 

AF 

014790 

c 

014800 

c 

014810 

COMMON/DATA/  EQN,NASA, JIPR, J2PR, J3PR, J4PR, J5PR,AZER0,AF  ,NZERO 

014820 

INTEGER  EQN 

014830 

REAL  NZERO 

014840 

COMMON/RDATA/  MODEL, RETARD, PLSTRN, OVLD, SIGMAX, SIGMIN, AP, SMALLM 

014850 

INTEGER  RETARD, 

PLSTRN 

014860 

COMMON/MDATA/  MAT I D (1 8) ,  C, AN, CARRAY (100), SNARAY (100), KSUBC, 

014870 

/ 

AKC,SIGMAY,  KTH,RMULT,RCUT,OLMAX 

014880 

REAL  KSUBC 

014890 

COMMON/LDATA/ 

S(  20,10),SMIN(  20, 10),CYCLS(  20, 10) , NLYRS( 10) , 

014900 

/  NBLKS,IBLKS(  50  ),ISEGS<  50  ),NSEGS 

014910 

COMMON/STEPS/ 

K , J 1 , J2 , J3 , I , J5 , ISTOP , NORTRD 

014920 

COMMON/CORF AC/ 

ISURF, RATIO, PHI, THICK, IBETA( 10), BETA( 10), NPTS, 

014930 

/ 

AOVERB ( 1 00 ) , STABLE ( 1 00 ) , NPTS2 , A0VRB2 (100), 

014940 

/ 

BTABL2( 100) , ASTART ( 10) , ASTOP ( 10) 

014950 

COMMON/PAR IS/Cl 

,SNl,DKCaM,C2,SN2 

014960 

COMMON  /DIRECT/ 

NDADN 

014970 

014980 

014990 

015000 

015010 

015020 

015030 

015040 

015050 

015060 

015070 

015080 


o  o  o  o  o  o  o  r>  c.i  o  r>  n  o  o  o 


L*  -v'  "i.  •wl 


FUNCTION  FOR  BC  LT  OR  EQ  SC  INITIAL 
D0WN<Z1,Z2,Z3)  =  Z1  -  (Z1-SC3)*(Z2/Z3) 


FUNCTIONS  FOR  SC  GT  3C  INITIAL 
SCONE <SC3)  =  SC3CB0L 


FUNCTION  FOR  INCREASING  CLOSURE  STRESS 
NPREV(SC,SC3,SC11)  =1. +(NSAT-1.  ):|:(SC-SCl  1  )/<SC3-SCl  1) 


INITIALIZE  PARAMETERS 
I GROW  =  0 
MODE  =  0 
ITEM  =  1 
I GROM  =  1 

CYSUM  =  CYCLSd.K) 

NOL  =  0. 

ASTRT  =  A1 
KLU  =  1 
R  =  0. 

IF  (  S(I,K)  .NE.  0.  )  R  =  SMIN(I,K)/S(I,K) 
SMNGR  =  SMINd.K) 

IF  (  SPEAK  .NE.  0.  )  GO  TO  30 
SINITL  =  Sd,K) 

IF  (  SINITL  .LE.  0.)  SINITL  =  0.05»SIGMAY 
R  =  SMINd,K)/SINlTL 
CALL  BETAS < S (I , K ) , A 1 , ALP , QE ) 

XK  =  SINITL*SaRT(PI*Ani:ALP 
G1  =  SINITL 
G2  =  R 

IF  (  Q2  .LT.  -1.)  G2  =-l. 

SC  =  CLOSE  i:G1,Q2) 

SCI  =  SC 
SC2  =  SC 
SC3  =  SC 
SPEAK  =  SINITL 
PRVMX  =  SINITL 
PRVMN  =  SMINdpK) 

APEAK  =  A1 

AP  =  A1  +  RY(XK,PLSTRN) 

0MGA2  =  AP 
30  CONTINUE 


START  ANALYSIS 
100  ISUM  =0. 


Q1  =  S(I,K) 
G2  =  R 
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0150'30 

015100 

015110 

015120 

015130 

015140 

015150 

015160 

015170 

015180 

015190 

015200 

015210 

015220 

015230 

015240 

015250 

015260 

015270 

015280 

015290 

015300 

015310 

015320 

015330 

015340 

015350 

015360 

015370 

015380 

015390 

015400 

015410 

015420 

015430 

015440 

015450 

015460 

015470 

015480 

015490 

015500 

015510 

015520 

015530 

015540 

015550 

015560 

015570 

015580 

015590 

015600 
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IF  <  t32  +  1.  )  104,105,105 

015610 

104 

S2  =  -1. 

015620 

105 

SC3  =  CL0SE(G1,G2) 

015630 

IF(SMIN(I,K)  -  PRVMN)  5003,60,60 

015640 

c 

MINIMUM  STRESS  ADJUSTMENT 

015650 

5003 

CONTINUE 

015660 

PRVMN  =  SMIN(I,K) 

015670 

G1  =  SPEAK 

015680 

G2  =  SMIN<I,K)/SPEAK 

015690 

IF  (Q2  ♦  1.)  5005,5006,5006 

015700 

5005 

G2  =  -1. 

015710 

5006 

SCI  =  CL0SE(G1,G2) 

015720 

G1  =  PRVMX 

015730 

G2  =  SMIN(I,K)/PRVMX 

015740 

IF  <Q2  +  1.)  57,58,58 

015750 

57 

Q2  =  -1. 

015760 

58 

SC3  =  CLOSE (G1,G2) 

015770 

G3  =  ASTRT  -  APEAK 

015780 

Q4  *  AP  -  APEAK 

015790 

IF(G3/G4  .LT.  0.)  G3=0. 

015800 

SC2T  =  D0UN<SC1,63,64) 

015810 

IF(SC2T  -  SC2)  59,60,60 

015820 

59 

SC2  =  SC2T 

015830 

60 

SC  =  SC2 

015840 

IF(S(I,K)-SC2)  450,450,5009 

015850 

5009 

SMNGR  =  SC2 

015860 

IF  (SMIN(I,K)-SC2)  5010,5011,5011 

015870 

5010 

SMNGR  =  SC2 

015880 

soil 

CALL  BETAS(S(I,K),A1,ALP,QE) 

015890 

XKEFF=  ( S  ( I ,  K )  -SMNGR )  4:SQRT  (P I  »A  1 )  TALP 

015900 

CKTH  =  KTH*(1.-RMULT1:R) 

015910 

IF  CR  .LT.  0.  )  CKTH  =  KTH 

015920 

IF<:XKEFF*(1.-R)/(1.-CF)  .LE.  CKTH)  GO  TO  450 

015930 

5012 

G1  =  S(I,K) 

015940 

G2  =  R 

015950 

IF  (Q2  +  1.)  93,94,94 

015960 

93 

Q2  =  -1. 

015970 

94 

SC3  =  CL0SE(Q1,G2) 

015980 

IF  (  Sa,K)  -  SC3  )  95,96,96 

015990 

95 

SC3  =  S(I,K) 

016000 

96 

CONTINUE 

016010 

IF(S(I,K)  -  SPEAK)  90,5013,5013 

016020 

5013 

KLU  =  3 

016030 

c 

016040 

C 

INITIALIZATION  FOR  INTEGRATION  ROUTINE 

016050 

C  016060 


c 

MODE  =  1 

IF 

CYCLS(I,K) 

.LT. 

20  AND  SC3 

.LE.  SC 

016070 

c 

MODE  *  2 

IF 

CYCLS(I,K) 

.LT. 

20  AND  SC3 

.GT.  SC 

016080 

c 

MODE  »  3 

IF 

CYCLS<I,K) 

.GE. 

20,  SC3.GT. 

SC  AND  N.LT.NSAT 

016090 

c 

MODE  =  4 

IF 

CYCLS(I,K) 

.  GE. 

20  AND  SC3 

.LT.  SC 

016100 

c 

MODE  =  5 

IF 

CYCLSCI,K) 

.GF. 

20  AND  SC3 

*  SC 

016110 

c 

INTEGRATION 

PERFORMED  FOR 

MODE 

=  4  AND  5 

016120 

o  o  o  o  o 


90 

MODE  =  5 

016130 

IF  (  CYCLS(I,K)  -  20.  )  201,203,203 

016140 

201 

MODE  =  1 

016150 

IF  (  SC  -  SC3  )  202,205,205 

016160 

202 

MODE  =  2 

016170 

GO  TO  205 

016180 

203 

IF  f.  SC  .EQ.  SC3  )  GO  TO  205 

016190 

MODE  =  3 

016200 

IF  (  SC3  -  SC  )  204,205,205 

016210 

204 

MODE  =  4 

016220 

205 

CONTINUE 

016230 

016240 

016250 

DCYC  =  CYCLS(I,K)/30000. 

016260 

IF  (  DCYC  .LT.  1.0  )  DCYC  =  1.0 

016270 

NNCYC  =  CYCLS<I,K)/DCYC 

016280 

016290 

START  CYCLE  -  BY  -  CYCLE  ANALYSIS 

016300 

016310 

DO  310  J  =  1, NNCYC 

016320 

IF  (  MODE  -  3  )  210,210,6000 

016330 

210 

IF  (SC3  -  SC)  92,92,111 

016340 

92 

GO  TO  (101,160,101,101),KLU 

016350 

101 

IF  (OMQAl)  5014,5015,5015 

016360 

5014 

OMGAl  =  0, 

016370 

5015 

IF  (0M6A1/0MGA2  -  1.)  5017,5017,5016 

016380 

5016 

OMGAl  =  0MGA2 

016390 

5017 

IF  (0M6A1/0MQA2  .LT.  0.)  OMGAl  =0. 

016400 

SC  =  D0WN(SC2,0MGA1,0MGA2) 

016410 

111 

IF(SM1N(I,K)  -  SC)  5023,160,160 

016420 

5022 

SMNGR  =  SC 

016430 

160 

CONTINUE 

015440 

GO  TO  4050 

016450 

4051 

CONTINUE 

016460 

I SUM  =  ISUM  +  DCYC 

016470 

CYC  =  CYC  +  DCYC 

016480 

CYSUM  =  CYSUM  -  DCYC 

016490 

DN  =  CYSUM 

016500 

GO  TO  (138, 501, 501, 501), ITEM 

016510 

138 

CONTINUE 

016520 

IF  <  MODE  -  2  )  143,144,144 

016530 

143 

OMQAl  =  A1  -  ASTRT 

016540 

IF  (A1  -  AP)  300,141,141 

016550 

141 

KLU  =  2 

016560 

SC  *  SC3 

016570 

AP  =  A1 

016580 

GO  TO  300 

016590 

144 

sen  =  SCONE  (SC3) 

016600 

NPR  =  NPREV CSC, SC3, sell) 

016610 

IF  (  NPR  )  145,146,146 

016620 

145 

NPR  =  0. 

016630 

146 

NOL  =  NPR  ♦  DCYC 

016640 

o  o  o 


SC  =  sen  +  csc3-scii)*i:nol-i 

.)/(NSAT-l.) 

016650 

IF  (  SC3  -  sc  )  147,147,300 

0 1 6660 

147 

SC  =  SC3 

016670 

IF  (  MODE  -  2  )  300,300,148 

016680 

148 

MODE  =  5 

016690 

300 

CONTINUE 

016700 

310 

CONTINUE 

016710 

311 

CONTINUE 

016720 

PRVMX  =  S(I,K) 

016730 

SC2  =  SC 

016740 

80  TO  (501, 400, 400, 400), KLU 

016750 

400 

CALL  BETAS(S(I,K),A1,ALP,QE) 

016760 

XK  =  S(I,K)1:SQRT<PI*A1)*ALP 

016770 

AP  =  A1  ♦  RY(XK,PLSTRN) 

016780 

APEAK  =  A1 

016790 

SPEAK  =  S(I,K) 

016800 

PRVMN  =  SMIN(I,K) 

016810 

IF  (SC3  -  SC)  420,420,430 

016820 

420 

SCI  =  SC3 

016830 

GO  TO  501 

016840 

430 

SCI  =  SC 

016850 

GO  TO  501 

016860 

450 

ISUM  =  ISUM  +  CYCLS(I,K) 

016870 

CYC  =  CYC  +  CYCLS<I,K) 

016880 

DN  =  0. 

016890 

501 

CONTINUE 

016900 

016910 

GO  TO (500, 600, 620, 503), ITEM 

016920 

503 

I STOP  =  2 

016930 

A1  =  ATRANS/(2.4!RATI0) 

016940 

WRITE(6, 1)  A1,CYC 

016950 

1 

FORMATdHO,  70(1H*)/  5X,55HTRANSITI0N  TO  A  THRU 

CRACK 

OF  EFFECTIVE016960 

1  LENGTH,  AEFF  =  ,F9.5,4H  AT  ,F12 

.2,7H  CYCLES  / 

,1X,  70dH*)) 

016970 

GO  TO  50 

016980 

016990 

017000 

END  OF  CYCLE  BY  CYCLE  ROUTINE 

017010 

017020 

500 

CONTINUE 

017030 

GO  TO  50 

017040 

A1  EXCEEDED  AF  (AMAX) 

017050 

600 

ISTOP  =  1 

017060 

GO  TO  50 

017070 

620 

ISTOP  =  1 

017080 

WRITE(6,7) 

017090 

7 

FORMATdHO,  70dH*)/  5X,46HKMAX  APPLIED  EXCEEDS 

KSUBQ, 

PROBLEM 

TER017100 

IMINATED/IX,  70dH*)/lH0,  1X,26HLAST  CALCULATED  VALUES 

ARE///) 

017110 

DELTAK  =  XK*d.-R) 

017120 

WRITE (6, 8)  J1 , J2, I , IFLT, K, CYC, A1 , XK, XKEFF, DELTAK, DADN 

017130 

8 

FORMAT!  5X,18HBL0CK  IN  SPECTRUM  , 

14/ 

017140 

/  5X, 18HSEGMENT  NUMBER  , 

14/ 

017150 

/  SX, 18HMISSI0N  NUMBER  , 

14/ 

017160 

I  lf  i_iq  I  I.  ^  M .  1.^  L'  ll  j.^ 


/ 

5X, 18HFLIGHT  NUMBER  ,14/ 

017170 

/ 

5X,18HLAYER  IN  MISSION  ,14/ 

017180 

5 

5X, ISHACOJMULATED  CYCLES, E16. 8/ 

017190 

& 

5X,18HCRACK  LENGTH  ,E16.8/ 

017200 

7 

5X,18HKMAX  APPLIED  ,E16.8/ 

017210 

8 

5X,18HKMAX  EFFECTIVE  ,E16.8/ 

017220 

9 

5X,18HDELTA  K  ,E16.8/ 

017230 

/ 

5X,18HDA/DN  ,£16.8) 

017240 

017250 

017260 

50 

CONTINUE 

017270 

XKA  =  XK 

017280 

RETURN 

017290 

017300 

017310 

GROWTH  CALCULATIONS 

017320 

4050 

CONTINUE 

017330 

CHECK  FOR  XK  (KMAXA)  .GE.  AKC  (KSUBQ) 

017340 

CALL  BETAS(S(I,K),A1,ALP,QE) 

017350 

XK  =  S(I,K)  *SQRT(PI*A1)*ALP 

017360 

IF(XK  -  AKC)  4062,4060,4060 

017370 

4060 

ITEM  =  3 

017380 

DADN  =  0. 

017390 

GO  TO  (502, 502, 502, 4062, 4062), MODE 

017400 

4062 

XKEFF  =  (S(I,K)  -  SMNQR)*SQRT(PI*Al)*ALP 

017410 

XKEFFN  =  XKEFF/<1.-CF) 

017420 

IF  (  FGN  .EQ.  1  :i  GO  TO  9000 

017430 

IF  (  EQN  .EQ.  4  )  GO  TO  9010 

017440 

IF  (  EQN  -  2  )  8000,8000,8010 

017450 

8000 

CONTINUE 

017460 

BI-LINEAR  PARIS  EQUATION 

017470 

C  =  Cl 

017480 

AN  =  SNl 

017490 

IF  (  XKEFFN  .GE.  DKCOM  )  C  =  C2 

017500 

IF  (  XKEFFN  .GE.  DKCOM  )  AN  =  SN2 

017510 

DADN  =  C*XKEFFN*»AN 

017520 

GO  TO  3020 

017530 

017540 

8010 

CONTINUE 

017550 

TABULAR  RATE  VALUES 

017560 

DADN  =  TBLKUP(CARRAY,SNARAY,NDADN, 100, XKEFFN) 

017570 

DADN  =  10.4;:«DADN 

017580 

8020 

CONTINUE 

0175’90 

IF(RETARD.NE.O)DADNPR  =  DADN 

017600 

A1  =  A1  DADN»DCYC 

017610 

IF  (A1  -  AF)4064, 4063, 4063 

017620 

4063 

ITEM  =  2 

017630 

GO  TO  502 

017640 

4064 

IF  (  ISURF  .EQ.  0  )  GO  TO  502 

017650 

ATRANS  =  THICK  -  ( ( (XK/SIGMAY)*»2.  )/(2.l:PI) ) 

017660 

IF  (  A1  .LT.  ATRANS)  GO  TO  502 

017670 

ITEM  =  4 

017680 

123  - 


502  GO  TO  (4051,6060), I GROW 

017690 

017700 

017710 

END  OF  CRACK  GROWTH  CALCULATIONS 

017720 

017730 

017740 

017750 

INTEGRATION  ROUTINE 

017760 

6000  CONTINUE 

017770 

IF  (  MODE  -  4  )  6010,6010,6020 

017780 

6010  Q(l)  =  A1 

017790 

Q(2)  =  AP 

017800 

Z2  =  At  -  ASTRT 

017810 

Z3  =  AP  -  ASTRT 

017820- 

IF  <  Z2/Z3  .LT.  0. )  Z2  =0. 

017830 

QQ(1)  =  D0UN(SC2,Z2,Z3) 

017840 

QQ(2)  =  SC3 

017850 

IF  (  ISURF  .EQ.  0  )  GO  TO  6030 

017860 

IF  (  QC2)  .LT.  THICK  )  GO  TO  6030 

017870 

Q(2)  =  THICK 

017880 

Z2  =  THICK  -  ASTRT 

017890 

IF  (Z2/Z3  .LT.  0.)  Z2  =0. 

017900 

□Q(2)  =  D0WN<SC2,Z2,Z3) 

017910 

I STOP  =  2 

017920 

GO  TO  6030 

017930 

6020  Qd)  =  A1 

017940 

Q(2)  =  1.104:A1 

017950 

QQ<1)  =  SC3 

017960 

QQ(2)  =  SC3 

017970 

IF  (  ISURF  .EQ.  0  )  GO  TO  6030 

017980 

IF  <  Q(2)  .LT.  THICK  )  GO  TO  6030 

017990 

Q(2)  =  THICK 

018000 

ISTOP  =  2 

018010 

6030  CONTINUE 

018020 

SC  =  QQ(1) 

018030 

DO  6100  KK  =  1,2 

013040 

A1  =  Q(KK) 

018050 

SMNGR  =  SMIN<I,K) 

018060 

IF  <  SMIN(I,K)  -  QQCKK)  )  6045,6050,6050 

018070 

6045  SMNGR  =  QQCKK) 

018080 

£050  I GROW  =  2 

018090 

GO  TO  4050 

018100 

6060  GO  TO  (6070,6065,6065,6065), ITEM 

018110 

6065  GO  TO  (501, 6068 ),KK 

018120 

6068  ITEM  =  1 

018130 

ISTOP  =  0 

018140 

6070  KKK  =  KK  ♦  2 

018150 

KKKK  =  KK  +  4 

018160 

Q(KKK)  =  XKEFF 

018170 

QCKKKK)  =  DADN 

018130 

6100  CONTINUE 

018190 

Q2  =  (Q(3)-Q(4))/(Q(l)-Q(2)) 

018200 

01  =  0(3)  -  0240(1) 

018210 

04  =  (AL0G(Q(5)/0(6)))/(ALQG<Q(3)/0(4))) 

018220 

03  =  Q(5)/i:Q(3):4*04) 

018230 

05  =  1.  -  04 

018240 

DELTN  =  (0(4)1:4:05  -  0(3)1:4:05) /•:Q2:4:Q3*05) 

018250 

IF  (  DELTN  -  CYSUM  )  6200,6150,6250 

018260 

6150 

CONTINUE 

018270 

6200 

A1  =  0(2) 

018280 

IF  (  AP  -  A1  )  6210,6210,6220 

018290 

6210 

AP  =  A1 

018300 

KLU  =  2 

018310 

MODE  =  5 

018320 

6220 

I SUM  =  rSUM  +  DELTN 

018330 

CYC  =  CYC  +  DELTN 

018340 

CYSUM  =  CYSUM  -  DELTN 

018350 

DN  =  CYSUM 

018360 

IF  <  CYSUM  -  .01  )  6260,6240,6240 

018370 

6240 

IF  (  MODE  -  4  )  6010,6010,6020 

018380 

6250 

0(4)  =  (0(3)4:1:Q5+CYSUM1tQ2l:Q34:05)1:4:(l./05) 

018390 

0(2)  =  (0(4)  -Ql)/02 

018400 

DELTN  =  CYSUM 

018410 

I STOP  =  0 

018420 

GO  TO  6150 

018430 

6260 

SC  =  SC3 

018440 

6270 

22  =  A1  -  ASTRT 

018450 

23  =  AP  -  ASTRT 

018460 

IF  <23  .LE.  0.  )  23  =  22 

018470 

IF  <22/23  .LT.  0.)  22  =0. 

018480 

SC  =  DaWN<SC2,22,23) 

018490 

6280 

GO  TO  311 

018500 

018510 

018520 

END  OF  INTEGRATION  ROUTINE 

018530 

018540 

018550 

DIAGNOSTICS  FOR  IMPROPER  EON 

018560 

9000 

WRITE<6,9005) 

018570 

9005 

FORMAT <1H0,  1 X , 44HCL0SURE  MODEL  CAN  NOT  ACCEPT 

FORMAN 

EQUATION  / 

018580 

1  40X,20HEXECUTI0N  SUPPRESSED  ) 

018590 

GO  TO  9020 

018600 

9010 

MRITE<6,9015) 

018610 

9015 

FORMATdHO,  1 X , 44HCL0SURE  MODEL  CAN  NOT  ACCEPT 

WALKER 

EQUATION  / 

018620 

1  40X,20HEXECUTI0N  SUPPRESSED  ) 

018630 

9020 

ISTOP  =  1 

018640 

XKA  =  XK 

018650 

RETURN 

018660 

018670 

018680 

END 

018690 

SUBROUTINE  RUNKUT<X, Y, DX, F) 

018700 

018710 

EXTERNAL  F 

018720 

rr 

% 


% 


% 


% 


« 


L 


COMMON/STEPS/ ISEG, J1 , J2, J3, J4, J5, ISTOP, NORETRD 
10  X0=X 
X=X+DX 
H=DX 

20  IF(ABS(H).i3T.ABS(X-X0))  H=X-XO 
30  Y0=Y 
HT=H 
XT=X0 

RMAXP=1.E37 
40  YT=YO 

ASSIGN  50  TO  K 
GO  TO  100 
50  CONTINUE 
60  YP=Y 
70  HT=0.5*H 

ASSIGN  80  TO  K 
GO  TO  100 
80  CONTINUE 
90  YT=Y 
XT=XO+HT 
ASSIGN  150  TO  K 
100  CALL  F(XT,YT,PO) 

IF( ISTOP. EQ.O)  GO  TO  110 
X  =  XT 
Y  =  YT 
RETURN 

no  Y=YT4.0.5»:HT»:P0 

CALL  F<:XT+0.5*HT,Y,P1) 

IF( ISTOP. EQ.O)  GO  TO  120 

X  =  XT+O.SiHT 

RETURN 

120  Y=YT+0.5CHT1:P1 

CALL  F(XT<-0.5l:HT,Y,P2) 

IF C ISTOP. EO.O)  GO  TO  130 

X  =  XT+0.5*HT 

RETURN 

130  Y=YT-t-HT*P2 

CALL  FaT+HT,Y,P3) 

IF( ISTOP. EQ.O)  GO  TO  140 

X  =  XT+HT 

RETURN 

140  Y=YT-t-HT*(P0+2.1:(Pl+P2)+P3)/6. 

GO  TO  K, (50,80, 150) 

150  RI1AX=0. 

160  RMAX=AMAX 1 (RMAX , 0. 07*ABS( (Y-YP) /Y) ) 

IFi:(RMAX,QT.l,E-07).AND.  (RMAX.LT.RMAXP))  GO  TO  170 
XO=XO+H 

IFCXO.EQ.X)  RETURN 

IF((RMAX.LT. l.E-08).  OR. (RMAX.GT.RMAXP) )  H=H+H 
GO  TO  20 
170  H=HT 
XT=XO 


018730 
018740 
018750 
018760 
018770 
018780 
018790 
018800 
018810 
018820 
018830 
018840 
018850 
018860 
018870 
018880 
018890 
018900 
018910 
018920 
018930 
018940 
018950 
018960 
018970 
018980 
018990 
019000 
019010 
019020 
019030 
019040 
019050 
019060 
019070 
019080 
0 I 9090 
019100 
019110 
019120 
019130 
019140 
019150 
019160 
019170 
019180 
019190 
019200 
019210 
019220 
019230 
019240 


126  - 


180  YP=YT 
190  YT=Y0 

RMAXP*RMAX 
GQ  TG  70 
END 


Sawolt  Input 


TITLE 

I 

L0W-MED-HI6H  SPECIMEN  84-502,  EE3,  W»«EL£R  M=6.0,  1  CYCLE=20  SEC 

EQUATION 

SIGMOID 

MATERIAL 


INCO  718 

MSE  FIT 

DA/DT  TIMES  20 

21.0 

272.73 

-5.6942677-1.1 

1.8 

300.0 

120.0 

THRESHOLD 

21.0 

1.0 

LIMITS 

.4087 

1.1330 

0.  0.0 

ANALYSIS 

RETARD 

1.0 

0.0  6.0 

1.0 

BETA 

9.0 

1.5736  .394 

END 

LOADS 

1  0  PROOF  TEST  SPECTRUM 

MAX-MIN  SUSTAINED  LOAD  EQUIVALENT  OF  1  CYCLE  =  20  SEC  R«.4 
2.740  1.096  1. 

END 

MAX-MIN  OVERLOAD  OF  20  PERCENT  EQUIVALENT  CYCLE 
3.288  0.0  1. 


END 

END  LOADS 
SPECTRUM 
7  0 

72  1 

1  2 

1573  1 
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It  was  concluded  that  the  system  had  to  be  recalibrated 
after  each  overload  due  to  a  sudden  advancement  in  crack 
length.  — 

The  retardation  models  were  found  to  require  empirical 
parameters  that  depend  upon  the  stress  intensity  level  for 
each  overload  application.  Using  relationships  developed 
for  these  parameters,  the  CRACKS  program  using  the  Wheeler 
model  was  found  to  be  capable  of  predicting  the 
t  ime-to-fai lure  for  sustained-loading  with  periodic 
overloads  within  20  percent  of  test  data.  The  Willenborg 
model  was  found  to  be  inapplicable  to  this  problem  because 
it  depends  solely  on  stress  ratio  which  has  no  physical 
meaning  for  susta i ned-1 oad ing .  The  Wheeler  model,  on  the 
other  hand,  could  generally  be  applied  to  sustained-load 
crack  growth  using  equivalent  fatigue  cycles.  In 
conjunction  with  the  CRACKS  computer  program,  this  could 
provide  a  powerful  new  method  for  evaluating  crack  growth 
under  general  engine  mission  spectra  including  the  effects 
of  overloads . 
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